|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Evgenij Masherov 2:5020/175.2 05 May 2003 15:37:35
To : Vladislav Terehov
Subject : Апроксимация функций: ди/поли гаммы, errf
--------------------------------------------------------------------------------
Wed Apr 30 2003 00:28, Vladislav Terehov wrote to All:
VT> Hадобны сабжевые алгоpитмы (пpименительно к полю вещественных чисел),
VT> лучше, если они будут офоpмлены в виде кода на c.
VT> (errf - функция веpоятности ошибок).
Hу, скажем, дигамма.
Для начала заметим, что
Psi(x)=Psi(x+1)-1/(x+1)
так что аргумент можно нарастить до достаточно большой величины.
Затем берем разложение
Psi(x)=ln(x)+1/(2*x)-1/(12*x^2)+1/(120*x^4)-1/(252*x^6)
Погрешность, меньшая чем 1/(240*x^8)
Для отрицательных
Psi(-x)=Pi*ctg(Pi*x)+Psi(x-1)
Евгений Машеров АКА СанитарЖеня
--- ifmail v.2.15dev5
* Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
Апроксимация функций: ди/поли гаммы, errf 