Текст подпрограммы и версий ( Фортран )
qsj1r.zip 
Тексты тестовых примеров ( Фортран )
tqsj1r.zip 
Текст подпрограммы и версий ( Си )
qsj1r_c.zip 
Тексты тестовых примеров ( Си )
tqsj1r_c.zip 
Текст подпрограммы и версий ( Паскаль )
qsj1r_p.zip 
Тексты тестовых примеров ( Паскаль )
tqsj1r_p.zip 

Подпрограмма:  QSJ1R

Назначение

Вычисление определенного интеграла по формуле Гаусса - Якоби.

Математическое описание

 Подпрограмма 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