d_trapz_ab_sample Function

function d_trapz_ab_sample(f, a, b) result(sum)

Arguments

Type IntentOptional Attributes Name
real(kind=8) :: f(:)
real(kind=8) :: a
real(kind=8) :: b

Return Value real(kind=8)


Source Code

function d_trapz_ab_sample(f,a,b) result(sum)
  real(8) :: f(:)
  real(8) :: a,b,dh
  real(8) :: sum
  integer :: i,L
  L=size(f)
  dh=(b-a)/(L-1)/2d0
  sum=0d0
  do i=1,L-1
     sum = sum+(f(i+1)+f(i))*dh
  enddo
end function d_trapz_ab_sample