xxxxxxxxxx
def _(N=[2*10^3,5*10^3,10^4]):
L=[CC(0,0)]
c=[colormaps.cool((N^(-1)).n()*i)[:3]
for i in [1..N]]
z=lambda n:exp(I*pi*n*ln(n/4)*2^.5).n()
for n in [1..N]: L.append(L[n-1]+CC(z(n)))
p=sum([line([L[i],L[i+1]],color=c[i])
for i in [200..N-1]])
p.show(figsize=(5,5))
No comments:
Post a Comment