Текст подпрограммы и версий ( Фортран )
sf36r.zip  sf36d.zip  		Тексты тестовых примеров ( Фортран )
tsf36r.zip  tsf36d.zip 
Текст подпрограммы и версий ( Си )
sf36r_c.zip  sf36d_c.zip  		Тексты тестовых примеров ( Си )
tsf36r_c.zip  tsf36d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf36r_p.zip , sf36e_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf36r_p.zip , tsf36e_p.zip

Функция:  SF36R

Назначение

Вычисление бета - функции распределения вероятностей   Ix (a, b).

Математическое описание

SF36R вычисляет вероятность того, что случайная величина, подчиненная закону бета распределения с параметрами a и b, меньше или pавна значению  x.

Другими словами, SF36R вычисляет отношение 

      Ix (a, b) = Bx (a, b) / B (a, b) =
                 x                                                                                                             
                ∫ t a - 1 (1 - t) b - 1 dt 
               0                                 
          =  ---------------------------      ,
                 1                                                                            
                ∫ t a - 1 (1 - t) b - 1 dt
               0

где Bx (a, b) - неполная бета - функция,
      B (a, b) - бета - функция (полная),
      a и b - принимают вещественные значения (a > 0, b > 0).

M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, New York, Dover Publications, INC., 1965.

Использование 

    REALL  FUNCTION  SF36R (X, A, B, IERR) 
   Параметры
X - заданное значение аргумента  x (тип: вещественный);
A, B - заданные значения параметров бета распределения (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=65 - когда значение аргумента X не принадлежит отрезку (0, 1); значение функции полагается равным 3.4E38;
IERR=66 - когда один или оба параметра бета pаспределения меньше либо равны 0; значение функции полагается равным 3.4E38. 
   Версии
SF36D - вычисление бета - функции распределения вероятностей  Ix (a, b) с повышенной точностью. 
   Вызываемые подпрограммы
UTSF12 - подпрограмма выдачи диагностических сообщений при работе функции SF36R.
UTSF13 - подпрограмма выдачи диагностических сообщений при работе функции SF36D. 
   Замечания по использованию
  Для функции SF36D параметры X, A и B должны иметь тип DOUBLE PRECISION и при  IERR ≠ 0 значение SF36D полагается равным 1.7D308.

Пример использования 

         X = 0.99
         A = 5.5
         B = 0.5
         Y = SF36R (X, A, B, IERR)

Результаты:

       Y  =  0.7451499080,     IERR  =  0