OFFSET
0,1
EXAMPLE
Equals 0.5985017924490618300179988376616175749146889655229...
MAPLE
evalf(Int(exp(-x)/(1 + erf(sqrt(x))), x=0..infinity, method=_Gquad), 110);
MATHEMATICA
RealDigits[NIntegrate[1/(E^z (1 + Erf[Sqrt[z]])), {z, 0, Infinity}, PrecisionGoal -> 20]][[1]]
(* Second program: *)
digits = 105; dz = 50;
f[z0_] := f[z0] = NIntegrate[1/(E^z (1 + Erf[Sqrt[z]])), {z, 0, z0}, WorkingPrecision -> digits+5] + 1/(2 E^z0) // RealDigits[#, 10, digits][[1]]&; f[z0 = dz]; f[z0 = z0 + dz];
While[Print["z0 = ", z0]; f[z0] != f[z0 - dz], z0 = z0 + dz];
f[z0] (* Jean-François Alcover, Jun 18 2019 *)
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Peter Luschny, Jun 08 2019
EXTENSIONS
More digits from Jean-François Alcover, Jun 18 2019
STATUS
approved