OFFSET
0,2
COMMENTS
Denominators of expansion of Integral_{t=0..x} exp(-(t^2)/2) dt = sqrt(Pi/2)*erf(x/sqrt(2)) in powers x^(2*n+1), n >= 0. Numerators are (-1)^n. - Wolfdieter Lang, Jun 29 2007
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..100
C. Nicholson, The probability integral for two variables, Biometrika 33 (1943), 59-72.
Eric Weisstein's World of Mathematics, Normal Distribution Function.
FORMULA
E.g.f.: (1+2x)/(1-2x)^2.
G.f.: G(0)/(2*x) - 1/x, where G(k)= 1 - 2*x+ 1/(1 - 2*x*(k+1)/(2*x*(k+1) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
From Amiram Eldar, Jul 31 2020: (Start)
Sum_{n>=0} 1/a(n) = sqrt(Pi/2) * erfi(1/sqrt(2)).
Sum_{n>=0} (-1)^n/a(n) = sqrt(Pi/2) * erf(1/sqrt(2)). (End)
V(h, q) = -h/(q*sqrt(2*Pi)) + Sum_{k>=0} (-1)^k*h*q^(2*k-1)*(q^2+(2*k+1))/(a(k)*sqrt(2*Pi)) = (h/2)*erf(q/sqrt(2)) + h*(exp(-q^2/2) - 1)/(q*sqrt(2*Pi)), where V is Nicholson's V-function. V(h, q) = Integral_{x=0..h} Integral_{y=0..q*x/h} phi(x)*phi(y) dydx, where phi(x) is the standard normal density exp(-x^2/2)/sqrt(2*Pi). - Thomas Scheuerle, Jan 21 2025
MATHEMATICA
a[n_]:=2^n*n!*(2*n+1); Array[a, 18, 0] (* Stefano Spezia, Jan 03 2025 *)
PROG
(Magma) [2^n*Factorial(n)*(2*n+1): n in [0..50]]; // Vincenzo Librandi, Apr 25 2011
(Haskell)
a014481 n = a009445 n `div` a001147 n -- Reinhard Zumkeller, Dec 03 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved