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

Подпрограмма:  sf81r_c

Назначение

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

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

Подпрограмма sf81r_c вычисляет значение  x, для которого вероятность того, что случайная величина, имеющая распределение Парето с параметрами  A,  x0, меньше или pавна  x, pавно заданному значению  P 

                         x
        P = A/x0  ∫ (x0 / Z)A - 1 dZ  ,      x ≥ x0 > 0  ,   A > 0
                       x0

Справочник по специальным функциям. Под ред. M.Абрамовица, И.Стеган. M.: "Hаука", 1979.

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

    int sf81r_c (real *p, real *a, real *x0, real *x,
            integer *ierr)

   Параметры
p - заданное значение аргумента обратной функции распределения (тип: вещественный);
a, x0 - заданные значения параметров распределения Парето (тип: вещественный);
x - вещественная переменная, содержащая вычисленное значение обратной функции распределения Парето;
ierr - целая переменная, служащая для сообщения об ошибках; при этом:
ierr=65 - когда  a ≤ 0 либо  x0 ≤ 0,
ierr=66 - когда  p < 0 либо  p > 1; в обоих случаях значение  x полагается равным  x0. 
   Версии:  нет 
   Вызываемые подпрограммы
utsf16_c - подпрограмма выдачи диагностических сообщений при работе подпрограммы sf80r_c. 
   Замечания по использованию:  нет

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

int main(void)
{
    /* Local variables */
    extern int sf81r_c(float *, float *, float *, float *, int *);
    static float x;
    static int ierr;
    static float c_b1 = .75f;
    static float c_b2 = 1.f;

    sf81r_c(&c_b1, &c_b2, &c_b2, &x, &ierr);

    printf("\n %16.7e %5i \n",fx,ier);
    return 0;
} /* main */


Результаты:    fx  =  4.0 ,   ier  =  0