
clear;

v0=[0; 0]; 
maxit=50;
tol=10^(-10);

[v_hist, step_hist, residual_hist]=newton_sistemi(v0,maxit,tol);