% Metodo de Newton para sistema de equaçoes nao lineares
% Resolução do exercicio 2 da ficha nº 7
%

x=[1 0.5]';
k=1;
max_it=20;
tol=1e-6;
erro=1;

while (erro>tol & k<=max_it)
k=k+1
F=[x(1)^2-2*sin(x(2))-1; 3*x(1)^2*x(2)-pi];
J=[2*x(1) -2*cos(x(2)); 6*x(1)*x(2) 3*x(1)^2];
d=J\(-F);
x=x+d
erro=norm(d)
end

