%openf94, newtons metod för system
x=[1 1]';
h=[1 1]';
result=[];
while norm(h)>1e-6
    f=[2*x(1)+x(2)-x(1)*x(2)/2-2;
       x(1)+2*x(2)-(cos(x(2)))/2-1.5];
    jac=[2-x(2)/2 1-x(1)/2;
         1 2+(sin(x(2)))/2];
    h=jac\f;
    normh=norm(h);
    x=x-h;
    result=[result; [normh x']];
end
result
    
    
     