Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 68 / 4

The tableau is the three-stage Lobatto IIIA method. To check its order of a Runge-Kutta method directly, write $e=(1,1,1{)}^T$, $c=Ae=(0,1/2,1{)}^T$, and $C=\mathrm{diag}(c)$. The eight fourth-order conditions for a Runge-Kutta method evaluate toPowers of $c$ here mean componentwise powers. For example, $Ac=(0,1/8,1/2{)}^T$, $A{c}^2=(0,1/24,1/3{)}^T$, and $A^{2}c=(0,1/48,1/6{)}^T$, which make the last three checks immediate. These Butcher order conditions establish order at least four for general nonlinear ordinary differential equations.

It is not order five: part (b)'s stability function has expansionThe fifth coefficient already fails on the Dahlquist test equation. The order is exactly four.

For a Runge-Kutta method, the stability function is $R(z)=1+z{b}^T(I-zA{)}^{-1}e$. Solving the stage equations on the Dahlquist test equation givesThe denominator has roots of a polynomial $3\pm i\sqrt{3}$, both in the open right half-plane. If $z=x+iy$, direct expansion givesSince $|z{|}^2+12>0$, the modulus of $R$ is at most one exactly when $x\le 0$. ThusThe method is A-stable. Its stability function has unit modulus on the imaginary axis and tends to one as $z\to -\infty$, so it is not L-stable.

The criterion for algebraic stability of a Runge-Kutta method is $b_i\ge 0$ and positive semidefiniteness ofThe Runge-Kutta contractivity identity implies that this criterion is sufficient for B-stability, namely contractivity for a dissipative vector field, whenever the stage equations are defined. The implication from algebraic stability to B-stability is the relevant nonlinear theorem; A-stability alone is a linear property.

All weights here are positive, butA positive semidefinite matrix cannot have this negative quadratic value. Hence the method is not algebraically stable. Failure of this sufficient criterion alone would not be a proof of failure of B-stability; no such converse is needed for the question.