|
Текст подпрограммы и версий ( Фортран ) qsl1r.zip qsl1d.zip |
Тексты тестовых примеров ( Фортран ) tqsl1r.zip tqsl1d.zip |
|
Текст подпрограммы и версий ( Си ) qsl1r_c.zip qsl1d_c.zip |
Тексты тестовых примеров ( Си ) tqsl1r_c.zip tqsl1d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) qsl1r_p.zip qsl1e_p.zip |
Тексты тестовых примеров ( Паскаль ) tqsl1r_p.zip tqsl1e_p.zip |
Вычисление определенного интеграла по формулам Лобатто.
QSL1R вычисляет интеграл
B
∫ f (x) dx
A
с заданной абсолютной погрешностью E по формулам Лобатто по семи узлам (формула 11 - й степени точности) с автоматическим выбором шага.
О.В.Брушлинская, Л.Г.Васильева. Набор стандартных программ приближенного вычисления однократных интегралов с автоматическим выбором шага. Сб. "Численный анализ на ФОРТРАНе", вып. 8, Изд-во МГУ, 1974.
SUBROUTINE QSL1R (RINT, A, B, F, E)
Параметры
| RINT - | переменная, содержащая вычисленное значение интеграла (тип: вещественный); |
| A, B - | заданные нижний и верхний пределы интегрирования (тип: вещественный); |
| F - | имя вещественной подпрограммы - функции вычисления подинтегральной функции f (x); |
| E - | заданная абсолютная погрешность вычисления интеграла (тип: вещественный). |
Версии
| QSL1D - | вычисление с удвоенной точностью определенного интеграла по формулам Лобатто. |
Вызываемые подпрограммы: нет
Замечания по использованию
|
Подпрограмма не гарантирует заданную погрешность, хотя в большинстве случаев интегралы будут вычислены с тpебуемой точностью. В подпрограмме QSL1D параметры RINT, A, B, F, E имеют тип DOUBLE PRECISION. |
FUNCTION F(X)
IF (X) 1, 2, 1
2 F = 0.
GO TO 3
1 F = 1./SQRT(X)
3 RETURN
END
EXTERNAL F
REAL RINT
A = 0.
B = 1.
E = 0.00000001
CALL QSL1R (RINT, A, B, F, E)
Результат: RINT = 1.9999963224