cp's OEIS Frontend

A253095 Moments of 4-step random walk in 4 dimensions.

%I A253095 #25 Apr 05 2024 10:10:43
%S A253095 1,4,22,148,1144,9784,90346,885868,9115276,97578688,1079676448,
%T A253095 12285725632,143204046496,1704422018992,20660609113186,
%U A253095 254522834851516,3180935346538684,40269426101933392,515743456513546072,6675036087017279056,87221496402779437696,1149701868292524559744
%N A253095 Moments of 4-step random walk in 4 dimensions.
%H A253095 J. M. Borwein, <a href="https://scholarship.claremont.edu/jhm/vol6/iss1/7">A short walk can be beautiful</a>, Journal of Humanistic Mathematics, Volume 6 Issue 1 (January 2016), pages 86-109.
%H A253095 J. M. Borwein, <a href="https://carmamaths.org/resources/jon/OEIStalk.pdf">Adventures with the OEIS: Five sequences Tony may like</a>, Guttmann 70th [Birthday] Meeting, 2015, revised May 2016.
%H A253095 J. M. Borwein, <a href="/A060997/a060997.pdf">Adventures with the OEIS: Five sequences Tony may like</a>, Guttmann 70th [Birthday] Meeting, 2015, revised May 2016. [Cached copy, with permission]
%H A253095 Jonathan M. Borwein, Armin Straub and Christophe Vignat, <a href="https://arminstraub.com/pub/dwalks">Densities of short uniform random walks, Part II: Higher dimensions</a>, Preprint, 2015.
%p A253095 W := proc(n,nu,twok)
%p A253095     option remember;
%p A253095     local k;
%p A253095     k := twok/2 ;
%p A253095     if n = 2 and nu = 1 then
%p A253095         binomial(2*k+2,k+1)/(k+2) ;
%p A253095     else
%p A253095         add( procname(n-1,nu,2*j)*binomial(k,j)*(k+nu)!*nu!/(k-j+nu)!/(j+nu)!,j=0..k) ;
%p A253095         simplify(%,GAMMA) ;
%p A253095     end if;
%p A253095 end proc:
%p A253095 A253095 := proc(n)
%p A253095     W(4,1,n) ;
%p A253095 end proc:
%p A253095 seq(A253095(2*n),n=0..25) ; # _R. J. Mathar_, Jun 14 2015
%t A253095 W[n_, nu_, twok_] := W[n, nu, twok] = Module[{k}, k = twok/2; If[n == 2 && nu == 1, Binomial[2k+2, k+1]/(k+2), Sum[W[n-1, nu, 2j]*Binomial[k, j]*(k+nu)!*nu!/(k-j+nu)!/(j+nu)!, {j, 0, k}]]];
%t A253095 A253095[n_] := W[4, 1, n];
%t A253095 Table[A253095[2n], {n, 0, 25}] (* _Jean-François Alcover_, Apr 16 2023, after _R. J. Mathar_ *)
%K A253095 nonn
%O A253095 0,2
%A A253095 _N. J. A. Sloane_, Feb 16 2015