rk8_10 Module Procedure

module procedure rk8_10 module subroutine rk8_10(me, t, x, h, xf)

Take one Runge Kutta 8 integration step: t -> t+h (x -> xf) This is Formula (8-10) from Reference [1].

Reference

E. B. Shanks, "Higher Order Approximations of Runge-Kutta Type", NASA Technical Note, NASA TN D-2920, Sept. 1965.

Arguments

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

Calls

Help