%I #24 Jun 20 2021 02:46:12
%S 1,8,5,4,1,0,5,9,6,7,9,2,1,0,2,6,4,3,2,7,4,8,3,7,0,7,1,8,4,1,0,2,9,3,
%T 2,4,5,4,2,9,2,3,2,6,7,5,0,2,7,2,6,1,9,3,0,8,4,6,9,7,5,1,0,8,4,6,8,8,
%U 0,6,2,1,2,4,8,7,3,2,6,1,6,6,5,5,9,2,4,0,3,3,6,6,1,7,0,6,8,2,4,3,8,8,0
%N Decimal expansion of x such that x^x = Pi.
%H G. C. Greubel, <a href="/A030437/b030437.txt">Table of n, a(n) for n = 1..10000</a>
%e x = 1.8541059679210264327483707184102932454292... .
%p x^x=Pi; solve(%,x); evalf(%, 140); # solution is log(Pi)/LambertW(log(Pi)), where LambertW is the Omega function.
%t x=Pi; RealDigits[Log[x]/ProductLog[Log[x]],10,6! ][[1]] (* _Vladimir Joseph Stephan Orlovsky_, Feb 11 2010 *)
%o (PARI) solve(x=1, 2, x^x-Pi) \\ _Michel Marcus_, Jan 14 2015
%o (PARI) exp(lambertw(log(Pi))) \\ _Charles R Greathouse IV_, Nov 11 2017
%Y Cf. A000796 (Pi), A100947 (continued fraction), A073243 (reciprocal).
%K nonn,cons
%O 1,2
%A James L. Dean (csvcjld(AT)nomvs.lsumc.edu)
%E More terms from _Simon Plouffe_
%E Better name from _Jon E. Schoenfield_, Dec 30 2014
|