ode45.h File Reference

#include "../utils.h"
#include "ode.h"
#include <functional>
#include <vector>

Go to the source code of this file.

Functions
ODEResult	ODE45 (const ODE &fun, const std::vector< double > &tInterval, const Matrix< double > &y0, double h=0.0)

Detailed Description

ODE Solver implementation based on Dormand-Prince Method to solve ordinary differential equations https://en.wikipedia.org/wiki/Dormand%E2%80%93Prince_method

Current state:

• implemented 7 function evaluations
• using 6/7 evaluations for 5th order RK

TODO: use calculated error of 5th order - 4th order to alter step width h dynamically to fully implement dp-method.

Requires:

#include <math/numerics/ode/ode45.h>

Function Documentation

ODE45()

ODEResult ODE45	(	const ODE &	fun,
		const std::vector< double > &	tInterval,
		const Matrix< double > &	y0,
		double	h = 0.0
	)

Implementation of 5th order Runge-Kutta-Method

Parameters

fun	ode to approximate
tInterval	interval to perform approximation on
y0	start value
h	step width for time values

Returns

approximated values

Examples

numerics/ode/TestODE45.cpp.