Proving Quadratic Convergence via Taylor Expansion

Scootertaj
Messages
97
Reaction score
0

Homework Statement



The following is a modification of Newton's method:

xn+1 = xn - f(xn) / g(xn) where g(xn) = (f(xn + f(xn)) - f(xn)) / f(xn)

Homework Equations



We are supposed to use the following method:

let En = xn + p where p = root → xn = p + En

Moreover, f(xn) = f(p + En) = f(p) + f'(p)En + O(En2) = f'(p)En + O(En2) since f(p) = 0 at root

The Attempt at a Solution



Change the equation to xn+1 - xn = - f(xn) / g(xn) where g(xn) = (f(xn + f(xn)) - f(xn)) / f(xn)
Then, LHS: En+1 - En
RHS broken up into parts:
f(xn) = f(p + En) = f'(p)En + O(En2)
f(xn + f(xn)) = f(p + En + f(p + En)) = f(p + En + f'(p)En + O(En2)) = f(p) + f'(p)(En + Enf'(p) + O(En2)) + O(En2) = f'(p)(En + Enf'(p) + O(En2)) + O(En2)

So, f(xn + f(xn)) - f(xn) = [ f'(p)(En + Enf'(p) + O(En2)) + O(En2) ] - f'(p)En + O(En2)

= [ f'(p)(Enf'(p) + O(En2)) + O(En2) ] - O(En2)

Also, f(xn)2 = (f'(p)En + O(En2))2

So, En+1 - En = (f'(p)En + O(En2))2 / ( [ f'(p)(Enf'(p) + O(En2)) + O(En2) ] - O(En2) )

So, from this somehow I'm supposed to get En+1 = O(En2)
 
Physics news on Phys.org
Just a bump
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top