Текст подпрограммы и версий ( Фортран )
sf37r.zip  sf37d.zip 
Тексты тестовых примеров ( Фортран )
tsf37r.zip  tsf37d.zip 
Текст подпрограммы и версий ( Си )
sf37r_c.zip  sf37d_c.zip 
Тексты тестовых примеров ( Си )
tsf37r_c.zip  tsf37d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf37r_p.zip , sf37e_p.zip
Тексты тестовых примеров ( Паскаль )
tsf37r_p.zip , tsf37e_p.zip

Функция:  SF37R

Назначение

Вычисление неполной бета - функции  Bx (a, b) вещественных аргументов.

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

SF37R вычисляет неполную бета - функцию вещественных аргументов 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  FUNCTION  SF37R (X, A, B, IERR) 

Параметры

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

Версии

SF37D - вычисление неполной бета - функции  Bx (A, B) вещественных аргументов с повышенной точностью.

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

SF22R - подпрограмма вычисления значений бета - функции.
SF36R - подпрограмма вычисления бета - функции распределения вероятностей  Ix (A, B).
SF22D - подпрограмма вычисления значений бета - функции с повышенной точностью.
SF36D - подпрограмма вычисления бета - функции распределения вероятностей  Ix (A, B) с повышенной точностью.
UTSF12 - подпрограмма выдачи диагностических сообщений при работе функциии SF37R.
UTSF13 - подпрограмма выдачи диагностических сообщений при работе функции SF37D.

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

  Для функции SF37D параметры X, A и B должны иметь тип DOUBLE PRECISION и при IERR ≠ 0 значение SF37D полагается равным 1.7D308.

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

         X = 1.
         A = 2.
         B = 3.
         Y = SF37R (X, A, B, IERR)

Результаты:

       Y  =  0.0833333333, 
       IERR  =  0