beta = yp[1:] / dydx