Midpoint Method

Evaluates a function $g$ at each midpoint between timesteps:

$$y_{i+1}\approx y_i+\Delta t\cdot g\left(t_i+\frac{\Delta t}{2},y_i+g(t_i,y_i)\frac{\Delta t}{2}\right)$$