trapz Function

  1. korc_synthetic_camera.f90
  2. korc_synthetic_camera
  3. trapz

private function trapz(x, f)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in), DIMENSION(:):: x
real(kind=rp), intent(in), DIMENSION(:):: f

Return Value real(kind=rp)


Contents

Source Code

trapz

Source Code

  FUNCTION trapz(x,f)
    IMPLICIT NONE
    REAL(rp), DIMENSION(:), INTENT(IN) :: x
    REAL(rp), DIMENSION(:), INTENT(IN) :: f
    REAL(rp) :: trapz
    INTEGER :: N

    N = SIZE(x)

    trapz = 0.5_rp*SUM( (x(2:N) - x(1:N-1))*(f(1:N-1) + f(2:N)) )
  END FUNCTION trapz