d_trapz_dh_sample Function

function d_trapz_dh_sample(f, dh) result(sum)

Arguments

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

Return Value real(kind=8)


Source Code

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