function [x] = mysqrt1(a,x0,n, f, f1) %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,f e f1 %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, a)/f1(x0); x0=x; end end