#17328: incomplete gamma function bugs for certain arguments
-------------------------------------+-------------------------------------
Reporter: kcrisman | Owner:
Type: defect | Status: needs_review
Priority: major | Milestone: sage-6.4
Component: calculus | Resolution:
Keywords: incomplete gamma | Merged in:
function | Reviewers:
Authors: Peter Bruin | Work issues:
Report Upstream: N/A | Commit:
Branch: | 3cae7c8469d6c2f5dd3d15a872965f528f637c5c
u/pbruin/17328-incgam | Stopgaps:
Dependencies: |
-------------------------------------+-------------------------------------
Comment (by kcrisman):
Seems to give correct answers in various precisions.
{{{
plot(lambda t: gamma(-1, -log(t)).real() -
numerical_integral(1/ln(x)^2,2,t)[0],2,3,color='red')
}}}
is a very nice straight line ''just'' below zero now. And the imaginary
part seems to be consistently pi for t>1. And the fix is a good one.
> We should also add a doctest for the original bug.
Did you want to add
{{{
N(integral(1/log(x)^2,(x,2,3)))
}}}
as well then, or is the example you have with the actual antideriv and
floats sufficient?
--
Ticket URL: <http://trac.sagemath.org/ticket/17328#comment:8>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" 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-trac.
For more options, visit https://groups.google.com/d/optout.
