|
Текст подпрограммы и версий ( Фортран ) 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