60 G4double reducedMass = mass*nucleusMass/(mass+nucleusMass);
74 G4double yMod = std::sqrt(y[0]*y[0]+y[1]*y[1]+y[2]*y[2]);
88 dydx[3] = yMod == 0 ? 0 : -deriv*y[0]/yMod*
c_light;
89 dydx[4] = yMod == 0 ? 0 : -deriv*y[1]/yMod*
c_light;
90 dydx[5] = yMod == 0 ? 0 : -deriv*y[2]/yMod*
c_light;