clear all
syms x y z
f1=x*y*z-x^2+y^2-1.34
f2=x*y-z^2-0.09
f3=exp(x)-exp(y)+z-0.41
df1x=diff(f1,x);
df1y=diff(f1,y);
df1z=diff(f1,z);
df2x=diff(f2,x);
df2y=diff(f2,y);
df2z=diff(f2,z);
df3x=diff(f3,x);
df3y=diff(f3,y);
df3z=diff(f3,z);
J=[df1x df1y df1z;df2x df2y df2z;df3x df3y df3z];
B=[f1;f2;f3];
sol=[x;y;z]-inv(J)*B;
x=1;
y=1;
z=1;
tolerancia=0.00001;
error= 20;
disp('------------x------- y--------z')
while( tolerancia
y1=eval(sol(2));
z1=eval(sol(3));
error= sqrt(power(x1-x,2)+power(y1-y,2)+power(z1-z,2));
fprintf('%15.5f %15.5f %15.5f %15.5f \n',x1,y1, z1,error)
x=x1;
y=y1;
z=z1;
end