• Remix
  • Share
  • New Entry

on 24 Oct 2021
  • 9
  • 32
  • 1
  • 0
  • 275
n=285;
s=zeros(1,n);
z=-1;
f=winter(n);
figure('Color','k','Renderer','painters');
hold; axis off equal;
Current plot held
for i=1:n-1
t=s(i)-i;
for j=1:i+1
if s(j)==t||t<0
t=s(i)+i;
break;
end
end
s(i+1)=t;
r=abs(s(i+1)-s(i))/2;
m=(s(i)+s(i+1))/2;
x=linspace(-r+m,r+m,1e3);
y=z^i*sqrt(r^2-(x-m).^2);
plot(x,y,'Col',f(i,:));
end
Remix Tree