count=0
h=input('assegna h')
np=input('assegna numero di iterate')
ym=input('minimo di y per finestra')
yp=input('massimo di y per finestra')
mu=h^2
clear x y
hold off
plot([-pi 0 pi 0],[0 0 0 pi],'*m')
axis([-pi pi ym yp])
hold on
for iii=1:30
[x0 y0]=ginput(1)

x(1)=x0;
y(1)=y0;

plot(x0,y0,'og')

for k=2:np
  y(k)=y(k-1)-mu*sin(x(k-1));
  x(k)=x(k-1)+y(k);
% riduzione al toro; qui -pi < x < pi ; 0 , y < 2*pi
  x(k)=mod(x(k),2*pi);
  if(x(k)>pi)
     x(k)=x(k)-2*pi;
  end
  y(k)=mod(y(k),2*pi);
%  if (y(k)>pi)
%     y(k)=y(k)-2*pi;
%  end
end
count=count+1
cc=mod(count,6)
if cc == 1
  ccol='.r'
elseif cc==2
  ccol='.b'
elseif cc==3
  ccol='.y'
elseif cc==4
  ccol='.c'
elseif cc==5
  ccol='.g'
elseif cc==0
  ccol='.k'
end
plot(x,y,ccol)
end





