% questo programma calcola e mostra in una figura
% la prima orbita della mappa standard del pendolo
mu=h^2
clear x y
x(1)=x0  % queste condizioni iniziali devono essere assegnate 
y(1)=y0  % prima di chiamare questo programma
% dare hold off  solo per la prima orbita; se si vogliono tracciare
         % diverse orbite sulla stesa figura, basta ridare nuovi
         % valori a x0,y0 e richiamare lo stesso programma
plot(x0,y0,'or') % segnare con un circoletto la condizione iniziale
hold on
for k=2:np % il numero np di iterate deve essere assegnato
  y(k)=y(k-1)-mu*sin(x(k-1));
  x(k)=x(k-1)+y(k);
% riduzione al toro, -pi < x < pi, -pi < y < 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
plot(x,y,'.g') % plot con dei punti verdi; per cambiare il colore,
               % usare plot(x,y,'.r') per rosso, ecc. 

  % per un programma che gestisce automaticamente i cambiamenti di colore 
  % e di finestra, usare stmapauto.



