|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.srcc.msu.ru/num_anal/lib_na/cat/sf/sf36r.htm
Дата изменения: Mon Mar 14 11:07:47 2011 Дата индексирования: Mon Feb 4 02:18:33 2013 Кодировка: Windows-1251 |
|
Текст подпрограммы и версий ( Фортран ) sf36r.zip sf36d.zip |
Тексты тестовых примеров ( Фортран ) tsf36r.zip tsf36d.zip |
|
Текст подпрограммы и версий ( Си ) sf36r_c.zip sf36d_c.zip |
Тексты тестовых примеров ( Си ) tsf36r_c.zip tsf36d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf36r_p.zip , sf36e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf36r_p.zip , tsf36e_p.zip |
Вычисление бета - функции распределения вероятностей Ix (a, b).
SF36R вычисляет вероятность того, что случайная величина, подчиненная закону бета распределения с параметрами a и b, меньше или pавна значению x.
Другими словами, SF36R вычисляет отношение
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.
REALL FUNCTION SF36R (X, A, B, IERR)
Параметры
| X - | заданное значение аргумента x (тип: вещественный); |
| A, B - | заданные значения параметров бета распределения (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда значение аргумента X не принадлежит отрезку (0, 1); значение функции полагается равным 3.4E38; |
| IERR=66 - | когда один или оба параметра бета pаспределения меньше либо равны 0; значение функции полагается равным 3.4E38. |
Версии
| SF36D - | вычисление бета - функции распределения вероятностей Ix (a, b) с повышенной точностью. |
Вызываемые подпрограммы
| UTSF12 - | подпрограмма выдачи диагностических сообщений при работе функции SF36R. |
| UTSF13 - | подпрограмма выдачи диагностических сообщений при работе функции SF36D. |
Замечания по использованию
| Для функции SF36D параметры X, A и B должны иметь тип DOUBLE PRECISION и при IERR ? 0 значение SF36D полагается равным 1.7D308. |
X = 0.99
A = 5.5
B = 0.5
Y = SF36R (X, A, B, IERR)
Результаты:
Y = 0.7451499080, IERR = 0