Документ взят из кэша поисковой машины. Адрес оригинального документа : http://num-anal.srcc.msu.su/lib_na/cat/sf_htm_c/sf37r_c.htm
Дата изменения: Tue Apr 28 08:29:42 2015
Дата индексирования: Sun Apr 10 01:45:59 2016
Кодировка: Windows-1251
БЧА НИВЦ МГУ. SF37R_C. Гамма-функция и связанные с ней функции
Текст подпрограммы и версий
sf37r_c.zip  sf37d_c.zip  		Тексты тестовых примеров
tsf37r_c.zip  tsf37d_c.zip 

Функция:  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