436.sagews

t = var("t")
a=4
b=1
tmini=0
tmaxi=10000*pi

a,b,tmini,tmaxi

(4, 1, 0, 10000*pi)

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

t = var("t")
a=1/3
b=1
tmini=0
tmaxi=10000*pi

a,b,tmini,tmaxi

(1/3, 1, 0, 10000*pi)

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

t = var("t")
a=4/3
b=1
tmini=0
tmaxi=10000*pi

a,b,tmini,tmaxi

(4/3, 1, 0, 10000*pi)

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

t = var("t")
a=pi
b=1
tmini=0
tmaxi=100000*pi

a,b,tmini,tmaxi

(pi, 1, 0, 100000*pi)

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