rkf45 Module Procedure

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

Fehlberg's 4(5) method.

This is Table III, RK4(5), Formula 2 in the reference.

References

E. Fehlberg, "Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems", NASA Technical Report R-315, July 1, 1969.

Arguments

Type	Intent	Optional	Attributes		Name
class(rkf45_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