xxxxxxxxxx
def _(N=[20,25,..,40],M=[20,40,..,100]):
f(l,k)=[l*sin(pi*k/M)+sin(2*pi*k/M),\
-l*(1+cos(pi*k/M)+cos(2*pi*k/M))]
L1=[f(j,i) for j in [1..N] for i in range(-2*M,2*M+1)]
L2=[f(j,i) for i in range(-2*M,2*M+1) for j in [1..N]]
l1=line(L1,color='#3636ff'); l2=line(L2,color='#ff3636')
graphics_array(([l1],[l2])).show(figsize=5)
No comments:
Post a Comment