Текст подпрограммы и версий ( Фортран )
sf23c.zip  sf23p.zip 
Тексты тестовых примеров ( Фортран )
tsf23c.zip  tsf23p.zip 
Текст подпрограммы и версий ( Си )
sf23c_c.zip  sf23p_c.zip 
Тексты тестовых примеров ( Си )
tsf23c_c.zip  tsf23p_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf23c_p.zip , sf23z_p.zip
Тексты тестовых примеров ( Паскаль )
tsf23c_p.zip , tsf23z_p.zip

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

Назначение

Вычисление значений гамма - функции и ee натурального логарифма от комплексного аргумента.

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

Подпрограмма SF23C вычисляет значение гамма-функции

                   
     Г (Z)  =  ∫ t z - 1 e - t dt
                   0 

и ее натурального логарифма для комплексных значений Z.

H.Kuki, Complex Gamma Function with Error Control, Communications ACM, 15, 4, 1972.

Использование

    SUBROUTINE  SF23C (ARGR, ARGI, ANSR, ANSI, EPS, L, IERR) 

Параметры

ARGR - заданное значение действительной части аpгумента Z (тип: вещественный);
ARGI - заданное значение мнимой части аргумента Z (тип: вещественный);
         ANSR -
         ANSI  
вещественные переменные, значения которых в pезультате работы подпрограммы полагаются равными действительной и мнимой части pезультата соответственно;
EPS - вещественная переменная, определяющая на входе абсолютную точность, с которой надо найти результат; в результате работы подпрограммы значение EPS полагается равным действительно достигнутой точности;
L - задает режим работы подпрограммы (тип: целый); при этом:
L = 0 - если вычисляется натуральный логарифм от гамма - функции;
L = 1 - если вычисляется гамма - функция;
IERR - целая переменая, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; при этом:
IERR=65 - когда мнимая часть аргумента pавна нулю, а вещественная часть аргумента больше 21.4;
IERR=66 - когда вычисляемая функция превосходит максимально допустимое в машине число;
IERR=67 - когда значение аргумента pавно нулю или любому целому отрицательному числу;
  при IERR ≠ 0 переменным ANSR и ANSI присваивается значение 3.4E38.

Версии

SF23P - вычисление значений гамма - функции и ее натурального логарифма в комплексной области с повышенной точностью.

Вызываемые подпрограммы

UTSF10 - подпрограмма печати диагностических сообщений при работе подпрограмы SF23C.
UTSF11 - подпрограмма печати диагностических сообщений при работе подпрограммы SF23P.

Замечания по использованию

 

Для подпрограммы SF23P параметры ARGR, ARGI, ANSR, ANSI имеют тип DOUBLE PRECISION и IERR полагается равным 65, если значение ARGI pавно нулю и значение ARGR больше 408.909. При IERR ≠ 0 значения ARGR и ARGI полагаются равными 1.7D308.

Пример использования

         ARGR = 6.7
         ARGI = 0.3
         EPS = 1.E-12
         L = 0
         CALL  SF23C (ARGR, ARGI, ANSR1, ANSI1, EPS, L, IERR)

Результаты:

       ANSR1  =  6.017194153159;
       ANSI1  =  0.547804744098;
       EPS       =  3.2725E-11;
       IERR     =  0.