Past exam of the mathematics course of the University of Cambridge / 2014 / ib / Paper 1 / 6C (Numerical Analysis)

Insert an exact smooth solution into the linear multistep method and define its unscaled local truncation error byIts Taylor expansion has coefficientsThe method has order at least $p$ exactly when $c_0=\cdots =c_p=0$, since then ${\mathcal{R}}_h[y]=O(h^{p+1})$ for every sufficiently smooth $y$. Necessity can be checked by inserting polynomials of successive degrees. This also corresponds to the usual truncation error ${\mathcal{R}}_h/h=O(h^p)$.

On the other hand, expansion of the characteristic polynomials of a linear multistep method givesThus the exponential-symbol order criterion for a multistep method isIf “order $p$” is meant to be the exact rather than guaranteed order, one additionally requires $c_{p+1}\ne 0$. The first two vanishing conditions are the usual consistency of a numerical method relations $\rho (1)=0$, ${\rho }^{\prime }(1)=\sigma (1)$.

The relevant one of the characteristic polynomials of a linear multistep method isThe coefficients on the right are chosen to give third order, so consistency of a numerical method is already available. By the Dahlquist equivalence theorem, convergence with convergent starting values is equivalent to zero-stability: every root of $\rho$ must satisfy $|w|\le 1$, and every root on the unit circle must be simple.

If $a<0$, the quadratic is negative at $w=-1$ and tends to $+\infty$ as $w\to -\infty$, so it has a root less than $-1$. If $a>1$, its root product is $a$, so at least one root has modulus greater than one. For $0\le a\le 1/4$, its real roots $(-1\pm \sqrt{1-4a})/2$ lie in $[-1,0]$. For $1/4<a\le 1$, its conjugate roots have modulus $\sqrt{a}\le 1$.

At $a=0$, the three roots are $1,0,-1$, with the unit-modulus roots simple. At $a=1$, they are $1,e^{2\pi i/3},e^{-2\pi i/3}$, again simple. The repeated quadratic root at $a=1/4$ is $-1/2$, strictly inside the unit circle, which the root condition for a multistep method permits. ThereforeThis is the root stability of a cubic multistep polynomial; oscillatory parasitic roots at the endpoints do not invalidate convergence when their starting amplitudes tend to zero.