Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/6.5.0/doc/ssclib/node16.html
Дата изменения: Tue Aug 16 01:58:13 2005
Дата индексирования: Sat Dec 22 08:06:19 2007
Кодировка:
Поисковые слова: integral

XMM-Newton Science Analysis System


ssclib (ssclib-3.24) [xmmsas_20050815_1803-6.5.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