I'm trying to use the DOZZ formula in Mathematica, but I'm not sure what to do about these special functions. I had hoped to use some built-in functions, or at least work with a nice series expansion, but things like the multiple gamma function seem tough to include. Have previous works on Liouville field theory used the full formula, or do they just rely on asymptotics of the upsilon function/integral representation?