Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=8) | :: | f(:) | ||||
real(kind=8) | :: | a | ||||
real(kind=8) | :: | b |
function c_trapz_ab_sample(f,a,b) result(sum) complex(8) :: f(:) real(8) :: a,b,dh complex(8) :: sum integer :: i,L L=size(f) dh=(b-a)/(L-1)/2.d0 sum=0d0 do i=1,L-1 sum = sum+(f(i+1)+f(i))*dh enddo end function c_trapz_ab_sample