using SpecialFunctions gamma(1/2) - sqrt(pi)