// activité 5

function y=h(x)
y=[x(1)^2+x(2)^2-2;x(1)^2-x(2)^2-1]
endfunction

function y=dh(x)
y=[2*x(1),2*x(2);2*x(1),-2*x(2)];
endfunction

function [x,nb,xx]=newton2(f,df,x0,eps)
err=1;
xt=x0;x=x0;xx=[x0];
nb=0;
while (err>eps)
          nb=nb+1;
          x=x-df(x)\f(x);
          err=norm(x-xt);
          xt=x;
          xx=[xx,x];
end
endfunction