All published worksheets from http://sagenb.org
t = var("t") a=5 b=1 tmini=0 tmaxi=5*pi
a,b,tmini,tmaxi
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)
+