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