Текст подпрограммы и версий ( Фортран ) 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 |
Вычисление значений гамма - функции и 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.