Текст подпрограммы и версий sf79r_c.zip |
Тексты тестовых примеров tsf79r_c.zip |
Вычисление обратной функции распределения вероятностей Коши.
sf79r_c вычисляет значение x такое, что F (X, λ ,A) = P, где F (Х, λ, А) - функция распределения Коши:
x F(X, λ, A) = 1/πλ ∫ dz / (λ2 + (A - Z)2) , λ > 0 -∞
Справочник по специальным функциям. Под ред. M.Абрамовица, И.Стеган. M.: "Hаука", 1979.
int sf79r_c (real *p, real *rl, real *a, real *x, integer *ierr)
Параметры
p - | заданное значение аргумента (тип: вещественный); |
rl, a - | параметры λ, A распределения Коши (тип: вещественный); |
x - | вещественная переменная, содержащая вычисленное значение функции; |
ierr - | целая переменная, служащая для сообщения об ошибках, при этом: |
ierr=65 - | когда λ ≤ 0; |
ierr=66 - | когда p < 0 либо p > 1; значение x в обоих случаях полагается равным нулю. |
Версии: нет
Вызываемые подпрограммы
utsf16_c - | подпрограмма выдачи диагностических сообщений при работе подпрограммы sf79r_c. |
Замечания по использованию: нет
int main(void) { /* Local variables */ extern int sf79r_c(float *, float *, float *, float *, int *); static float x; static int ierr; static float c_b1 = .5f; static float c_b2 = 1.f; sf79r_c(&c_b1, &c_b2, &c_b2, &x, &ierr); printf("\n %16.7e %5i \n",x,ierr); return 0; } /* main */ Результаты: x = 1.000 000 000 005 , ierr = 0