rkdp65 Module Procedure

module procedure rkdp65 module subroutine rkdp65(me, t, x, h, xf, xerr)

Dormand-Prince 6(5) method. This is RK6(5)8M from the reference.

Reference

P.J. Prince, J.R. Dormand, "High order embedded Runge-Kutta formulae", Journal of Computational and Applied Mathematics, Volume 7, Issue 1, March 1981, Pages 67-75.

Arguments

Type	Intent	Optional	Attributes		Name
class(rkdp65_class),	intent(inout)			::	me
real(kind=wp),	intent(in)			::	t	initial time
real(kind=wp),	intent(in),		dimension(me%n)	::	x	initial state
real(kind=wp),	intent(in)			::	h	time step
real(kind=wp),	intent(out),		dimension(me%n)	::	xf	state at time `t+h`
real(kind=wp),	intent(out),		dimension(me%n)	::	xerr	truncation error estimate for `x`

Calls

Help