Текст подпрограммы и версий
qtgcr_p.zip  qtgce_p.zip  		Тексты тестовых примеров
tqtgcr_p.zip  tqtgce_p.zip 

Подпрограмма:  QTGCR (модуль QTGCR_p)

Назначение

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

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

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

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

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

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


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

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

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

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

Unit TQTGCR_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, QTGCR_p;

function TQTGCR: String;

implementation

function TQTGCR: String;
var
N,I,_i :Integer;
H,X :Real;
F :Array [0..49] of Real;
RINT :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 := X+H;
 end;
QTGCR(RINT,H,F,N);
Result := Result + #$0D#$0A;
for _i:=0 to 49 do
 begin
  Result := Result + Format('%20.16f ',[RINT[_i]]);
  if ( ((_i+1) mod 4)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
UtRes('TQTGCR',Result);  { вывод результатов в файл TQTGCR.res }
exit;
end;

end.

Результаты:

      RINT(1)    =  0.00000
      RINT(2)    =  0.00205
      RINT(49)  =  1.99794
      RINT(50)  =  2.00000