clear all, close all, clc
xy = [1  3]';   % Startgissningar där xy = [x y]'

dxy = 1;   % Sätt värde för att komma in i while-loopen
while  max(abs(dxy)) > 1e-10
    x = xy(1);   y = xy(2);   % Omskrivning för enklare uttryck
    
    f1 = x^2*y + x^3*y^2 - 5;
    f2 = x + 3*x*y^2 - 13;
    f3 = x + y - 3;
    f  = [f1  f2  f3]';   % Funktionsvektor. OBS f(x, y) = 0
    
    J11 = 2*x*y + 3*x^2*y^2;
    J12 = x^2 + 2*x^3*y;
    J = [J11  J12 ; 1+3*y^2  6*x*y ; 1  1];
    
    dxy = J\f   % Löser J*dxy=f där dxy=[dx dy]'
    xy  = xy - dxy;
end
x = xy(1)
y = xy(2)