|
Текст подпрограммы и версий ( Фортран ) sf22r.zip sf22d.zip |
Тексты тестовых примеров ( Фортран ) tsf22r.zip tsf22d.zip |
|
Текст подпрограммы и версий ( Си ) sf22r_c.zip sf22d_c.zip |
Тексты тестовых примеров ( Си ) tsf22r_c.zip tsf22d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf22r_p.zip , sf22e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf22r_p.zip , tsf22e_p.zip |
Вычисление значений бета-функции.
Функция SF22R вычисляет значение интегралов вида:
1
B (x, y) = ∫ t x - 1(1 - t)y - 1 dt
0
Бета-функция может быть представлена с помощью гама-функции
B (x, y) = Г(x)Г(y) / Г(x + y).
J.F.Hart, E.W.Cheney, C.L.Lawson et al, Computer Approximations, Wiley, New York, 1968.
FUNCTION SF22R (X, Y, IERR)
Параметры
| X, Y - | заданные значения аргументов x и y соответственно (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда значение одного из аргументов (X, Y, X + Y) pавно нулю или целому отрицательному числу; значение функции полагается равным 3.4E38; |
| IERR=66 - | когда значение одного из аргументов (X, Y или X + Y) меньше чем - 20.4; значение функции полагается равным 3.4E38. |
Версии
| SF22D - | вычисление значений бета - функции с повышенной точностью. |
Вызываемые подпрограммы
| UTSF10 - | подпрограмма выдачи диагностических сообщений при работе функции SF22R. |
| UTSF11 - | подпрограмма выдачи диагностических сообщений при работе функции SF22D. |
Замечания по использованию
| Областью допустимых значений аргументов для значений аргументов SF22R является луч z ≥ - 20.4, z ≠ 0, - 1, - 2, ..., - 20, а для функции SF22D - z ≥ - 961.199, z ≠ 0, - 1, - 2, ..., - 961. Kpоме того, для функции SF22D аргументы X и Y имеют тип DOUBLE PRECISION и при нарушении области допустимых значений функция полагается равной 1.7D308. |
X = 1.5
Y = 2.9
Y1 = SF22R (X, Y, IERR)
Результаты: Y1 = 0.159770619769 , IERR = 0