Past exam of the mathematics course of the University of Cambridge / 2017 / iii / Paper 341 / 1

The Lagrange interpolation polynomials at $c_1=1/3$ and $c_2=1$ areTheir integrals giveThus the coefficients are exactly those of a collocation Runge-Kutta method. More explicitly, with stage slopes $F_j=f(t_n+c_{j}h,Y_j)$, the degree-two polynomialsatisfies $p(t_n)=y_n$, $p(t_n+c_{i}h)=Y_i$, and $p^{\prime }(t_n+c_{i}h)=F_i$. These are the collocation Runge-Kutta method equations, and $y_{n+1}=p(t_n+h)$. For an implicit Runge-Kutta method, the nearby stages exist uniquely for sufficiently small $h$ under Lipschitz continuity; see stage solvability of an implicit Runge-Kutta method.

Write $e=(1,1{)}^T$, $c=Ae$, and interpret powers of $c$ componentwise. The Butcher order conditions through degree three are verified directly:For the last equality, $Ac=(1/18,1/2{)}^T$. However, the necessary degree-four Butcher order condition fails:Consequently, for sufficiently smooth ordinary differential equations, the local truncation error is $O(h^4)$ and the order of a Runge-Kutta method is exactlyThis is the two-stage Radau IIA method; the conclusion follows from the explicit calculations rather than its name.

For the Dahlquist test equation, put $z=h\lambda$. Eliminating the stages of the Runge-Kutta method gives the stability functionIts poles are $2\pm i\sqrt{2}$, strictly in the right half-plane. To prove A-stability directly, write $z=-r+iy$ with $r\ge 0$. With $D(z)=1-2z/3+z^2/6$ and $N(z)=1+z/3$, elementary expansion givesTherefore $|R(z)|\le 1$ throughout the closed left half-plane; the inequality is strict except at $z=0$. The method is A-stable. Moreover $R(z)\to 0$ as $|z|\to \infty$, so the Radau IIA method is also L-stable.