%newton stazionario
clear all
close all
clc
%in realtà ci sarebbe anche le iterazioni massime da considerare
%qui non le ho messe perchè ci mette 17893 iterazioni
x=linspace(-pi/2,pi/2);
f=@(x) tan(x);
df=@(x) 1./cos(x)^2;
tol=sqrt(eps);

[c] = NSTAZ(f,df,x,tol)
fplot(f)
xline(c)

function [x0] = NSTAZ(f,df,x,tol)
    x0=x(2);
    l=df(x0);
    k=1;
    e(1)=abs(f(x0));

    while e(k)>tol
        x0= x0 - f(x0)/l;
        k=k+1;
        e(k)=abs(f(x0));
    end
end