format long
xi=2;
tolerance=1e-5;
err=1;
while (err>tolerance)
fxi=xi^3+4*xi^2-10;
dfdxi=3*xi^2+8*xi;

xr=xi-fxi/dfdxi;
err=abs((xr-xi)/xr)*100;
xi=xr
end
root=xr