function chaos2
% tato funkce animuje kmitavy pohyb bodu s linearne rostouci frekvenci
n = 99;
R = ones(n,1);
fi = linspace(0,2*pi,n)';
om = 1+nlinspace(n)';

I1 = 1:2:n;
I2 = 2:2:n;

figure
hold on;
tisk1 = plot(R(I1).*cos(fi(I1)),R(I1).*sin(fi(I1)),'marker','o','markerfacecolor','b','markersize',10,'linestyle','none');
tisk2 = plot(R(I2).*cos(fi(I2)),R(I2).*sin(fi(I2)),'marker','o','markerfacecolor','r','markersize',10,'linestyle','none');
axis equal;
axis([-2 2 -2 2]);
axis off


t = 0;
dt = 0.02;
while 1
    R = 1 + 0.5*sin(om*t);
    set(tisk1,'Xdata',R(I1).*cos(fi(I1)),'Ydata',R(I1).*sin(fi(I1)));
    set(tisk2,'Xdata',R(I2).*cos(fi(I2)),'Ydata',R(I2).*sin(fi(I2)));
    drawnow;
    t = t + dt;
end


function y = nlinspace(n)
x1 = linspace(0,1,n/2+1);
x2 = linspace(1,0,n/2);
y = [x1,x2];