脚本从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])