66 G4double momentum_mag_square = y[3]*y[3] + y[4]*y[4] + y[5]*y[5];
67 G4double inv_momentum_magnitude = 1.0 / std::sqrt( momentum_mag_square );
74 dydx[0] = y[3]*inv_momentum_magnitude;
75 dydx[1] = y[4]*inv_momentum_magnitude;
76 dydx[2] = y[5]*inv_momentum_magnitude;
78 dydx[3] = G[0]*cof1*cof2/
c_light;
79 dydx[4] = G[1]*cof1*cof2/
c_light;
80 dydx[5] = G[2]*cof1*cof2/
c_light;
84 dydx[7] = inverse_velocity;