потребовалось, к сожалению, перевести свои вычисления с Wolfram M. на более быстрый С
Однако, столкнулся с проблемой.
Использую ненормализированную неполную гамма функцию, то что в математике выглядит как
Код:
Gamma [a,z]
в библиотеке gsl, которую решил использовать, вызываю аналогичную (по описанию в хелпе):
Код:
double gsl_sf_gamma_inc (double a, double x)
при попытке посчитать, ответы у них разные, причем, серьезно.
возможно, кто-то уже сталкивался с данной проблемой? в чем разница между этими реализациями?
