- /
-
Collatz's Feather Negative
on 17 Oct 2022
- 73
- 177
- 1
- 5
- 220
% Credits for the code go to "taoizm"
c=@(n)(mod(n,2)*(2*n+1)+n)/2;
s=@cumsum;
figure('Color','w','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',[0 0 0 .1]);
end
axis equal off