Текст подпрограммы и версий ( Фортран ) sf40r.zip sf40d.zip |
Тексты тестовых примеров ( Фортран ) tsf40r.zip tsf40d.zip |
Текст подпрограммы и версий ( Си ) sf40r_c.zip sf40d_c.zip |
Тексты тестовых примеров ( Си ) tsf40r_c.zip tsf40d_c.zip |
Текст подпрограммы и версий ( Паскаль ) sf40r_p.zip , sf40e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf40r_p.zip , tsf40e_p.zip |
Вычисление неполной гамма - функции Гx (p) вещественного аргумента и вещественного параметра.
SF40R вычисляет неполную гамма - функцию Гx (p)
x Гx (p) = ∫ e - t t p - 1 dt , 0
где x и p - принимают вещественные значения, причем x ≥ 0 и p > 0.
M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions New York, Dover Publications, INC., 1965.
REAL FUNCTION SF40R (X, P, IERR)
Параметры
X - | заданное значение аргумента x (тип: вещественный); |
P - | заданное значение параметра p (тип: вещественный); |
IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
IERR=65 - | когда значение аргумента x меньше 0; значение функции полагается равным 3.4E38; |
IERR=66 - | когда значение параметра P неполной гамма - функции не принадлежит полуинтервалу (0, 21.4] ; значение функции полагается равным 3.4E38. |
Версии
SF40D - | вычисление неполной гамма - функции Гx (p) с повышенной точностью. |
Вызываемые подпрограммы
SF39R - | подпрограмма вычисления гамма - функции распределения вероятностей I (x, p); |
SF39D - | подпрограмма вычисления гамма - функции распределения вероятностей I (x, p) с повышенной точностью; |
SF23C - | подпрограмма вычисления гамма - функции и ее натурального логарифма от комплексного аргумента; |
SF23P - | подпрограмма вычисления гамма - функции и ее натурального логарифма от комплексного аргумента с повышенной точностью; |
UTSF10 - UTSF12 | подпрограммы выдачи диагностических сообщений при работе функций SF23C, SF39R, SF40R; |
UTSF11 - UTSF13 | подпрограммы выдачи диагностических сообщений при работе функций SF23P, SF39D, SF40D |
Замечания по использованию
Для функции SF40D параметры X и P (0 < P ≤ 408.909) должны иметь тип DOUBLE PRECISION и при IERR ≠ 0 значение SF40D полагается равным 1.7D308. |
X = 1. P = 1. Y = SF40R (X, P, IERR) Результаты: Y = 0.6321205588074 , IERR = 0