444.sagews

t = var("t")
a=.5
b=1
tmini=0
tmaxi=10*pi

a,b,tmini,tmaxi

(0.500000000000000, 1, 0, 10*pi)

h1 = lambda t: [ (a+b)*cos(t)+b*cos(((b+a)/b)*t), (a+b)*sin(t)-b*sin(((b+a)/b)*t) ]
q1 = parametric_plot(h1(t), tmini, tmaxi, rgbcolor=(1,0,0))
show(q1)