private static double BesselJCoeff(int m, double alpha)
 {
     return(Math.Pow(-1, m) / (MathFunctions.factorial(m) * SpecialFunctions.lanchosGamma(m + alpha + 1)));
 }