Текст подпрограммы и версий ( Фортран )
qtthr.zip  qtthd.zip  		Тексты тестовых примеров ( Фортран )
tqtthr.zip  tqtthd.zip 
Текст подпрограммы и версий ( Си )
qtthr_c.zip  qtthd_c.zip  		Тексты тестовых примеров ( Си )
tqtthr_c.zip  tqtthd_c.zip 
Текст подпрограммы и версий ( Паскаль )
qtthr_p.zip  qtthe_p.zip  		Тексты тестовых примеров ( Паскаль )
tqtthr_p.zip  tqtthe_p.zip 

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

Назначение

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

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

Вычисляются значения неопределенного двойного интеграла по прямоугольникам x1 ≤ x ≤ xi,  y1 ≤ y ≤ yj,   i = 1, 2, ..., N,   j = 1, 2, ..., M от табличной функции f (x, y), заданной на равномерной сетке 

     xi = x1 + (i - 1) h1,     yj = y1 + (j - 1) h2,

по формуле трапеций.

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

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

    SUBROUTINE  QTTHR (RINT, H1, H2, F, N, M) 
   Параметры
RINT - вещественный двумерный массив размера N на M, содержащий вычисленные значения неопределенного интеграла;
H1 - заданный шаг равномерной сетки по x (тип: вещественный);
H2 - заданный шаг равномерной сетки по y (тип: вещественный);
F - вещественный двумерный массив размера N на M, содержащий значения подинтегральной функции f (x, y);
N - заданное число узлов сетки по x (тип: целый);
M - заданное число узлов сетки по y (тип: целый). 
   Версии
QTTHD - вычисление с удвоенной точностью неопределенного двукратного интеграла от табличной функции, заданной на равномерной сетке, по формуле трапеций. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
 

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

В подпрограмме QTTHD параметры RINT, H1, H2, F имеют тип DOUBLE PRECISION.

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

          DIMENSION F(50, 40), RINT(50, 40)
          H1 = 3.14159 / 49.
          H2 = 3.14159 / 39.
          N = 50
          M = 40
          Y = 0.
          DO 2 J = 1, 40
          X = Y
          DO 1 I = 1, 50
          F(I, J) = SIN(X)
       1 X = X + H1
       2 Y = Y + H2
          CALL  QTTHR (RINT, H1, H2, F, N, M)

Результаты:

       RINT(1, 1)        =  0.00000
       RINT(2, 2)        =  0.00037
       RINT(N-1, M)  =  0.12803
       RINT(N, M)      =  0.00001