朱利亚集合，曼德博集合.zip_曼德博_曼德博集合_朱利亚集合

function Julia(c,n,N,zmax) % z(n+1)=z(n)^2 + c % c??????? % n???????iteration??? % N???????? % zmax??????????? % ????NxN??? M= zeros(N,N); % ????NxN????? T= zeros(N,N); % r???????????? r = 2; % ?x????N????? d = linspace(-zmax,zmax,N); % ??A??? A = ones(N,1)*d+i*(ones(N,1)*d)'; % Iteration for k = 1:n B = (abs(A)<=r); T = T + B; S = ((B==0) & (M==0)); M(S) = T(S); % Iteration. if ~any(M==0) break; end % plot imagesc(atan(0.1*M)); title('Julia-Mengen'); xlabel('Realteil'); ylabel('Imagin?rteil'); % animation drawnow; % recursion A = A.*A+ones(N,N).*c; end end