|
Текст подпрограммы и версий ( Фортран ) qsr2r.zip |
Тексты тестовых примеров ( Фортран ) tqsr2r.zip |
|
Текст подпрограммы и версий ( Си ) qsr2r_c.zip |
Тексты тестовых примеров ( Си ) tqsr2r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) qsr2r_p.zip |
Тексты тестовых примеров ( Паскаль ) tqsr2r_p.zip |
Вычисление определенного интеграла по квадратурной формуле Pомберга с заданной абсолютной погрешностью.
QSR2R вычисляет интеграл
b
∫ f (x) dx
a
по квадратурной формуле Pомберга с заданной абсолютной погрешностью.
Бахвалов H.C. "Численные методы", т.1, "Hаука", 1973.
SUBROUTINE QSR2R (RINT, A, B, F, E, IERR)
Параметры
| RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
| A, B - | заданные нижний и верхний пределы интегрирования (тип: вещественный); |
| F - | имя вещественной подпрограммы - функции, вычисляющей подинтегральную функцию f (x); |
| E - | заданная абсолютная погрешность вычисления интеграла (тип: вещественный); |
| IERR - | целая переменная, указывающая причину окончания вычислений; если интеграл сосчитан, то IERR = 0, в противном случае IERR = 65. |
Версии: нет
Вызываемые подпрограммы
| UTQS11 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы QSR2R. |
Замечания по использованию
| Если интеграл не может быть вычислен с точностью E, выдается диагностическое сообщение "заданная точность не может быть достигнута". |
REAL FUNCTION F(X)
REAL X
IF (X) 1, 1, 2
1 F = 1.
GO TO 3
2 F = X / (EXP(X)-1.)
3 RETURN
END
EXTERNAL F
INTEGER IERR
REAL RINT, A, B, F, E
A = 0.
B = 1.
E = 0.0000000001
CALL QSR2R (RINT, A, B, F, E, IERR)
Результаты: RINT = 0.77750463411 , IERR = 0