Absolute stability region
Stable region in red.

Coefficients

\begin{align} y_{n + 4} - y_{n + 3} = \frac{55}{24}h f(y_{n + 3}) - \frac{59}{24}h f(y_{n + 2}) + \frac{37}{24}h f(y_{n + 1}) - \frac{3}{8}h f(y_{n})\end{align}

Properties