- /
-
Wobbly spirograph
on 5 Oct 2021
- 5
- 37
- 2
- 0
- 213
% wobbly spirograph
t = linspace(0,100,15000);
r1 = 10 + 4 * sin(t/2000);
r2 = 8 + 3 * sin(t/1000);
r3 = 2 + sin(t/1000);
a = r1 - r2;
b = (r1 - r1 ./ r2 );
x = a .* cos(t)+ r3 .* cos(t.*b);
y = a .* sin(t)+ r3 .* sin(t.*b);
plot(x, y)
axis equal off