◆ Contx1_r()

Function Contx1_r	(	I As	Long,
		T As	Double,
		RCont() As	Double,
		ICont() As	Long
	)

Initial value problem of ordinary differential equations (Extrapolation method (GBS algorithm)) (Reverse communication version) (Interpolation for dense output)

NOTE - THIS PROGRAM IS DEPRECATED AND WILL BE REMOVED IN THE NEXT VERSION.

Purpose: This is the support routine which interpolates the solution for dense output when an initial value problem of a system of ordinary differential equations is solved by Odex_r.

Odex_r must be called with Iout = 1 to enable dense output.
When returned from Odex_r with IRev = 5 after each accepted step, the interpolated solution at any T within latest step interval can be computed by using this routine.

Returns: Double
Interpolated solution Y(I).

Parameters

[in]	I	Element number of solution to be interpolated. (0 <= I <= N - 1)
[in]	T	T at which the interpolated solution is computed. T must be in the interval of the last step. That is, T1 <= T <= T2 where T1 is the value of T on the previous return from Odex_r with IRev = 5, and T2 is the value on this time.
[in]	RCont()	Array RCont(LRCont - 1) (LRCont >= (2Km + 5)N) Control information for dense output.
[in]	ICont()	Array ICont(LICont - 1) (LICont >= N) Integer control information for dense output.

Reference: E. Hairer, S.P. Norsett and G. Wanner, "Solving Ordinary Differential Equations. Nonstiff Problems. 2nd edition", Springer Series in Computational Mathematics, Springer-Verlag (1993)

Example Program: See example program (2) of Odex_r.