|
Текст подпрограммы и версий sf43r_c.zip sf43d_c.zip |
Тексты тестовых примеров tsf43r_c.zip tsf43d_c.zip |
Вычисление функции F - распределения вероятностей (распределение Снедекора) с вещественными степенями свободы.
sf43r_c вычисляет вероятность того, что случайная величина, подчиненная закону F - распределения вероятностей (распределению Снедекора) с вещественными степенями свободы, меньше или равна значению f .
M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, New York, Dover Publications, INC., 1965.
real sf43r_c (real *f, real *an1, real *an2, integer *ierr)
Параметры
| f - | заданное значение аргумента f (тип: вещественный); |
|
an1 - an2 | заданные значения степеней свободы (тип: вещественный); |
| ierr - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| ierr=65 - | когда одна или обе степени свободы меньше или равны 0; значения функции полагаются равными 3.4e38; |
| ierr=66 - | когда заданный предел интегрирования меньше 0; значение функции полагается равным 3.4e38. |
Версии
| sf43d_c - | вычисление функции F - распределения вероятностей (распределение Снедекора) с вещественными степенями свободы с повышенной точностью. |
Вызываемые подпрограммы
| sf36r_c - | вычисление бета - функции распределения вероятностей Ix (p, q) (используется в sf43r_c). |
| sf36d_c - | вычисление бета - функции распределения вероятностей Ix (p, q) с повышенной точностью (используется в sf43d_c). |
| utsf12_c - | подпрограмма выдачи диагностических сообщений при работе функции sf43r_c. |
| utsf13_c - | подпрограмма выдачи диагностических сообщений при работе функции sf43d_c. |
Замечания по использованию
| В sf43d_c параметры f, an1 и an2 должны иметь тип double и при ierr ≠ 0 значение sf43d_c полагается равным 1.7e308. |
int main(void)
{
/* Local variables */
extern float sf43r_c(float *, float *, float *, int *);
static int ierr;
static float f, y, an1, an2;
f = .44444444444444442f;
an1 = 12.f;
an2 = 8.f;
y = (float)sf43r_c(&f, &an1, &an2, &ierr);
printf("\n %16.7e %5i \n",y,ierr);
return 0;
} /* main */
Результаты:
y = 0.099352575999
ierr = 0