mandelbrot.m

来自「mandelbrot 分型图的仿真与绘制。mandelbrot.m, mande」· M 代码 · 共 8 行

function W=mandelbrot(X,Y,c,a,n_iter)
C=X+i*Y;W=zeros(size(X));Z=c;
for k=1:n_iter
    Z=Z.^2+C;
    i0=find(abs(Z)>a);W(i0)=k;Z(i0)=NaN;
end
i0=find(W==0);W(i0)=NaN;