function newtF(x0,y0,n)

F = inline('[x^2-3*x*y ; 5*x*y]','x','y');
dF = inline('[2*x-3*y , -3*x; 5*y , 5*x]','x','y'); % dF = derivative of F

v=[x0 ; y0]
list = v;

% Newton iteration
for k=1:n
    A = dF(v(1),v(2));
    b = F(v(1),v(2));
    z = A\b;  % A\b solves the linear equation A*z=b
    v = v-z;
    list = [list,v];
end

list'  % output