OFFSET
1,2
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..20000
J.-P. Allouche and T. Baruchel, Variations on an error sum function for the convergents of some powers of e, arXiv preprint arXiv:1408.2206 [math.NT], 2014.
FORMULA
c = sqrt(e*Pi/2)*erf(1/sqrt(2)), or 2^(-1/2)*exp(1/2)*sqrt(Pi)*(1 - erfc(1/sqrt(2))). - Michael Kleber, Mar 21 2001
From Peter Bala, Feb 09 2024: (Start)
Generalized continued fraction expansion:
c = 1/(1 - 1/(4 - 3/(6 - 5/(8 - 7/(10 - 9/(12 - ... )))))). See A286286.
c/(1 + c) = Sum_{n >= 0} (2*n-1)!!/(A112293(n) * A112293(n+1)) = 1/(1*2) + 1/(2*7) + 3/(7*36) + 15/(36*253) + 105/(253*2278) + ... = 0.5851803411..., a rapidly converging series. (End)
Equals Sum_{n >= 0} ((n - 1)*(n + 1)!*2^(n + 1))/(2*n)!. - Antonio Graciá Llorente, Feb 13 2024
EXAMPLE
1.410686134642447997690824711419115041323478...
MATHEMATICA
RealDigits[ Sqrt[E*Pi/2] * Erf[1/Sqrt[2]], 10, 107] // First
(* or *) 1/Fold[Function[2*#2-1+(-1)^#2*#2/#1], 1, Reverse[Range[100]]] // N[#, 107]& // RealDigits // First (* Jean-François Alcover, Mar 07 2013, updated Sep 19 2014 *)
PROG
(PARI) { default(realprecision, 20080); x=2^(-1/2)*exp(1/2)*sqrt(Pi)*(1 - erfc(1/sqrt(2))); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b060196.txt", n, " ", d)); } \\ Harry J. Smith, Jul 02 2009
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Evan Michael Adams (evan(AT)tampabay.rr.com), Simon Plouffe, Mar 21 2001
EXTENSIONS
More terms from Vladeta Jovovic, Mar 27 2001
STATUS
approved