public static double ComplementedPoissonDistribution(int k, double x) { if (k < 0 || x < 0) return 0.0; return IncompleteGamma((double) (k + 1), x); }