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