integrate_range_re Function

public function integrate_range_re(x, y, a, b) result(r)

Constructs a piecewise hermitian cubic interpolation of an array y(x) from discrete numerical data, and then integrates the interpolation in the range (a, b). The mesh spacing does not have to be uniform.

Arguments

Type IntentOptional 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

Return Value real(kind=wp)

Integral ∫y(x)·dx