yp[ii + 1] = tau[ii] * beta[ii] * dydx[ii]