Текст подпрограммы и версий
qtg1r_p.zip  qtg1e_p.zip 
Тексты тестовых примеров
tqtg1r_p.zip  tqtg1e_p.zip 

Подпрограмма:  QTG1R (модуль QTG1R_p)

Назначение

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

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

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

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

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

procedure QTG1R(var RINT :Real; H :Real; var F :Array of Real;
                N :Integer);

Параметры

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

Версии

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

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

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

 

Требуется, чтобы   N ≥ 6.

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

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

Unit TQTG1R_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, QTG1R_p;

function TQTG1R: String;

implementation

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

end.

Результат:

       RINT  =  2.00000