Текст подпрограммы и версий sf36r_c.zip sf36d_c.zip |
Тексты тестовых примеров tsf36r_c.zip tsf36d_c.zip |
Вычисление бета - функции распределения вероятностей 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