na-lab-newton

code

clc;

x0=input("enter initial approximation");

syms x;

g=x^2-17;

f=inline(g);

dg=diff(g);

df=inline(dg);

n=input("enter no of iterations");

e=input("Enter tolerance");

for i=1:n

x1=x0-(f(x0)/df(x0));

fprintf("The root is %f\n",x1);

if(abs(x1-x0)<e)

break;

end

x0=x1;

end