%%%--------------------------------------------------------------------------%%%
%%%                                Item 1                                    %%%
%%%--------------------------------------------------------------------------%%%
clear
clc

x0   = [1; 1; 1]; % vetor de entrada x
tol  = 1e-8;      % tolerancia para o residuo (criterio de parada - metodo Newton)
nmax = 1000;      % numero maximo de iteracoes (criterio de parada - metodo Newton)

%%%--------------------------------------------------------------------------%%%
% item (b)
sol   = fsolve(@fun,x0); % calcula a solucao sol usando a funcao pronta fsolve

% item (c)
fsol  = fun(sol);        % calcula o resultado da funcao aplicada a solucao sol
r_sol = norm(fun(sol')); % calcula o residuo r_sol usando a norma

%%%--------------------------------------------------------------------------%%%
% item (d)
% Calcula a solucao x e o residuo res usando o metodo de Newton
[x res] = nonlinNewton(@fun,@Jac,x0,tol,nmax);
r_x = norm(fun(x));

%%%--------------------------------------------------------------------------%%%
printf("O resultado obtido com a funcao pronta fsolve e x = [%f %f %f], com residuo igual a %e.\n", sol(1), sol(2), sol(3), r_sol)

printf("O resultado obtido com o metodo de Newton e x = [%f %f %f], com residuo igual a %e.\n", x(1), x(2), x(3), r_x)
%%%--------------------------------------------------------------------------%%%