• Remix
  • Share
  • New Entry

on 12 Oct 2021
  • 6
  • 14
  • 0
  • 0
  • 276
p = [1/24 0 -0.5 0 1];
dp=polyder(p);
rt=roots(p);
f=@(q) polyval(p,q);
g=@(q) polyval(dp,q);
n=200;
s=0.5*linspace(-1,1,n);
[X,Y]=meshgrid(s,s);
Z=X+1i*Y;
it=zeros(n);
for m=1:numel(Z)
z = Z(m);
for k=1:99
z=z-f(z)/g(z);
if ~isempty(find(abs(rt-z)<1e-6))
it(m)=k;
break
end
end
end
imagesc(it)
axis square off
Remix Tree