Текст подпрограммы и версий
sf36r_c.zip  sf36d_c.zip 
Тексты тестовых примеров
tsf36r_c.zip  tsf36d_c.zip 

Функция:  sf36r_c

Назначение

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

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

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

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

      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.

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

    real sf36r_c (real *x, real *a, real *b, integer *ierr)

Параметры

x - заданное значение аргумента  x (тип: вещественный);
a, b - заданные значения параметров бета распределения (тип: вещественный);
ierr - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
ierr=65 - когда значение аргумента x не принадлежит отрезку (0, 1); значение функции полагается равным 3.4e38;
ierr=66 - когда один или оба параметра бета pаспределения меньше либо равны 0; значение функции полагается равным 3.4e38.

Версии

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

Вызываемые подпрограммы

utsf12_c - подпрограмма выдачи диагностических сообщений при работе функции sf36r_c.
utsf13_c - подпрограмма выдачи диагностических сообщений при работе функции sf36d_c.

Замечания по использованию

  Для функции sf36d_c параметры x, a и b должны иметь тип double и при  ierr ≠ 0 значение sf36d_c полагается равным 1.7e308.

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

int main(void)
{
    /* Local variables */
    static int ierr;
    extern float sf36r_c(float *, float *, float *, int *);
    static float a, b, x, y;

    x = .99f;
    a = 5.5f;
    b = .5f;
    y = (float)sf36r_c(&x, &a, &b, &ierr);

    printf("\n %16.7e %16.7e %16.7e \n",x,a,b);
    printf("\n %16.7e %5i \n",y,ierr);
    return 0;
} /* main */


Результаты:

       y  =  0.7451499080,     ierr  =  0