clear
cas_spusteni = clock();
a1=1;
a2=2;
T1=2;
T2=5;
faze=45/180*pi;
T=lcm(T1,T2);
n=250;
t=linspace(0,T,n);
x=a1*sin(2*pi/T1*t);
y=a2*sin(2*pi/T2*t+faze);
clearplot
axis([-a1, a1, -a2, a2])
gset nokey
gset grid
hold on
for i=2:n-1
plot(x(i:i+1),y(i:i+1),"-b");
plot(x(i-1:i),y(i-1:i),"1");
pause(0.05)
endfor
hold off
celkovy_cas = etime(clock(), cas_spusteni);
disp(['Program bezel ',num2str(celkovy_cas),' sekund.']);