Текст подпрограммы и версий
qtt2r_p.zip  qtt2e_p.zip 
Тексты тестовых примеров
tqtt2r_p.zip  tqtt2e_p.zip 

Подпрограмма:  QTT2R (модуль QTT2R_p)

Назначение

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

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

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

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

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

procedure QTT2R(var RINT :Real; var X :Array of Real;
                var F :Array of Real; N :Integer);

Параметры

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

Версии

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

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

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

 

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

В подпрограмме QTT2E параметры RINT, X, F имеют тип Extended.

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

Unit TQTT2R_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, QTT2R_p;

function TQTT2R: String;

implementation

function TQTT2R: String;
var
N,I :Integer;
H,X1,RINT :Real;
F :Array [0..49] of Real;
X :Array [0..49] of Real;
label
_1,_2;
begin
Result := '';  { результат функции }
H := 3.14159/49.0;
N := 50;
X[0] := 0.0;
X[1] := H/4.0;
FОR I:=3 to N do
 begin
_1:
  Х[I-1] := X[I-3]+H*2.0;
 end;
for I:=1 to N do
 begin
  X1 := X[I-1];
_2:
  F[I-1] := Sin(X1);
 end;
QTT2R(RINT,X,F,N);
Result := Result + Format(' %20.16f ',[RINT]) + #$0D#$0A;
UtRes('TQTT2R',Result);  { вывод результатов в файл TQTT2R.res }
exit;
end;

end.

Результат:

       RINT  =  1.99700357