alpha = yp[:-1] / dydx