% Métodos Numéricos - Carlos Balsa 2008/2009
% Cap. 4 Sistemas não-lineares
% Resolução do exercicio 2 da ficha prática nº6.
% Metodo de Newton para sistema de equaçoes nao lineares
%

clear all

% alineas a) e b)

x0=[1 1]';
maxit=20;
tol=1e-6;
fun=@(x,y) [x(1)^2+x(2)^2-1; x(1)+x(1)*x(2)-1]; % função
jfun=@(x,y) [2*x(1) 2*x(2); 1+x(2) x(1)];       % jacobiana

[x,res,nits]=nle_newtsys(fun,jfun,x0,maxit,tol)

% alinea c)

x0=[1 -1]';
[y,res,nits]=nle_newtsys(fun,jfun,x0,maxit,tol)




