- /
-
Freezing Mandelbrot Set
on 6 Oct 2021
- 5
- 23
- 0
- 0
- 242
p = 500;
x0 = -0.6;
r = 1.5;
k = 50;
x = linspace(x0-r, x0+r, p);
y = linspace(-r, r, p);
[X, Y] = meshgrid(x, y);
Z = X + 1i*Y;
m = k*ones(size(Z));
z = Z;
for n = 0:k
z = z.*z+Z;
h = abs(z)>2;
z(h) = NaN;
m(h) = n;
end
imagesc(m)
colormap(flipud(1-hot))
axis equal
axis off