code

> with(plottools);

> par := t -> [t,t^2]:

> rc := t -> (1+4*t^2)^(3/2)/(-2):

> nr := t-> 1/sqrt(1+4*t^2)*[2*t,-1]:

> ev := t-> rc(t)*nr(t)+par(t):

> parabola:=plot(x^2,x=-2..2,color=black):

> evolute :=plot([ ev(t)[1],ev(t)[2],t=-1..1],color=red):

> background :=plots[display]([parabola,evolute],scaling=constrained):

> pl3 := t -> plottools[circle](ev(t),rc(t),color=blue):

> pl4:= t->plot([[ev(t)[1],ev(t)[2]], par(t)],color=magenta,thickness=2 ):

> plots[display]([seq(plots[display](
[ background,pl3(i/20.),pl4(i/20.)]),
i=-10..10)],insequence=true,scaling=constrained);