|
|
A242813
|
|
Decimal expansion of the expected number of returns to the origin of a random walk on a 5-d lattice.
|
|
7
|
|
|
1, 1, 5, 6, 3, 0, 8, 1, 2, 4, 8, 4, 0, 2, 3, 1, 1, 7, 8, 7, 0, 7, 1, 3, 5, 1, 2, 1, 9, 3, 8, 5, 6, 6, 9, 8, 5, 5, 4, 5, 4, 2, 7, 3, 4, 8, 5, 0, 5, 1, 4, 2, 3, 8, 8, 2, 6, 9, 5, 6, 6, 0, 1, 1, 2, 1, 0, 0, 8, 7, 7, 0, 3, 4, 7, 0, 6, 8, 7, 3, 1, 1, 7, 2, 3, 6, 6, 5, 4, 3, 0, 4, 9, 5, 0, 9, 1, 7, 1, 6, 5, 2, 6, 7, 4, 3
(list;
constant;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
REFERENCES
|
Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 5.9 Polya's random walk constants, p. 323.
|
|
LINKS
|
|
|
FORMULA
|
m(d) = d/(2*Pi)^d*multipleIntegral(-Pi..Pi) (d-sum_(k=1..d) cos(t_k))^(-1) dt_1 dt_2 ... dt_d, where d is the lattice dimension.
m(d) = integral_(t>0) exp(-t)*BesselI(0,t/d)^d dt where BesselI(0,x) is the zeroth modified Bessel function.
|
|
EXAMPLE
|
1.1563081248...
|
|
MAPLE
|
m5:= int(exp(-t)*BesselI(0, t/5)^5, t=0..infinity):
s:= convert(evalf(m5, 120), string):
map(parse, subs("."=NULL, [seq(i, i=s)]))[]; # Alois P. Heinz, May 23 2014
|
|
MATHEMATICA
|
d = 5; d/Pi^d*NIntegrate[(d - Sum[Cos[t[k]], {k, 1, d}])^-1, Sequence @@ Table[{t[k], 0, Pi}, {k, 1, d}] // Evaluate] // RealDigits[#, 10, 10]& // First
|
|
PROG
|
|
|
CROSSREFS
|
Cf. A086230, A086231, A086232, A086233, A086234, A086235, A086236, A242812, A242814, A242815, A242816.
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|