Текст подпрограммы и версий qtt4r_p.zip qtt4e_p.zip |
Тексты тестовых примеров tqtt4r_p.zip tqtt4e_p.zip |
Вычисление определенного двухкратного интеграла от табличной функции, заданной на равномерной сетке, по формуле трапеций.
По формуле трапеций вычисляется значение определенного двойного интеграла от табличной функции f (x, y), заданной на равномерной сетке узлов (xi, yj), xi = x1 + (i - 1) h1, i = 1, ..., N, yj = y1 + (j - 1) h2, j = 1 , ..., M .
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
procedure QTT4R(var RINT :Real; H1 :Real; H2 :Real; N :Integer; M :Integer; var F :Array of Real);
Параметры
RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
H1 - | заданный шаг равномерной сетки по x (тип: вещественный); |
H2 - | заданный шаг равномерной сетки по y (тип: вещественный); |
N - | заданное число узлов сетки по x (тип: целый); |
M - | заданное число узлов сетки по y (тип: целый); |
F - | вещественный двумерный массив размера N на M, содержащий значения функции f (x, y). |
Версии
QTT4E - | вычисление с расширенной (Extended) точностью определенного двухкратного интеграла от табличной функции, заданной на равномерной сетке, по формуле трапеций. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы N ≥ 2, M ≥ 2 . В подпрограмме QTT4E параметры RINT, H1, H2, F имеют тип Extended. |
Unit TQTT4R_p; interface uses SysUtils, Math, { Delphi } Lstruct, Lfunc, UtRes_p, QTT4R_p; function TQTT4R: String; implementation function TQTT4R: String; var N,M,J,I :Integer; H1,H2,Y,X,RINT :Real; F :Array [0..1999] of Real; label _2,_1; begin Result := ''; { результат функции } H1 := 3.14159/49.0; H2 := 3.14159/39.0; N := 50; M := 40; Y := 0.0; for J:=1 to M do begin X := Y; for I:=1 to N do begin F[(I-1)+(J-1)*50] := Sin(X); _2: X := X+H1; end; _1: Y := Y+H2; end; QTT4R(RINT,H1,H2,N,M,F); Result := Result + Format('%20.16f ',[RINT]) + #$0D#$0A; UtRes('TQTT4R',Result); { вывод результатов в файл TQTT4R.res } exit; end; end. Результат: RINT = 0.00001