- /
-
Numerical Phoenix
on 12 Oct 2021
- 24
- 128
- 0
- 0
- 264
L=650;
v=linspace(-1,1,L);
[x,y]=meshgrid(1.5*v);
c=[.306+.02i,.36+.4i,.2+.6i];
for k=1:3
z=x+y*i;
for n=1:30
z=z.^2+c(k);
end
n=1;
for t=linspace(0,pi,L)
w=exp(t*i)*(x+y*i);
a=interp2(x,y,double(abs(z)<2),real(w),imag(w),'linear',0);
r(n,:,k)=sum(a,1)';
n=n+1;
end
end
r=r/max(r(:));
imshow(r);