- /
-
The NeverEnding Spiral
on 4 Oct 2021
- 58
- 122
- 1
- 0
- 259
% Based on Mandelbrot Chapter of Cleve Moler "Experiments in MATLAB"
% By Victoria A. Sablina
a=-1.74975914513036646;
b=-3.68513796e-9;
w=3.2e-12;
s=5.0e-15;
x=a-w:s:a+w;
y=x'-a+b;
n=numel(x);
e=ones(n,1);
z0=x(e,:)+i*y(:,e);
z=zeros(n);
c=z;
d=2048;
for k=1:d
z=z.^2+z0;
c(abs(z)<2)=k;
end
image(c);
axis image off;
colormap(flipud(jet(d)));
title('MATLAB Central 20^{th} anniversary!!!');