d_trapz_nonlin_sample Function

function d_trapz_nonlin_sample(f, x) result(sum)

Arguments

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

Return Value real(kind=8)


Source Code

function d_trapz_nonlin_sample(f,x) result(sum)
  real(8) :: f(:)
  real(8) :: x(size(f))
  real(8) :: a,b,dh
  real(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 d_trapz_nonlin_sample