function [x] = mysqrt(a,x0,n)
%approssimare la radice quadrata di un numero positivo a
%con il metodo delle tangenti con punto iniziale x0
%generando n termini della successione
%input a,x0,n
%output x approssimazione di sqrt(a)
f=@(x)x^2-a;
f1=@(x)2*x;
for k=1:n
    %x(k+1)=x(k)-f(x(k))/f'(x(k));
 %   x=x0 -(x0^2-a)/(2*x0);
 x=x0-f(x0)/f1(x0);
    x0=x;
end

end