void DC2DT6::BodyRHS(UserData<real_t>& ud) { eRHS(4) += _a3*ud.BodyForce(_x[3],_time); eRHS(5) += _a3*ud.BodyForce(_x[4],_time); eRHS(6) += _a3*ud.BodyForce(_x[5],_time); }