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