Текст подпрограммы и версий ( Фортран ) qts1r.zip qts1d.zip |
Тексты тестовых примеров ( Фортран ) tqts1r.zip tqts1d.zip |
Текст подпрограммы и версий ( Си ) qts1r_c.zip qts1d_c.zip |
Тексты тестовых примеров ( Си ) tqts1r_c.zip tqts1d_c.zip |
Текст подпрограммы и версий ( Паскаль ) 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.
SUBROUTINE QTS1R (RINT, H, F, N)
Параметры
RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
H - | заданный шаг равномерной сетки (тип: вещественный); |
F - | вещественный вектоp длины N, содержащий значения функции f (x); |
N - | заданное число узлов сетки (тип: целый). |
Версии
QTS1D - | вычисление с удвоенной точностью определенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы N ≥ 3. В подпрограмме QTS1D параметры 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 QTS1R (RINT, H, F, N) Результат: RINT = 2.00000