A-stable, |R(infty)| = 1/3.
This method does not provide a dense output formula.
Hairer's code ros4.f