- /
-
A Newton's Method Fractal Code
on 6 Oct 2021
- 1
- 24
- 4
- 0
- 274
N=100;
h=.001;
[xs,ys]=meshgrid(linspace(-1,1,800),linspace(-1,1,800));
A=xs(:)+i*ys(:);
ss=1:numel(xs);
ni=N*ones(numel(xs),1);
for i=1:N
k=A(ss);
A(ss)=k-(k.^2+1)./5.*(k).^3;
d=(k-A(ss));
cc=abs(d)<h;
prob=isnan(d);
ni(ss(cc))=i;
ni(ss(prob))=N+1;
ss(cc|prob)=[];
end
ni=reshape(ni,size(xs));
imagesc(ni);