angulo=input('hasta que ángulo deseas simular( 0..360)= ');
v=input('Deseas una visualización bidimensional ó tridimensional(2,3)=');
retardo=input('Dame un número de retardo en la visualización(1..10000)=');
r=1;
n=angulo/10+1;
c=zeros(n,1); s=zeros(n,1);
uno=zeros(n,2); dos=zeros(n,2); beta=zeros(n,2);
M=moviein(n);
% dos(alfa,:),uno(alfa,:), 1:alfa,     'r-', ...
% 1:alfa,     c(1:alfa),    s(1:alfa), 'y-');
for alfa=1:n
        c(alfa,1)=r*cos((alfa-1)*10*pi/180);
        s(alfa,1)=r*sin((alfa-1)*10*pi/180);
	%beta(alfa,:)=
	uno(alfa,:) = [ c(alfa,1)-s(alfa,1) c(alfa,1)+s(alfa,1)];
	dos(alfa,:) = [-c(alfa,1)-s(alfa,1) c(alfa,1)-s(alfa,1)];
	for i=1:retardo;
	m=retardo;
	end
	if v==2
        plot(dos(1:alfa,1), uno(1:alfa,1), 'b-', ...
	     dos(1:alfa,2), uno(1:alfa,2), 'g-', ...
             dos(alfa,:), uno(alfa,:),'r-');
	axis([-1.5 1.5 -1.5 1.5]);
	title('movimiento bidimensional');
	else
	eje=1:alfa;
	plot3(dos(1:alfa,1),  eje,uno(1:alfa,1), 'g-', ...
	      dos(alfa,1),    eje,uno(alfa,1),   'r-', ...
	      dos(1:alfa,2),  eje,uno(1:alfa,2), 'b-', ...
	      dos(alfa,2),    eje,uno(alfa,2),   'r-');
        axis([-1.5 1.5 0 n+1 -1.5 1.5]);
        end
	title('movimiento tridimensional');
	M(:,alfa)=getframe;
end

