Newton Ecs No Lineales.docx

Metodo de Newton para sistemas de ecuaciones no lineales %fecha: 20/03/2013 %autor: JMA clc clear %format long x=input('Dame el vector de valores iniciales: '); EPS=1e-3; MAXIT=input('Dame el número de iteraciones: '); k=1; while k<=MAXIT F=[x(1)^2+x(2)^2-1 ; x(1)^2-x(2)^2-1]; J=[2*x(1) 2*x(2) ; 2*x(1) -2*x(2)]; %[h,error]=funcion_gauss_pivoteo(J,-F); h=J\-F; %if error %disp('Los valores iniciales no son adecuados, elige otros') %break %else fprintf('la iteración %d es: ',k); xn=x+h if norm(h)< EPS fprintf('\nSe encontró una solución en: '),xn break end x=xn; k=k+1; end %end if k>MAXIT disp('En número máximo de iteraciones excedido') end