Текст подпрограммы и версий qts1r_p.zip qts1e_p.zip |
Тексты тестовых примеров tqts1r_p.zip tqts1e_p.zip |
Вычисление определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона.
Вычисляется значение определенного интеграла на отрезке [x1, xN] от табличной функции f (x), заданной на равномерной сетке x = xi + (i - 1) h, i = 1, ..., N, по квадратурной формуле, точной для многочленов второй степени.
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
procedure QTS1R(var RINT :Real; H :Real; var F :Array of Real; N :Integer);
Параметры
RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
H - | заданный шаг равномерной сетки (тип: вещественный); |
F - | вещественный вектоp длины N, содержащий значения функции f (x); |
N - | заданное число узлов сетки (тип: целый). |
Версии
QTS1E - | вычисление с расширенной (Extended) точностью определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы N ≥ 3. В подпрограмме QTS1E параметры RINT, H, F имеют тип Extended. |
Unit TQTS1R_p; interface uses SysUtils, Math, { Delphi } Lstruct, Lfunc, UtRes_p, QTS1R_p; function TQTS1R: String; implementation function TQTS1R: 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 := X+H; end; QTS1R(RINT,H,F,N); Result := Result + Format(' %20.16f ',[RINT]) + #$0D#$0A; UtRes('TQTS1R',Result); { вывод результатов в файл TQTS1R.res } exit; end; end. Результат: RINT = 2.00000