Matlab - newton-raphson

clear

x0=0.1;

Tol=1.e-6;

maxiter=100;

delta=1.e-6;

epsilon=1.e-10;

for k=1:maxiter

    x1=x0-((199000-(1849.24*(1-(1+x0)^-120))*x0^-1 /(1849.24*(x0^-2)*(1-(1+x0)^-120)-(221908.8*((1+x0)^-121)*(x0^-1));

    err=abs(x1-x0);

    relerr=2*err/(abs(x1)+delta);

    x0 = x1;

    y = (199000-(1849.24*(1-(1+x0)^-120))/x0)

    if (err < delta)|(relerr<delta)|(abs(y)< epsilon), break, end

end

x0

y