rkv98r Module Procedure

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

Verner's "most robust" Runge-Kutta (16:9(8)) pair.

Reference

J.H. Verner, "Explicit Runge-Kutta Methods with Estimates of the Local Truncation Error", SIAM Journal on Numerical Analysis, Vol. 15, Issue 4 (1978)
J.H. Verner, "Numerically optimal Runge--Kutta pairs with interpolants", Numerical Algorithms, 53, (2010) pp. 383-396. doi
Jim Verner's Refuge for Runge-Kutta Pairs
Coefficients

Note

The dense output (interpolant) is in rkv98r_dense_coefficients.

Arguments

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