Текст подпрограммы и версий ( Фортран )
qtt2r.zip  qtt2d.zip 
Тексты тестовых примеров ( Фортран )
tqtt2r.zip  tqtt2d.zip 
Текст подпрограммы и версий ( Си )
qtt2r_c.zip  qtt2d_c.zip 
Тексты тестовых примеров ( Си )
tqtt2r_c.zip  tqtt2d_c.zip 
Текст подпрограммы и версий ( Паскаль )
qtt2r_p.zip  qtt2e_p.zip 
Тексты тестовых примеров ( Паскаль )
tqtt2r_p.zip  tqtt2e_p.zip 

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

Назначение

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

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

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

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

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

    SUBROUTINE  QTT2R (RINT, X, F, N) 

Параметры

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

Версии

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

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

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

 

Требуется, чтобы  xi > xi - 1,   i = 2, ..., N;   N ≥ 2 .

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

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

          DIMENSION  F(50), X(50)
          H = 3.14159 / 49.
          N = 50
          X(1) = 0.
          X(2) = H / 4.
          DO 1 I = 3, N
       1 X(I) = X(I - 2) + 2.*H
          DO 2 I = 1, N
          X1 = X(I)
       2 F(I) = SIN(X1)
          CALL  QTT2R (RINT, X, F, N)

Результат:

       RINT  =  1.99700357