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