5th-order ODE Solver More...
#include <OdeSolver.h>
Public Member Functions | |
OdeSolver (const SCA epsilon=SCA(0.001)) | |
Constuctor. Argument sets absolute and relative error tolerance. | |
~OdeSolver (void) | |
Default Destructor. | |
AbstractSolutionPoint < my_alg_type > | solve (AbstractSolutionPoint< my_alg_type > &InitialPoint, POLYLIE< my_alg_type > &theEquationIn) |
Exponentiates theEquationIn and solves for initial point InitialPoint at time one. |
5th-order ODE Solver
Dormand-Prince method (a version of the Runge-Kutta-Fehlberg method RKF54) with adaptive step size control. Method as described in: J. C. Butcher: Numerical methods for ordinary differential equations, 2nd ed., Wiley, 2008. W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery: Numerical recipes, 3rd ed., Cambridge University Press, 2007.