//Example 5.29
//Newton Raphson Method
//Page no. 199
clc;clear;close;
deff('x=f(x)','x=x^3+2*x^2+10*x-20')
deff('x=f1(x)','x=3*x^2+4*x+10')
printf('n\txn\t\t\f(xn)\t\tf1(xn)\t\tXn+1\t\tError\n')
printf('-----------------------------------------------------------------------------------------------------\n')
x0=01;e=0.00001
for i=1:4
    x1=x0-f(x0)/f1(x0)
    e1=abs(x0-x1)
    printf(' %i\t%.10f\t%.10f\t%.10f\t%.10f\t%.10f\n',i-1,x0,f(x0),f1(x0),x1,e1)
    x0=x1;
    if abs(x0)<e then
        break;
    end
end
printf('\n\nThe solution of this equation after %i Iterations is %.10f',i,x1)