eulergamma x (1 + (eulergamma-1)x) +x^2(pi^2/6+0.251x+2.33x^2), where x=(2-sqrt pi)/(2.+sqrt pi)