xxxxxxxxxx
var('t'); n=10
f1=exp(cos(t)^2+sin(t))-3*cos(4*t)
f2=-floor(4*t)*abs(sin(4*t))/2
s1=sum([polar_plot(i*f1/n,(t,0,2*pi),
color=Color([i/n,0,1-i/n]))
for i in [1..n]])
s2=sum([polar_plot(i*f2/n,(t,0,2*pi),
color=Color([1-i/n,i/n,0]))
for i in [1..n]])
(s2+s1).show(axes=false)
No comments:
Post a Comment