yp[-1] = 2.0 * dydx[-1] - yp[-2]