14 throw itk::ExceptionObject(__FILE__, __LINE__,
"The Psi function is not defined in 0.",ITK_LOCATION);
16 double resVal = - emc;
17 for (
unsigned int i = 1;i < n;++i)
18 resVal += 1.0 / ((
double)i);
25 double resVal = std::tgamma(2*n + 1);
27 resVal /= std::pow(4.0,(
double)n);
28 resVal /= std::tgamma(n+1);
35 double resVal = std::pow(-4.0,(
double)n);
36 resVal *= std::tgamma(n+1);
38 resVal /= std::tgamma(2*n + 1);
double psi_function(unsigned int n, double emc)
double gammaHalfPlusN(unsigned int n)
double gammaHalfMinusN(unsigned int n)