clear x y
x(1)=x0
y(1)=y0
per=2

mu=h^2;

%iterazioni del metodo di Newton
nit=10
for n=1:nit

    xb=x(n);yb=y(n);
    An=eye(2);
% iterazioni della mappa, e jacobiana
    for p=1:per
        A=[1-mu*cos(xb) 1; -mu*cos(xb) 1];
        An=A*An;
        yb=yb-mu*sin(xb);
        xb=xb+yb;
        xb=xb-2*pi*floor(xb/(2*pi));
        yb=yb-2*pi*floor(yb/(2*pi));
%        if yb > pi 
%          yb=yb-2*pi;
%        end
    end
% metodo di newton 
    v=(An-eye(2))^(-1)*[xb-x(n) ; yb-y(n)];
    x(n+1)=x(n)-v(1);
    y(n+1)=y(n)-v(2);
    [x(n+1) y(n+1)]
end

% risultato
[x(nit+1) y(nit+1)]
% controllo di convergenza
dif1=x(nit+1)-x(nit) 
dif2=y(nit+1)-y(nit)
% matrice jacobiana
An
% suoi autovalori ed autovettori
[V,D]=eig(An)
% se convergente usa come punto iniziale
if abs(dif1) < 1e-12 & abs(dif2) <1e-12
   x0=x(nit+1)
   y0=y(nit+1)
end









