- /
-
Recamán's sequence
on 4 Oct 2021
- 3
- 22
- 0
- 0
- 252
hold on
a=0:1e-4:pi;
n=64;
v(1)=0;
i=1;
b=ones(999,1);
colormap jet
c=colormap;
for i=1:n
t=v(i)-i;
if t>0 & b(t)
v(i+1)=t;
b(t)=0;
h=1;
else
v(i+1)=t+2*i;
b(t+2*i)=0;
h=2;
end
plot((v(i)+(-1)^h*(i-i*cos(a))/2)*1,-(-1)^i*i*sin(a)/2,'Color',c(ceil(255*i/n),:),'LineWidth',2)
axis off
end