|
%I #47 Jan 08 2024 17:01:05
%S 8,8,5,6,0,3,1,9,4,4,1,0,8,8,8,7,0,0,2,7,8,8,1,5,9,0,0,5,8,2,5,8,8,7,
%T 3,3,2,0,7,9,5,1,5,3,3,6,6,9,9,0,3,4,4,8,8,7,1,2,0,0,1,6,5,8,7,5,1,3,
%U 6,2,2,7,4,1,7,3,9,6,3,4,6,6,6,4,7,9,8,2,8,0,2,1,4,2,0,3,5,9
%N Decimal expansion of real number y such that y = Gamma(x) is a minimum.
%H G. C. Greubel, <a href="/A030171/b030171.txt">Table of n, a(n) for n = 0..5000</a>
%H Simon Plouffe, <a href="http://www.plouffe.fr/simon/constants/minygam.txt">Minimal point of GAMMA(x)</a>
%H Simon Plouffe, <a href="http://www.worldwideschool.org/library/books/sci/math/MiscellaneousMathematicalConstants/chap68.html">Minimal y of GAMMA(x)</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/GammaFunction.html">Gamma Function</a>
%e 0.885603194410888700278815900582588733207951533669903448871200165875136...
%p Digits:=500; x0:=fsolve(Psi(x)=0, x); evalf(GAMMA(x0), 120) # _Iaroslav V. Blagouchine_, Feb 16 2016
%t First@ RealDigits[ FindMinimum[ Gamma[x], {x, 1.4}, WorkingPrecision -> 2^7][[1]]] (* _Robert G. Wilson v_, Aug 03 2010 *)
%t RealDigits[ Gamma[x /. FindRoot[ PolyGamma[0, x] == 0, {x, 1}, WorkingPrecision -> 100]]][[1]] (* _Jean-François Alcover_, Oct 23 2012 *)
%o (PARI) gamma(solve(x=1,2,psi(x))) \\ _Charles R Greathouse IV_, Apr 17 2015
%Y Cf. A030169 for value of x.
%K nonn,cons
%O 0,1
%A _Eric W. Weisstein_
|
