x=-10:0.01:10;
y=f(x);
plot(x,y)
yd=fd(x);
plot(x,y,x,yd)
x=10
x=x-f(x)/fd(x)
x=x-f(x)/fd(x)
x=x-f(x)/fd(x)
x=x-f(x)/fd(x)
x=x-f(x)/fd(x)
newton(10,5)
newton(10,5)
newton(10,10)
format long
newton(10,10)
newton2(10,10^(-5))
newton2(10,10^(-5))
newton2(10,10^(-5))
newton2(10,10^(-30))
2023024^(1/101)
newton(2,10) 
newton(2,10)
newton(2,10) 
newton(2,10) 
newton(2,10)
2023024^(1/101)
newton(1,10)
newton(1.5,10)
newton(1.5,100)