I haven't yet heard from Numerical Recipes about my translation of their
Brent's method routine, so as a stop gap here is a little program which
implements a fail-safe Newton's method:
load(fsnewton);
fsnewton(x-tan(x),3.5,4.5);
It seems to be pretty robust, assuming that the initial function is
continuous and differentiable over the given interval.
-Alasdair
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fsnewton.mac
Type: application/octet-stream
Size: 2351 bytes
Desc: not available
Url : http://www.math.utexas.edu/pipermail/maxima/attachments/20070516/9f1de5cf/attachment.obj