Документ взят из кэша поисковой машины. Адрес
оригинального документа
: 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 )