This thread stimulated me into playing with a few ideas,
I pieced together this little newton-raphson procedure...
> newton := proc( f, # the function x0 # the initial guess n, # step count limit tol # error tolerance)
local x, g, k;
g := D(f); # the derivative f’(x)
x[0] := evalf(x0); # initialize the iteration
for k from 1 to n do # loop for newton’s iteration
# Newton’s iteration formula