Newton's method
- To: mathgroup at smc.vnet.net
- Subject: [mg26231] Newton's method
- From: "drek" <drek1976 at yahoo.com>
- Date: Wed, 6 Dec 2000 02:16:18 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
Hi all,
I am trying to use the Newton-Raphson method to find the root to an
equation.
The formula looks like this:
newton[f_, x_, x0_, n_, opts___] :=
With[{df = D[f, x]}, FixedPointList[(x - f/df) /. x -> # &, N[x0], n,
opts]]
with f=Sqrt[x^2 - 2584] * Coth[0.00128 * Sqrt[x^2 - 2584]]+ Sqrt[x^2 -
1013].
When I set
newton[f, x, 31, 20]
I end up getting values which do not converge. However, if I were to plot
the function using the Plot function in Mathematica (between x values of 32
and 51), it seems like the root is somewhere near 33.
I would thus like to know if perhaps there is something wrong with this
formulation for the Newton's method, or that there is some quirk in the Plot
function and the function, f, in fact do not have a root at all.
Thanks.
Derek