Текст подпрограммы и версий
sf48r_c.zip  		Тексты тестовых примеров
tsf48r_c.zip 

Функция:  sf48r_c

Назначение

Вычисление обратной функции χ2 - распределения вероятностей.

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

sf48r_c вычисляет значение  X такое, что 

          χ2 (X; AN) = P ,

где  P   - заданное значение вероятности (0 < P < 1) ,
      AN - заданное значение степени свободы χ2 - распределения
               вероятностей (0.5 < AN < 200000).

M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, Dover Publications, INC., 1965.

Использование 

    real sf48r_c (real *p, real *an, integer *ierr)

   Параметры
p - заданное значение вероятности (тип: вещественный);
an - заданное значение степени свободы (тип: вещественный);
ierr - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
ierr=65 - когда значение вероятности  p  не принадлежит интервалу (0, 1); значение функции полагается равным 3.4e38;
ierr=66 - когда заданное значение степени свободы меньше 0.5; значение функции полагается равным 3.4e38. 
   Версии:  нет 
   Вызываемые подпрограммы
sf47r_c - вычисление функции χ2 - распределения вероятностей с вещественными степенями свободы;
sf51r_c - вычисление обратной функции нормального pаспределения вероятностей;
utsf12_c - подпрограмма выдачи диагностических сообщений в ходе работы функции sf48r_c. 
   Замечания по использованию:  нет

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

int main(void)
{
    /* Local variables */
    static int ierr;
    extern float sf48r_c(float *, float *, int *);
    static float p, y, an;

    an = 2.f;
    p = .99f;
    y = (float)sf48r_c(&p, &an, &ierr);

    printf("\n %16.7e %16.7e \n",p,an);
    printf("\n %16.7e %5i \n",y,ierr);
    return 0;
} /* main */


Результаты:

       y  =  9.210340350925
       ierr  =  0