|
Текст подпрограммы и версий ( Фортран ) sf60r.zip sf61r.zip |
Тексты тестовых примеров ( Фортран ) tsf60r.zip tsf61r.zip |
|
Текст подпрограммы и версий ( Си ) sf60r_c.zip sf61r_c.zip |
Тексты тестовых примеров ( Си ) tsf60r_c.zip tsf61r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf60r_p.zip , sf61r_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf60r_p.zip , tsf61r_p.zip |
Вычисление функции распределения вероятностей Пуассона.
SF60R вычисляет вероятности того, что целочисленная случайная величина, распределенная по закону Пуассона с параметром RLАМ, примет значение I
P {ξ = I} = T(I + 1) = (RLAM) I exp (- RLAM) / I ! , I = 0, 1, ..., K
а также вероятность P = P {ξ ≤ K} , т.е.
K
P = ∑ T(I + 1)
I=0
M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, New York, Dover Publications, Inc., 1965.
SUBROUTINE SF60R (K, RLAM, T, P)
Параметры
| K - | максимальное значение I, для которого вычисляется вероятность P {ξ = I} (тип: целый); |
| RLAM - | заданный параметр распределения вероятностей Пуассона (тип: вещественный); |
| T - | вещественный вектоp длины K + 1, содержащий вычисленные значения вероятностей P {ξ = I} , 0 ≤ I ≤ K; |
| P - | вещественная переменная, содержащая вычисленное значение вероятности P {ξ ≤ K}. |
Версии
| SF61R - | вычисление функции распределения вероятностей Пуассона P {ξ ≤ K} (см. замечания по использованию). |
Вызываемые подпрограммы: нет
Замечания по использованию
|
Если P {ξ = I} меньше наименьшего представимого на данной ЭВМ вещественного числа , то программа полагает T(I) = 0. Если RLAM ≤ 0 , то P {ξ ≤ K} = 1. , T(1) = 1. , I(J) = 0, 2 ≤ J ≤ K. Подпрограмма SF61R имеет следующий список формальных параметровSUBROUTINE SF61R (K, RLAM, P) , так как производится вычисление только вероятности P {ξ ≤ K}. Смысл остальных параметров остается прежним. |
REAL T(2)
K = 1
RLAM = 0.6
CALL SF60R (K, RLAM, T, P)
Результаты:
T(1) = 0.54881163609
T(2) = 0.32928698166
P = 0.87809861775