Matlab - Iteraciones con ciclo

clc;

disp('METODO DE NEWTON PARA SISTEMAS DE ECUACIONES NO LINEALES');

syms x y f g;

funcion_x= input('Introduce la primera Función: ');

funcion_g= input('Introduce la segunda Función: ');

x0= input('Introduce x0: ');

y0= input('introduce y0: ');

%DERIVADAS

der_funcion_fx=diff(funcion_x,x);

der_funcion_fy=diff(funcion_x,y);

der_funcion_gx=diff(funcion_g,x);

der_funcion_gy=diff(funcion_g,y);

%SUSTITUYENDO VALORES

%F

f1=subs(funcion_x,x,x0);

f=subs(f1,y,y0);

%G

f2=subs(funcion_g,x,x0);

g=subs(f2,y,y0);

%SUSTITUYENDO EN DERIVADAS

%FX

f3=subs(der_funcion_fx,x,x0);

fx=subs(f3,y,y0);

%FY

f4=subs(der_funcion_fy,x,x0);

fy=subs(f4,y,y0);

%GX

f5=subs(der_funcion_gx,x,x0);

gx=subs(f5,y,y0);

%GY

f6=subs(der_funcion_gy,x,x0);

gy=subs(f6,y,y0);

%DELTA X y DELTA Y

format long

delta_x=(-f*gy+g*fy)/(fx*gy-fy*gx);

delta_y=(-g*fx+f*gx)/(fx*gy-fy*gx);

delta_x=double(delta_x);

delta_y=double(delta_y);

disp('Delta x= ');

disp(delta_x);

disp('Delta y= ');

disp(delta_y);

x1=x0+delta_x;

y1=y0+delta_y;

disp('Nuevo Valor de x0');

disp(x1);

disp('Nuevo Valor de y0');

disp(y1);

%CALCULANDO ERROR

error=sqrt((x1-x0)^2+(y1-y0)^2);

disp('Error= ');

disp(error);

>> newton_no_lineal7

METODO DE NEWTON PARA SISTEMAS DE ECUACIONES NO LINEALES

Introduce la primera Función: x^2-10*x+y^2+8

Introduce la segunda Función: x*y^2+x-10*y+8

Introduce x0: 0

introduce y0: 0

Delta x=

   0.800000000000000

Delta y=

   0.880000000000000

Nuevo Valor de x0

   0.800000000000000

Nuevo Valor de y0

   0.880000000000000

Error=

   1.189285499785481

Delta x=

   0.191787221105863

Delta y=

   0.111711737096164

Nuevo Valor de x0

   0.991787221105863

Nuevo Valor de y0

   0.991711737096164

Error=

   0.221950107872359