 Solution

Expand the second stage by the Taylor theorem, writing $f=f(t_n,y_n)$ and $D_{y}f$ for the derivative in the state variable:The update therefore agrees with the exact solution's Taylor series:Thus the one-step local truncation error is $O(h^3)$, and the method has order at least two under the usual smoothness and Lipschitz continuity assumptions.

For $f=\lambda y$, the two-stage Runge-Kutta method givesOn the real axis, $R(z)=\frac{1}{2}(z+1{)}^2+\frac{1}{2}\ge \frac{1}{2}$, so $|R(z)|\le 1$ is equivalent to $(z+1{)}^2\le 1$. HenceThe decay convention gives $(-2,0)$. Since negative real values such as $z=-3$ lie outside the linear stability domain, the method is not A-stable. Its order is exactly two, since $R(z)$ lacks the cubic term of $e^z$.