Текст подпрограммы и версий ( Фортран )
qttvr.zip  qttvd.zip 
Тексты тестовых примеров ( Фортран )
tqttvr.zip  tqttvd.zip 
Текст подпрограммы и версий ( Си )
qttvr_c.zip  qttvd_c.zip 
Тексты тестовых примеров ( Си )
tqttvr_c.zip  tqttvd_c.zip 
Текст подпрограммы и версий ( Паскаль )
qttvr_p.zip  qttve_p.zip 
Тексты тестовых примеров ( Паскаль )
tqttvr_p.zip  tqttve_p.zip 

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

Назначение

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

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

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

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

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

    SUBROUTINE  QTTVR (RINT, X, F, N) 

Параметры

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

Версии

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

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

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

 

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

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

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

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

Результаты:

       RINT(1)    =  0.00000
       RINT(2)    =  0.00012

       RINT(49)  =  1.99610
       RINT(50)  =  1.99700