c_trapz_nonlin_sample Function

function c_trapz_nonlin_sample(f, x) result(sum)

Arguments

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

Return Value complex(kind=8)


Source Code

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