function [x, err, it] = my_newt_rec(a, x0, tol, maxit)
%metodo delle tangenti per il calcolo del reciproco
%condizione di arresto: abs(a*x-1)<tol o it>=maxit
n=length(a);
err=norm(a*x0-eye(n));
it=0;
while(err>=tol & it<maxit)
    x=2*x0-x0*a*x0;
    it=it+1;
    x0=x;
    err=norm(a*x0-eye(n))
end
end