xxxxxxxxxx
var('t,i,a,b,c,d'); tmin,tmax,n=-pi,pi,4
f1(t,i,a,b,c,d)=(a+.9*cos(b*t+pi*i/n))*(1+.1*cos(c*t+pi*i/n))
f2(t,i,a,b,c,d)=(1+.05*cos(d*t+pi*i/n))*(1+sin(t+pi*i/n))
[a,b,c,d]=[randint(8,11),randint(12,24),
randint(25,81),randint(216,256)]
sum([polar_plot(f1(t,i,a,b,c,d)*f2(t,i,a,b,c,d),(tmin,tmax),
thickness=.3,axes=False,color=hue(i/(3*n)))
for i in [1..6*n]])
No comments:
Post a Comment