Cumulative distribution function for the Inverse-gamma distribution gnuplot source under GPL: dinvgamma(x,a,b) = exp(a*log(b) - lgamma(a) - (a+1)*log(x)...
