- /
-
Julia set - circular shape
on 20 Oct 2021
- 37
- 41
- 0
- 0
- 252
n=5000;
[R,I]=meshgrid(linspace(-2,2,n));
c=-1.3+0.005i;
It=zeros(n);
m=colormap([hot; 1 1 0]);
for ind=1:n*n
z=R(ind)+I(ind)*i;
it=0;
while abs(z)<2 && it<=257 && ~isinf(abs(z))
z=cosh(z^4)./(pi+log(abs(z)))-0.59;
it=it+1;
end
It(ind)=it;
end
imagesc(uint8(It))
axis equal off