Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/8.0.0/doc/ssclib/node16.html
Дата изменения: Wed Jul 2 03:52:39 2008
Дата индексирования: Fri Sep 5 18:38:06 2008
Кодировка:

XMM-Newton Science Analysis System


ssclib (ssclib-3.33.4) [xmmsas_20080701_1801-8.0.0]

Meta Index / Home Page / General-purpose array-processing utilities / Array or vector scalar

Trapezoid-rule, 1-D, numerical integral

Given two vectors $x$ and $y$ containing samples of a function $F$, this subroutine returns the `trapezoidal-rule' estimate of the integral of $F$. In other words, $F$ is replaced by the set of piecewise-continuous line segments defined by $x$ and $y$. Note that $x$ is assumed to be monotonically increasing - if it isn't, you'll get strange results.

  interface trapezoidSum
    function trapezoidSumSingle(x, y) result(summ)
      real(single), intent(in) :: x(:), y(size(x))
      real(single)             :: summ
    end function trapezoidSumSingle

    function trapezoidSumDouble(x, y) result(summ)
      real(double), intent(in) :: x(:), y(size(x))
      real(double)             :: summ
    end function trapezoidSumDouble
  end interface