julia1.m

matlab的大量实例代码

function W=julia1(X,Y,c,n_iter,a)
Z=X+i*Y; W=zeros(size(X));
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;