Текст подпрограммы и версий ( Фортран )
sf86r.zip  		Тексты тестовых примеров ( Фортран )
tsf86r.zip 
Текст подпрограммы и версий ( Си )
sf86r_c.zip  		Тексты тестовых примеров ( Си )
tsf86r_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf86r_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf86r_p.zip

Подпрограмма:  SF86R

Назначение

Вычисление функции распределения вероятностей Эрланга (распределения Пирсона 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   )