- /
- 
        Recamán's sequence (dark-colored)
        on 24 Oct 2021
        
        
 
    - 9
- 45
- 1
- 0
- 275
n=285;
s=zeros(1,n);
z=-1;
f=winter(n);
figure('Color','k','Renderer','painters');
hold; axis off equal;
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


 

 
             
             
            