Текст подпрограммы и версий ( Фортран ) qttvr.zip qttvd.zip |
Тексты тестовых примеров ( Фортран ) tqttvr.zip tqttvd.zip |
Текст подпрограммы и версий ( Си ) qttvr_c.zip qttvd_c.zip |
Тексты тестовых примеров ( Си ) tqttvr_c.zip tqttvd_c.zip |
Текст подпрограммы и версий ( Паскаль ) qttvr_p.zip qttve_p.zip |
Тексты тестовых примеров ( Паскаль ) tqttvr_p.zip tqttve_p.zip |
Вычисление неопределенного однократного интеграла от табличной функции, заданной на неравномерной сетке, по формуле трапеций.
Вычисляются значения неопределенного интеграла на отрезках [x1, xi] i = 1, ..., N, от табличной функции f (x), заданной на неравномерной сетке xi, i = 1, ..., N, по формуле трапеций.
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
SUBROUTINE QTTVR (RINT, X, F, N)
Параметры
RINT - | вещественный вектоp длины N, содержащий вычисленные значения неопределенного интеграла; |
X - | вещественный вектоp длины N, содержащий узлы неравномерной сетки; |
F - | вещественный вектоp длины N, содержащий значения функции f (x); |
N - | заданное число узлов сетки (тип: целый). |
Версии
QTTVD - | вычисление с удвоенной точностью неопределенного однократного интеграла от табличной функции, заданной на неравномерной сетке, по формуле трапеций. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы xi > xi - 1, i = 2, ..., N; N ≥ 2 . В подпрограмме QTTVD параметры RINT, X, F имеют тип DOUBLE PRECISION. |
DIMENSION F(50), X1(50), RINT(50) H = 3.14159 / 49. N = 50 X1(1) = 0.0 X1(2) = H / 4. DO 1 I = 3, N, 2 X1(I) = X1(I - 2) + 2.*H 1 X1(I + 1) = X1(I - 1) + 2.*H DO 2 I = 1, N X = X1(I) 2 F(I) = SIN(X) CALL QTTVR (RINT, X1, F, N) Результаты: RINT(1) = 0.00000 RINT(2) = 0.00012 RINT(49) = 1.99610 RINT(50) = 1.99700