|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.srcc.msu.ru/num_anal/lib_na/cat/sf/sf86r.htm
Дата изменения: Mon Mar 14 11:07:51 2011 Дата индексирования: Mon Feb 4 02:21:37 2013 Кодировка: Windows-1251 |
|
Текст подпрограммы и версий ( Фортран ) sf86r.zip |
Тексты тестовых примеров ( Фортран ) tsf86r.zip |
|
Текст подпрограммы и версий ( Си ) sf86r_c.zip |
Тексты тестовых примеров ( Си ) tsf86r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf86r_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf86r_p.zip |
Вычисление функции распределения вероятностей Эрланга (распределения Пирсона 3 - го типа)
Подпрограмма SF86R вычисляет функцию распределения Эрланга (распределения Пирсона 3 - го типа) для i - го интервала группировки [xi - 1, xi] ( i = 1, 2, ..., N) при заданном x0 = x0 по формуле
x i x i-1
P = т f(x) dx - т f(x) dx ,
0 0
где f(x) = [ (x - x0)k - 1 / (k - 1)! ] (k/m)k e - k/m (x - x0) ;
m - среднее значение распределения Эрланга;
k = (m - x0) 2 / vars;
x0 - нижняя граница первого интервала группировки;
vars - дисперсия распределения Эрланга;
N - количество интервалов группировки;
xi = xi ;
xi-1 = xi-1
Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979
SUBROUTINE SF86R (RM, VARS, N, X, X0, P)
Параметры
| RM - | заданное среднее значение M распределения Эрланга (тип: вещественный); |
| VARS - | заданная дисперсия vars распределения Эрланга (тип: вещественный); |
| N - | заданное количество интервалов группировки (тип: целый); |
| X - | вещественный одномерный массив длины N, содержащий заданные границы xi интервалов группировки, i = 1, 2, ..., N; |
| X0 - | заданная нижняя граница первого интервала группировки (тип: вещественный); |
| P - |
вещественный одномерный массив длины N + 1, содержащий вычисленные значения
pi (i = 1, 2, ..., N) функции
распределения вероятностей Эрланга для i - го интервала группировки, причем
N
P(N+1) = 1 - е pi ;
i=1
|
Версии: нет
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION X(10), P(11)
DATA X /0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95/
RM = 2.0
VARS = 1.0
N = 10
X0 = 0.0
CALL SF86R (RM, VARS, N, X, X0, P)
Результаты:
P = ( 0.3841996E-5, 0.2620249E-3, 0.1485732E-2,
0.4001807E-2, 0.7705259E-2, 0.1228324E-1,
0.1735353E-1, 0.2254702E-1, 0.2754702E-1,
0.3210832E-1, 0.8747022 )