||
脚本从MATLAB相应程序移植而来,用到了一些复数数组的计算。
Z=+0.04+.020211004j col=9 m=1200 cx=-1 cy=-0.5 l=2 x=linspace(cx-l,cx+l,m) y=linspace(cy-l,cy+l,m) [X,Y]=meshgrid(x,y) C=X+1j*Y for k in range(col): Z=conj(Z)**(Z)+C W=exp(-abs(Z)) axes(aspect=1.5, axis=False) imshow(rot90(W)) axis([100, 1200, 200, 1000])
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2023-3-25 16:05
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社