Implicit midpoint rule (source code)

= Implicit midpoint rule
{wiki=Midpoint_method}

The implicit midpoint rule advances an autonomous ordinary differential equation by
$$
y_{n+1}=y_n+h f\left(\frac{y_n+y_{n+1}}2\right).
$$
It is a second-order, <A-stable>, <time-symmetric numerical method>.