Текст подпрограммы и версий
qtt3r_p.zip  qtt3e_p.zip  		Тексты тестовых примеров
tqtt3r_p.zip  tqtt3e_p.zip 

Подпрограмма:  QTT3R (модуль QTT3R_p)

Назначение

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

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

По формуле трапеций вычисляется значение тройного определенного интеграла от табличной функции 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.

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

procedure QTT3R(var RINT :Real; H1 :Real; H2 :Real; H3 :Real;
                var F :Array of Real; N :Integer; M :Integer;
                L :Integer);


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

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

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

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

Unit TQTT3R_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, QTT3R_p;

function TQTT3R: String;

implementation

function TQTT3R: String;
var
N,M,L,K,J,I :Integer;
H1,H2,H3,Z,Y,X,RINT :Real;
F :Array [0..1330] of Real;
label
_3,_2,_1;
begin
Result := '';  { результат функции }
H1 := 0.1;
H2 := 0.1;
H3 := 0.1;
N := 11;
M := 11;
L := 11;
Z := 0.0;
for K:=1 to L do
 begin
  Y := 0.0;
  for J:=1 to M do
   begin
    X := 0.0;
    for I:=1 to N do
     begin
      F[(I-1)+(J-1)*11+(K-1)*121] := X+Y+Z;
_3:
      X := X+H1;
     end;
_2:
    Y := Y+H2;
   end;
_1:
  Z := Z+H3;
 end;
QTT3R(RINT,H1,H2,H3,F,N,M,L);
Result := Result + Format('%20.16f ',[RINT]) + #$0D#$0A;
UtRes('TQTT3R',Result);  { вывод результатов в файл TQTT3R.res }
exit;
end;

end.

Результат:

       RINT  =  1.50000