Julia Fractal

Ciro Bermudez on 19 Oct 2021
  • 106
  • 63
  • 9
  • 268
n = 500;
[R, I] = meshgrid(linspace(-2, 2, n));
c= -0.795 + 0.15i;
It = zeros(n, n);
m = colormap([hot; 1 1 0]);
for ind = 1: n*n
z = R(ind) + 1i*I(ind);
it = 0;
while abs(z)<2 && it<=size(m, 1) && ~isinf(abs(z))
z = z^2 + c;
it = it + 1;
end
It(ind) = it;
end
imagesc(uint8(It-1))
axis('off')

Remix tree

Load full remix tree