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 sage-support mailing list, this is probably a bug in Maxima, and it has been reported here

Tested in sage-5.8

Depends on #13973
Depends on #16224

Upstream: Fixed upstream, but not in a stable release.

CC: @sagetrac-JGuzman

Component: calculus

Keywords: maxima, gamma, limit

Author: Peter Bruin

Branch/Commit: cab4585

Reviewer: Karl-Dieter Crisman

Issue created by migration from https://trac.sagemath.org/ticket/15033