naginterfaces.library.inteq.volterra_​weights

naginterfaces.library.inteq.volterra_weights(method, iorder, nomg, nwt)[source]

volterra_weights computes the quadrature weights associated with the Adams’ methods of orders three to six and the Backward Differentiation Formulae (BDF) methods of orders two to five. These rules, which are referred to as reducible quadrature rules, can then be used in the solution of Volterra integral and integro-differential equations.

For full information please refer to the NAG Library document for d05bw

https://www.nag.com/numeric/nl/nagdoc_27.3/flhtml/d05/d05bwf.html

Parameters
methodstr, length 1

The type of method to be used.

For Adams’ type formulae.

For Backward Differentiation Formulae.

iorderint

The order of the method to be used. The number of starting weights, is determined by and .

If , .

If , .

nomgint

The number of convolution weights, .

nwtint

, the number of columns in the starting weights.

Returns
omegafloat, ndarray, shape

Contains the first convolution weights.

lenswint

The number of rows in the weights .

swfloat, ndarray, shape

contains the weights , for , for , where is as defined in Notes.

Raises
NagValueError
(errno )

On entry, .

Constraint: or .

(errno )

On entry, .

Constraint: .

(errno )

On entry, .

Constraint: .

(errno )

On entry, and .

Constraint: if , .

(errno )

On entry, and .

Constraint: if , .

(errno )

On entry, , and .

Constraint: if , .

(errno )

On entry, , and .

Constraint: if , .

(errno )

On entry, , and , .

Constraint: if , .

(errno )

On entry, , and , .

Constraint: if , .

Notes

volterra_weights computes the weights and for a family of quadrature rules related to the Adams’ methods of orders three to six and the BDF methods of orders two to five, for approximating the integral:

with , for , for some given constant .

In (1), is a uniform mesh, is related to the order of the method being used and , are the starting and the convolution weights respectively. The mesh size is determined as , where and is the chosen number of convolution weights , for . A description of how these weights can be used in the solution of a Volterra integral equation of the second kind is given in Further Comments. For a general discussion of these methods, see Wolkenfelt (1982) for more details.

References

Lambert, J D, 1973, Computational Methods in Ordinary Differential Equations, John Wiley

Wolkenfelt, P H M, 1982, The construction of reducible quadrature rules for Volterra integral and integro-differential equations, IMA J. Numer. Anal. (2), 131–152