yp[ii] = tau[ii] * alpha[ii] * dydx[ii]