1.3 Speed

In numerical analysis, we typically know that a sequence converges, unlike Intro to Real Analysis, so instead we are concerned with the speed of convergence.

We have a sequence pn with limnpn=p.

Order of Convergence of a Sequence

pn converges with order α1 if

limn|pn+1p||pnp|α=λ

for some real number λ>0

αln|pn+2ppn+1p|ln|pn+1ppnp|.

Order and Significant Digits

Relating this to significant digits we get:

|pn+1p|λ|pnp|α|pn+1pp|λ|pnpp|α|p|α1log|pn+1pp|log|pnpp|αlog(λ|p|α1)d(pn+1)αd(pn)log(λ|p|α1).
  1. for linear convergence (α=1),d(pn+1)d(pn)logλ, so each term has a fixed number more significant digits of accuracy (approximately equal to logλ) than the previous;
  2. for quadratic convergence (α=2),d(pn+1)2d(pn)log(λ|p|), so each term has double the number of significant digits of accuracy of the previous, give or take some;
  3. for cubic convergence (α=3),d(pn+1)3d(pn)log(λ|p|2), so each term has triple the number of significant digits of accuracy of the previous, give or take some

and so on

It can be shown for α>1 that

d(pn0+k)(dn0C)αk+C

where

C=log(λ|p|α1)α1.

Rate of Convergence

The sequence pn converges to p with the rate of convergence O(bn) if

bn converges to 0 and

|pnp|λ|bn|

for some real number λ>0 and all sufficiently large n

usual comparison:

bn=1nP or 1an