Module Owl_maths_quadrature

Numerical Integration

Integration functions
val trapz : ?n:int -> ?eps:float -> ( float -> float ) -> float -> float -> float

trapz f a b computes the integral of f on the interval [a,b] using the trapezoidal rule, i.e. :math:`\int_a^b f(x) dx`.

Parameters: * f: function to be integrated. * n: the maximum allowed number of steps. The default value is 20. * eps: the desired fractional accuracy. The default value is 1e-6. * a: lower bound of the integrated interval. * b: upper bound of the integrated interval.

Returns: * y: the integral of f on [a, b].

val simpson : ?n:int -> ?eps:float -> ( float -> float ) -> float -> float -> float

simpson f a b computes the integral of f on the interval [a,b] using the Simpson's rule, i.e. :math:`\int_a^b f(x) dx`.

Parameters: * f: function to be integrated. * n: the maximum allowed number of steps. The default value is 20. * eps: the desired fractional accuracy. The default value is 1e-6. * a: lower bound of the integrated interval. * b: upper bound of the integrated interval.

Returns: * y: the integral of f on [a, b].

val romberg : ?n:int -> ?eps:float -> ( float -> float ) -> float -> float -> float

romberg f a b computes the integral of f on the interval [a,b] using the Romberg method, i.e. :math:`\int_a^b f(x) dx`. Note that this algorithm is much faster than trapz and simpson.

Parameters: * f: function to be integrated. * n: the maximum allowed number of steps. The default value is 20. * eps: the desired fractional accuracy. The default value is 1e-6. * a: lower bound of the integrated interval. * b: upper bound of the integrated interval.

Returns: * y: the integral of f on [a, b].

val gaussian_fixed : ?n:int -> ( float -> float ) -> float -> float -> float

gaussian_fixed f a b computes the integral of f on the interval [a,b] using the Gaussian quadrature of fixed order. Note that this algorithm is much faster than others due to cached weights.

Parameters: * f: function to be integrated. * n: the order of polynomial. The default value is 10. * a: lower bound of the integrated interval. * b: upper bound of the integrated interval.

Returns: * y: the integral of f on [a, b].

val gaussian : ?n:int -> ?eps:float -> ( float -> float ) -> float -> float -> float

gaussian f a b computes the integral of f on the interval [a,b] using adaptive Gaussian quadrature of fixed tolerance.

Parameters: * f: function to be integrated. * n: the maximum order. The default value is 50. * eps: the desired fractional accuracy. The default value is 1e-6. * a: lower bound of the integrated interval. * b: upper bound of the integrated interval.

Returns: * y: the integral of f on [a, b].

Helper functions
val trapzd : ( float -> float ) -> float -> float -> int -> float

The function computes the nth stage of refinement of an extended trapezoidal rule. It is the workhorse of several integration functions including trapz, simpson, and romberg.

Parameters: * f: function to be integrated. * a: lower bound of the integrated interval. * b: upper bound of the integrated interval. * n: the nth stage.

Returns: * y: the integral of f on [a, b].

val gauss_legendre : ?eps:float -> ?a:float -> ?b:float -> int -> float array * float array

Given the lower and upper limits of integration a and b, and order n, the function computes the abscissas and weights of the Gauss-Legendre n-point quadrature formula.