Текст подпрограммы и версий ( Фортран )
qtt3r.zip  qtt3d.zip  		Тексты тестовых примеров ( Фортран )
tqtt3r.zip  tqtt3d.zip 
Текст подпрограммы и версий ( Си )
qtt3r_c.zip  qtt3d_c.zip  		Тексты тестовых примеров ( Си )
tqtt3r_c.zip  tqtt3d_c.zip 
Текст подпрограммы и версий ( Паскаль )
qtt3r_p.zip  qtt3e_p.zip  		Тексты тестовых примеров ( Паскаль )
tqtt3r_p.zip  tqtt3e_p.zip 

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

Назначение

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

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

По формуле трапеций вычисляется значение тройного определенного интеграла от табличной функции f (x, y, z), заданной на равномерной сетке (xi, yj, zk),
xi = x1 + (i - 1) h1,    i = 1, ..., N,     yj = y1 + (j - 1) h2,   j = 1, ..., M,
zk = z1 + (k - 1) h3,   k = 1, ..., L .

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

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

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

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

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

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

          DIMENSION  F(11, 11, 11)
          H1 = 0.1
          H2 = 0.1
          H3 = 0.1
          N = 11
          M = 11
          L = 11
          Z = 0.
          DO 1 K = 1, L
          Y = 0.
          DO 2 J = 1, M
          X = 0.
          DO 3 I = 1, N
          F(I, J, K) = X + Y + Z
       3 X = X + H1
       2 Y = Y + H2
       1 Z = Z + H3
          CALL  QTT3R (RINT, H1, H2, H3, F, N, M, L)

Результат:

       RINT  =  1.50000