Текст подпрограммы и версий ( Фортран ) qtt2r.zip qtt2d.zip |
Тексты тестовых примеров ( Фортран ) tqtt2r.zip tqtt2d.zip |
Текст подпрограммы и версий ( Си ) qtt2r_c.zip qtt2d_c.zip |
Тексты тестовых примеров ( Си ) tqtt2r_c.zip tqtt2d_c.zip |
Текст подпрограммы и версий ( Паскаль ) qtt2r_p.zip qtt2e_p.zip |
Тексты тестовых примеров ( Паскаль ) tqtt2r_p.zip tqtt2e_p.zip |
Вычисление определенного однократного интеграла от табличной функции, заданной на неравномерной сетке, по формуле трапеций.
Вычисляется значение определенного интеграла на отрезке [x1, xN] от табличной функции f (x), заданной на неравномерной сетке xi, i = 1, ..., N, по формуле трапеций.
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
SUBROUTINE QTT2R (RINT, X, F, N)
Параметры
RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
X - | вещественный вектоp длины N, содержащий узлы неравномерной сетки; |
F - | вещественный вектоp длины N, содержащий значения функции f (x); |
N - | заданное число узлов сетки (тип: целый). |
Версии
QTT2D - | вычисление с удвоенной точностью определенного однократного от табличной функции, заданной на неравномерной сетке, по формуле трапеций. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы xi > xi - 1, i = 2, ..., N; N ≥ 2 . В подпрограмме QTT2D параметры RINT, X, F имеют тип DOUBLE PRECISION. |
DIMENSION F(50), X(50) H = 3.14159 / 49. N = 50 X(1) = 0. X(2) = H / 4. DO 1 I = 3, N 1 X(I) = X(I - 2) + 2.*H DO 2 I = 1, N X1 = X(I) 2 F(I) = SIN(X1) CALL QTT2R (RINT, X, F, N) Результат: RINT = 1.99700357