double gammapdf(double value, double alpha, double beta) {
return (std::pow(beta, alpha)*std::pow(value, (alpha-1))*std::pow(M_E, (-1*beta*value)))/tgamma(alpha);}
Это реализует расчет pdf гамма-распределения, как указано в википедии.Немного покопался.Cdf немного сложнее, и я не уверен, как решить эту проблему.