NumberForm[NIntegrate[e^(-x^2/2) ln(e^x+1),{x,-infty,infty}],20]