yp[1:-1] = (dydx[:-1] + dydx[1:]) / 2.0