Текст подпрограммы и версий ( Фортран ) qtg1r.zip qtg1d.zip |
Тексты тестовых примеров ( Фортран ) tqtg1r.zip tqtg1d.zip |
Текст подпрограммы и версий ( Си ) qtg1r_c.zip qtg1d_c.zip |
Тексты тестовых примеров ( Си ) tqtg1r_c.zip tqtg1d_c.zip |
Текст подпрограммы и версий ( Паскаль ) 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.
SUBROUTINE QTG1R (RINT, H, F, N)
Параметры
RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
H - | заданный шаг равномерной сетки (тип: вещественный); |
F - | вещественный вектоp длины N, содержащий значения функции f (x); |
N - | заданное число узлов сетки (тип: целый). |
Версии
QTG1D - | вычисление с удвоенной точностью определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Грегори. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы N ≥ 6. В подпрограмме QTG1D параметры RINT, H, F имеют тип DOUBLE PRECISION. |
DIMENSION F(50) H = 3.14159 / 49. N = 50 X = 0. DO 1 I = 1, N F(I) = SIN(X) 1 X = X + H CALL QTG1R (RINT, H, F, N) Результат: RINT = 2.00000