%Lösning av icke-linjärt ekvationssystem
%med Newtons metod
%2x1+x2-x1x2/2-2=0, x1+2x2-cos(x2)/2-1.5=0
x=[1 0]';h=1;
while h>0.5e-4
    f=[2*x(1)+x(2)-x(1)*x(2)/2-2;
       x(1)+2*x(2)-cos(x(2))/2-1.5];
    J=[2-x(2)/2  1-x(1)/2;
           1   2+sin(x(2))/2]; 
    delx=-J\f;
    x=x+delx
    h=norm(delx,inf);
end
disp(' Ekvationens rot= '),x
