- /
-
Spirograph Collection
on 12 Oct 2021
- 2
- 35
- 0
- 0
- 253
clear all; close all; clc
r1 = 1;
t = linspace(0,100,1000);
i = 1;
v = linspace(0.2,0.8,4);
for r2 = v
for r3 = r2 * v
subplot(4,4,i), i = i +1;
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
end
end