Páginas

newton raphson ejercicio nro 08

clc
clear all
syms x0 y0
f1=((x0-y0)*(x0))/((2-x0-y0)*(1-x0))-2.6
f2=((2*y0^2)/((2-x0-y0)*(x0-y0)))-3.1
df1x=diff(f1,x0);
df1y=diff(f1,y0);
df2x=diff(f2,x0);
df2y=diff(f2,y0);
error =20;
A=[df1x df1y;df2x df2y];
B=[f1;f2];
sol=[x0;y0]-inv(A)*B;
tol=0.001;
x0=0.8;
y0=0.4;
fprintf(' x0= %15.5f y0 = %15.5f \n',x0,y0)


disp(' x y error ')
disp('-----------------------------------------------------')
while error>tol
x=eval(sol(1));
y=eval(sol(2));
error=(((x-x0)^2)+(y-y0)^2)^0.5;
fprintf(' %15.5f %15.5f %15.5f\n',x,y,error)
x0=x;
y0=y;
end