c_trapz_dh_sample Function

function c_trapz_dh_sample(f, dh) result(sum)

Arguments

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

Return Value complex(kind=8)


Source Code

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