• Remix
  • Share
  • New Entry

  • Philipp

  • /
  • Exploring the Complex Plane

on 8 Oct 2021
  • 8
  • 58
  • 1
  • 0
  • 274
clear all
m=1e3
m = 1000
[r,i]=meshgrid(-linspace(0.34224,0.343853,1e3),-linspace(0.638614,0.640529,1e3));
C=complex(r,i);
Z=C;
p=1:m*m;
d(p)=m+1;
for i = 1 : m
Z(p)=Z(p).*Z(p)+C(p);
U(p)=abs(Z(p));
a=p(find(2.0<U(p)));
d(a)=i;
p=p(find(U(p)<=2.0));
end
imagesc(log(reshape(d,m,m)))
c=colormap('jet');
axis off
Remix Tree