x = [-0.7; -2.5];
for k=1:12
    f = [x(1)^2*x(2) + 1; exp(x(1)) + x(2)]
    fder = [2*x(1)*x(2) x(1)^2; exp(x(1)) 1];
    x = x - inv(fder)*f
end