f(d) = (1/(a (CapitalGamma(x+1))))((x-m)/a)^(CapitalGamma(x+1)) e^-((x-m)/a) dx