Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: Newton's method.
Replies: 1   Last Post: Nov 12, 1996 1:22 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Gandalf

Posts: 1
Registered: 12/7/04
Newton's method.
Posted: Nov 9, 1996 1:55 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

Can someone provide a little help on this problem?:

Suppose I wanted to write a sequence to perform the Newton inversion of
a polynomial which inverts a series 1 + a(1)x + a(2)x^2+... through a
particular degree. The sequence would look like the following:

(* Newton's inversion of a polynomial *)
guess = 1;

precision = 1;

While[precision < degree,
precision *=2;
If[precision > degree, precision = degree];
temp = iterate[poly, precision - 1];
guess = guess + guess*(1 - temp*guess); (*The Newton iteration*)
guess = iterate[guess, precision -1];
];

(* Now guess is 1 / poly through degree of 'degree' *)

----------------
iterate would be a function, iterate[f, n] that returns a polynomial f
taken through degree n. This sequence effectively doubles the precision
each pass through it. I'm a little stumped on how to implement the
iterate funtion. Any help greatly appreciated. Thanks.






Date Subject Author
11/9/96
Read Newton's method.
Gandalf
11/12/96
Read Re: Newton's method.
Daniel Lichtblau

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.