Текст подпрограммы и версий qtg1r_p.zip qtg1e_p.zip |
Тексты тестовых примеров tqtg1r_p.zip tqtg1e_p.zip |
Вычисление определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Грегори.
Вычисляется значение определенного интеграла на отрезке [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