- /
-
Collatz Feather
on 7 Oct 2021
- 77
- 724
- 6
- 0
- 226
c=@(n)(mod(n,2)*(2*n+1)+n)/2;
s=@cumsum;
figure('Color','#02021b','Renderer','painters')
hold on
for a=3:9:20216
while a(1)>1
a=[c(a(1)) a];
end
t=s(rescale(mod(a,2),-.09,.08));
plot(s(cos(t)),s(sin(t)),'Color',[1 1 1 .1]);
end
axis equal off