34 #ifndef Tsitouras_RK45
35 #define Tsitouras_RK45
43 G4int numberOfVariables = 6,
G4int IntegratorOrder() const
G4double DistChord() const
G4TsitourasRK45 & operator=(const G4TsitourasRK45 &)
G4TsitourasRK45 * fAuxStepper
G4double * fLastInitialVector
void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
G4double * fLastFinalVector
void interpolate(const G4double yInput[], const G4double dydx[], G4double yOut[], G4double Step, G4double tau)
G4TsitourasRK45(G4EquationOfMotion *EqRhs, G4int numberOfVariables=6, G4bool primary=true)
void SetupInterpolation()
void Interpolate(const G4double yInput[], const G4double dydx[], const G4double Step, G4double yOut[], G4double tau)