- /
-
Julia Set
on 5 Oct 2021
- 4
- 35
- 1
- 0
- 271
Bound = 1.5;
C=0.285 + 0.01i;
N = 100;
temp = linspace(-Bound, Bound, 500);
R = repmat(temp,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(temp,temp,atan(0.1*M)); axis xy; drawnow;
colormap jet
end