Функция: sf37r_c
Назначение
Вычисление неполной бета - функции
Bx (a, b) вещественных аргументов.
Математическое описание
sf37r_c вычисляет неполную бета - функцию вещественных
аргументов a и b (a > 0, b > 0)
x
Bx (a, b) = ∫ t a - 1 (1 - t) b - 1 dt (0 ≤ x ≤ 1)
0
M.Abramowitz, I.A.Stegun, Handbook of Mathematical
Functions, New York, Dover Publications, INC., 1965.
Использование
real sf37r_c (real *x, real *a, real *b, integer *ierr)
Параметры
x -
|
заданное значение верхнего предела
интегрирования x (тип: вещественный);
|
a, b -
|
заданные значения аргументов a и b (тип: вещественный);
|
ierr -
|
целая переменная, служащая для сообщения об
ошибках, обнаруженных в ходе работы функции; при этом:
|
ierr=65 -
|
когда значение аргумента x не
придадлежит отрезку (0, 1); значение функции
полагается равным 3.4e38;
|
ierr=66 -
|
когда один или оба параметра неполной
бета - функции меньше либо равны 0;
значение функции полагается равным 3.4e38.
|
Версии
sf37d_c -
|
вычисление неполной бета - функции
Bx (a, b)
вещественных аргументов с повышенной точностью.
|
Вызываемые подпрограммы
sf22r_c -
|
подпрограмма вычисления значений бета - функции.
|
sf36r_c -
|
подпрограмма вычисления бета - функции
распределения вероятностей Ix (a, b).
|
sf22d_c -
|
подпрограмма вычисления значений бета - функции
с повышенной точностью.
|
sf36d_c -
|
подпрограмма вычисления бета - функции
распределения вероятностей Ix (a, b)
с повышенной точностью.
|
utsf12_c -
|
подпрограмма выдачи диагностических сообщений
при работе функциии sf37r_c.
|
utsf13_c -
|
подпрограмма выдачи диагностических сообщений
при работе функции sf37d_c.
|
Замечания по использованию
|
Для функции sf37d_c параметры x, a и b должны иметь тип
double и при
ierr ≠ 0 значение sf37d_c
полагается равным 1.7e308.
|
Пример использования
int main(void)
{
/* Local variables */
static int ierr;
extern float sf37r_c(float *, float *, float *, int *);
static float a, b, x, y;
x = 1.f;
a = 2.f;
b = 3.f;
y = (float)sf37r_c(&x, &a, &b, &ierr);
printf("\n %16.7e %5i \n",y,ierr);
return 0;
} /* main */
Результаты:
y = 0.0833333333,
ierr = 0