51 G4int numberOfVariables):
54 unsigned int noVariables=
std::max(numberOfVariables,8);
84 yTemp[7] = yOut[7] = yIn[7];
86 for( i = 0; i < numberOfVariables; i++ )
88 yTemp[i] = yIn[i] + h*dydx[i] ;
93 for( i = 0; i < numberOfVariables; i++ )
95 yOut[i] = yIn[i] + 0.5 * h * ( dydx[i] +
dydxTemp[i] );
T max(const T t1, const T t2)
brief Return the largest of the two arguments
void DumbStepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[])
G4ImplicitEuler(G4EquationOfMotion *EqRhs, G4int numberOfVariables=6)
void RightHandSide(const double y[], double dydx[]) const
G4int GetNumberOfVariables() const