integral 1 - (γ(eulergamma,x))^2/Gamma(eulergamma)^2 from 0 to infty