# logGamma(x)=log(gamma(x)); logGamma(1):=0; logGamma(oo):=oo;