Держи через аппроксимацию Ланшольца гамма функции #include
#include
#include
#define EPSILON 1e-07
double g = 7;
double n = 9;
double p[] = {
0.99999999999980993,
676.5203681218851,
-1259.1392167224028,
771.32342877765313,
-176.61502916214059,
12.507343278686905,
-0.13857109526572012,
9.9843695780195716e-6,
1.5056327351493116e-7};
std::complex drop_imag(std::complex z) {
if (abs(z.imag())