This function constructs a piecewise hermitian cubic interpolation of an array y(x) based on discrete numerical data, and subsequently evaluates the integral of the interpolation in the range (a,b). Note that the mesh spacing of x does not necessarily have to be uniform.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(:) | :: | x | Variable x |
|
real(kind=wp), | intent(in), | dimension(size(x)) | :: | y | Function y(x) |
|
real(kind=wp), | intent(in) | :: | a | Left endpoint |
||
real(kind=wp), | intent(in) | :: | b | Right endpoint |
Integral ∫y(x)·dx