#15033: Wrong limit value of expression involving gamma function
-------------------------+-------------------------------------------------
Reporter: JGuzman | Owner:
Type: defect | Status: new
Priority: major | Milestone: sage-5.12
Component: calculus | Keywords: maxima, gamma, limit
Merged in: | Authors:
Reviewers: | Report Upstream: Reported upstream. No
Work issues: | feedback yet.
Dependencies: | Branch:
| Stopgaps:
-------------------------+-------------------------------------------------
When computing the limit of the following function:
{{{
sage: f(x) = gamma(x+1/2)/gamma(x)/sqrt(x)
sage: limit(f,x=oo)
}}}
'''Sage returns 0''' in stead of '''one'''. See for example:
{{{
sage: plot(y,x,1,50)
}}}
or how the value approaches to one:
{{{
sage : for x in range(1, 20, 2):
sage : print("f(%2d)=%f" %(i, f(i).n() ) )
}}}
{{{
f( 1)=0.886227
f( 3)=0.959369
f( 5)=0.975350
f( 7)=0.982316
f( 9)=0.986214
f(11)=0.988705
f(13)=0.990433
f(15)=0.991703
f(17)=0.992675
f(19)=0.993443
}}}
As described in [https://groups.google.com/forum/#!topic/sage-
support/y9UnyWbUagY sage-support mailing list], this is probably a bug in
Maxima, and it has been reported
[http://sourceforge.net/p/maxima/bugs/2621/ here]
Tested in sage-5.8
--
Ticket URL: <http://trac.sagemath.org/ticket/15033>
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/groups/opt_out.
