Proving Quadratic Convergence via Taylor Expansion

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 1K views
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