When using the following values...
double a = 0;
double x = 0.2247;
double ans = MathNet.Numerics.SpecialFunctions.GammaUpperIncomplete(a, x);
...I am getting an answer of infinity.
...gives me the correct answer of 1.128
When I numerically integrate I also get an answer of approx. 1.128
What am I doing wrong when using your incomplete gamma function? The problem centers around using a value for a of zero.
Thanks! Love your math library!