• Remix
  • Share
  • New Entry

on 23 Oct 2021
  • 6
  • 33
  • 0
  • 0
  • 279
n=285;
s=zeros(1,n);
z=-1;
figure('Renderer','painters','InvertHardcopy','off');
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',[.3 .3 .3]);
end
Image
Remix Tree