xxxxxxxxxx
var('t,i'); tmin,tmax,n=-pi,pi,12
f1(t,i)=(randint(5,10)+.9*cos(randint(11,20)*t+pi*i/n))
f2(t,i)=(1+.1*cos(randint(25,100)*t+pi*i/n))
f3(t,i)=(1+.01*cos(randint(200,300)*t+pi*i/n))
f(t,i)=f1(t,i)*f2(t,i)*f3(t,i)*(1+sin(t+pi*i/n))
sum([polar_plot(f(t,i),(tmin,tmax),
thickness=.5,color=hue(i/n))
for i in [1..2*n]]).show(axes=False,figsize=6)
No comments:
Post a Comment