[C] Гамма - функция

[Sergey K]

потребовалось, к сожалению, перевести свои вычисления с Wolfram M. на более быстрый С

Однако, столкнулся с проблемой.

Использую ненормализированную неполную гамма функцию, то что в математике выглядит как

Gamma [a,z]

в библиотеке gsl, которую решил использовать, вызываю аналогичную (по описанию в хелпе):

 double gsl_sf_gamma_inc (double a, double x)

при попытке посчитать, ответы у них разные, причем, серьезно.

возможно, кто-то уже сталкивался с данной проблемой? в чем разница между этими реализациями?

[Vince Diesel]

C си реализациями не сталкивался. А компилирование в самой математике не устраивает? Это порядочно ускоряет процесс. В частности, гамма-функция там компилируется. Можно подключать для этого сторонний компилятор СИ.

[Sergey K]

по мимо скорости, там к сожалению проблема в том, что одно из разросшихся финальных выражений - подсчет интеграла численно - а в него входит эта гамма функция - она не хочет считать вообще никак(

[Vince Diesel]

Можно сравнить с Wolfram Alpha. Например, она выдает

Gamma[2,3]=4/E^3=0.1991482734714558

В математике то же самое. Это правильный ответ — интеграл по частям берется.

[Евгений Машеров]

А там верхняя или нижняя неполная гамма?