Текст подпрограммы и версий ( Фортран ) qtscr.zip qtscd.zip |
Тексты тестовых примеров ( Фортран ) tqtscr.zip tqtscd.zip |
Текст подпрограммы и версий ( Си ) qtscr_c.zip qtscd_c.zip |
Тексты тестовых примеров ( Си ) tqtscr_c.zip tqtscd_c.zip |
Текст подпрограммы и версий ( Паскаль ) qtscr_p.zip qtsce_p.zip |
Тексты тестовых примеров ( Паскаль ) tqtscr_p.zip tqtsce_p.zip |
Вычисление неопределенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона.
Вычисляются значения неопределенного интеграла на отрезках [x1, xi] i = 1, ..., N, от табличной функции f (x), заданной на равномерной сетке xi = x1 + (i - 1) h, i = 1, ..., N, по квадратурной формуле, точной для многочленов второй степени.
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
SUBROUTINE QTSCR (RINT, H, F, N)
Параметры
RINT - | вещественный вектоp длины N, содержащий вычисленные значения неопределенного интеграла; |
H - | заданный шаг равномерной сетки (тип: вещественный); |
F - | вещественный вектоp длины N, содержащий значения функции f (x); |
N - | заданное число узлов сетки (тип: целый). |
Версии
QTSCD - | вычисление с удвоенной точностью неопределенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы N ≥ 3. В подпрограмме QTSCD параметры RINT, H, F имеют тип DOUBLE PRECISION. |
DIMENSION F(50), X(50), RINT(50) H = 3.14159 / 49. N = 50 X = 0. DO 1 I = 1, N F(I) = SIN(X) 1 X = X + H CALL QTSCR (RINT, H, F, N) Результаты: RINT(1) = 0.00000 RINT(2) = 0.00205 RINT(49) = 1.99794 RINT(50) = 2.00000