Текст подпрограммы и версий ( Фортран ) qsj1r.zip |
Тексты тестовых примеров ( Фортран ) tqsj1r.zip |
Текст подпрограммы и версий ( Си ) qsj1r_c.zip |
Тексты тестовых примеров ( Си ) tqsj1r_c.zip |
Текст подпрограммы и версий ( Паскаль ) qsj1r_p.zip |
Тексты тестовых примеров ( Паскаль ) tqsj1r_p.zip |
Вычисление определенного интеграла по формуле Гаусса - Якоби.
Подпрограмма QSJ1R вычисляет определенный интеграл B ∫ ( 1 - x )α ( 1 + x ) β f (x) dx A по квадратурной формуле Гаусса - Якоби.
SUBROUTINE QSJ1R (RINT, A, B, ALF, BTA, F, N)
Параметры
RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
A, B - | заданные нижний и верхний пределы интегрирования, соответственно (тип: вещественный); |
ALF,BTA - | заданные параметры α и β соответственно, в весовом множителе (1 - x) α (1 + x) β (тип: вещественный); |
F - | имя вещественной подпрограммы - функции, вычисляющей подинтегральную функцию f (x); |
N - | заданное число узлов интегрирования (2 ≤ N ≤ 512, тип: целый); |
Версии: нет
Вызываемые подпрограммы
QSJAC - | подпрограмма, вычисляющая узлы и веса квадратурной формулы Гаусса - Якоби. |
Замечания по использованию
Число узлов интегрирования N на заданном отрезке [A; B] может быть задано лишь как целое из интервала 2 ≤ N ≤ 512. |
DIMENSION X(512), W(512) EXTERNAL F ALF = 1.0 BTA = 1.0 A = - 1.0 B = 1.0 N = 4 CALL QSJ1R (RINT, A, B, ALF, BTA, F, N) REAL FUNCTION F(X) REAL X F = X RETURN END Результат: RINT = 0.0000000000