Here too -- would it be better to use the exact gamma(2n,z) function, rather than the asymptotic form, since you have it now?  We don't actually need this to be quick, since it's just in the SersicInfo startup function, so the time normally gets amortized over many uses.  This would probably help us be more confident in increasing the range of n values (in #325).