OFFSET
1,1
COMMENTS
Solutions were provided by Bob Hanlon and Peter Pein.
EXAMPLE
4.331292424499713465838941491042338081138561546026782297287499557485752...
MAPLE
read("transforms3") ; Digits := 120 ; x := 4.0 ; for l from 1 to 10 do x := x-(1-40/GAMMA(x+1))/Psi(x+1) ; x := evalf(x) ; end do; CONSTTOLIST(x) ; # R. J. Mathar, Mar 23 2010
MATHEMATICA
FindRoot[(x + 1)! == 40, {x, 4, 5}, WorkingPrecision -> 50]( * Bob Hanlon (hanlonr(AT)cox.net) and Peter Pein, edited by M. F. Hasler, Jun 16 2024 *)
PROG
(PARI) A129624_first(N=99, t=40, L=log(t)*.55)=localprec(N+5); digits(solve(x=L+1.5, L+2.5, gamma(x+1)-t)\10^-N) \\ M. F. Hasler, Jun 13 2024
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Roger L. Bagula, May 30 2007
EXTENSIONS
More digits from R. J. Mathar, Mar 23 2010
STATUS
approved