julia's bugs

Peter Stampfli
on 6 Oct 2021
- 4
- 13
- 4
- 0
- 192
s=1000;
I=ones(s,s);
w=2.2;
d=w/s;
m=200;
r=1.4;
a=-0.225+0.05*i;
n=6;
for h=1:s
for k=1:s
z=(h-s/2)*d+i*(k-s/2)*d;
c=0;
while (abs(z)<r)&(c<m)
c=c+1;
z=z^n+a/z;
end
if (c<m)&(c>0)
I(h,k)=min(2,c^0.4);
end
end
end
imshow(I,[])