Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8) | :: | f(:) | ||||
real(kind=8) | :: | a | ||||
real(kind=8) | :: | b |
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