Текст подпрограммы и версий ( Фортран )
qtt1r.zip  qtt1d.zip 
Тексты тестовых примеров ( Фортран )
tqtt1r.zip  tqtt1d.zip 
Текст подпрограммы и версий ( Си )
qtt1r_c.zip  qtt1d_c.zip 
Тексты тестовых примеров ( Си )
tqtt1r_c.zip  tqtt1d_c.zip 
Текст подпрограммы и версий ( Паскаль )
qtt1r_p.zip  qtt1e_p.zip 
Тексты тестовых примеров ( Паскаль )
tqtt1r_p.zip  tqtt1e_p.zip 

Подпрограмма:  QTT1R

Назначение

Вычисление определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле трапеций.

Математическое описание

Вычисляется значение определенного интеграла на отрезке [x1, xN] от табличной функции f (x), заданной на равномерной сетке xi = x1 + (i - 1) h,   i = 1, ..., N, по формуле трапеций.

Н.С.Бахвалов. Численные методы, "Наука", M., 1975.

Использование

    SUBROUTINE  QTT1R (RINT, H, F, N) 

Параметры

RINT - вещественная переменная, содержащая вычисленное значение интеграла;
H - заданный шаг равномерной сетки (тип: вещественный);
F - вещественный вектоp длины N, содержащий значения функции f (x);
N - заданное число узлов сетки (тип: целый).

Версии

QTT1D - вычисление с удвоенной точностью определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле трапеций.

Вызываемые подпрограммы: нет

Замечания по использованию

 

Требуется, чтобы  N ≥ 2.

В подпрограмме QTT1D параметры RINT, H, F имеют тип DOUBLE PRECISION.

Пример использования

          DIMENSION  F(50)
          H = 3.14159 / 49.
          N = 50
          X = 0.
          DO 1 I = 1, N
          F(I) = SIN(X)
       1 X = X + H
          CALL  QTT1R (RINT, H, F, N)

Результат:

       RINT  =  1.99931