r=1;n=37;
c=zeros(n);
s=zeros(n);
M=moviein(n);
for alfa=1:n
        c(alfa)=r*cos(alfa*10*pi/180);
        s(alfa)=r*sin(alfa*10*pi/180);
        if alfa>1
        plot(c(1:alfa),    s(1:alfa),    'b-', ...
             c(alfa-1:alfa), s(alfa-1:alfa), 'r-' );
        title('movimiento');
        axis([-1.5 1.5 -1.5 1.5]);
        else
        plot(c(1:alfa),      s(1:alfa),      'b-');
        title('movimiento2');
        axis([-1.5 1.5 -1.5 1.5]);
        end
        M(:,alfa)=getframe;
end

