On Saturday, May 3, 2014 7:17:04 AM UTC-4, [email protected] wrote: > > I am new to Sage; trying to explore open source alternatives to > Mathematica. > > However, I seem to be having trouble with the incomplete gamma function. > Here are two difficulties. First, in trying to evaluate the incomplete > gamma function at a point where the result should be very small, I just get > zero even if I increase the precision arbitrarily. In particular consider > > numerical_approx(gamma(9, 10^(-3))-gamma(9), digits=40) > > the value of this number is approximate -1.1 \times 10^(-28), but I just > get 0.00000 ...0 for the input above. > > I believe mpmath (in Sage) should also be able to do this fine, though I haven't tested your particular case:
http://mpmath.googlecode.com/svn/trunk/doc/build/functions/expintegrals.html#gammainc -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sage-support. For more options, visit https://groups.google.com/d/optout.
