|
Текст подпрограммы и версий ( Фортран ) sf39r.zip sf39d.zip |
Тексты тестовых примеров ( Фортран ) tsf39r.zip tsf39d.zip |
|
Текст подпрограммы и версий ( Си ) sf39r_c.zip sf39d_c.zip |
Тексты тестовых примеров ( Си ) tsf39r_c.zip tsf39d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf39r_p.zip , sf39e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf39r_p.zip , tsf39e_p.zip |
Вычисление гамма - функции распределения вероятностей I (x, p).
SF39R вычисляет вероятность того, что случайная величина, подчиненная гамма - распределению с параметром p, меньше или равна заданному значению x.
Другими словами, SF39R вычисляет отношение I(x, p)= Гx(p) / Г(p) =
x ∞
= ∫ e - t t p - 1 dt / ∫ e - t t p - 1 dt ,
0 0
где Гx (p) - неполная гамма - функция,
Г (p) - гамма - функция (полная),
x и p - принимают вещественные значения, причем x ≥ 0 и p > 0.
M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, New York, Dover Publications, INC., 1965.
REAL FUNCTION SF39R (X, P, IERR)
Параметры
| X - | заданное значение аргумента x (тип: вещественный); |
| P - | заданное значение параметра p (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда значение аргумента X меньше 0; значение функции полагается равным 3.4E38; |
| IERR=66 - | когда значение параметра P гамма - распределения меньше или pавно 0; значение функции полагается равным 3.4E38. |
Версии
| SF39D - | вычисление гамма - функции распределения вероятностей I (x, p) с повышенной точностью. |
Вызываемые подпрограммы
| UTSF12 - | подпрограмма выдачи диагностических сообщений при работе функции SF39R. |
| UTSF13 - | подпрограмма выдачи диагностических сообщений при работе функции SF39D. |
Замечания по использованию
| Для функции SF39D параметры X и P должны иметь тип DOUBLE PRECISION и при IERR ≠ 0 значение SF39D полагается равным 1.7D308. |
X = 1.
P = 1.
Y = SF39R (X, P, IERR)
Результаты:
Y = 0.632120558812
IERR = 0