A Newton's Method Fractal Code - MATLAB Mini Hack
  • Remix
  • Share
  • New Entry

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);
Image
Remix Tree
Load full remix tree
Loading...
Go to top of page