1/pi*(integrate gamma(1+i*x) gamma(1-i*x) from x=-inf to inf)