- /
-
Julia Set
on 5 Oct 2021
- 3
- 25
- 0
- 0
- 251
B = 1.5;
C=-0.7269 + 0.1889i;
N = 100;
T = linspace(-B, B, 500);
R = repmat(T,500,1);
I = transpose(R);
Z = R + (1i*I);
M = repmat(N,500,500);
for n = (1 : N)
Z = Z.^2 + C;
M((abs(Z) > 2)) = n;
Z((abs(Z) > 2)) = NaN;
imagesc(T,T,atan(0.1*M)); axis xy; drawnow;
colormap jet
end