/** * Return value of beta function. * * @author Matt Steigart ? * * @param x * @param y * @return double */ double Statistics::beta(double x, double y) { if(x <= 0.0 || y <= 0.0){ throw FDistributionException(); } return(exp(gamma_log(x) + gamma_log(y) - gamma_log(x + y))); }
long double lgamma(long double z) { double v = z; return gamma_log(&v); }
float lgamma(float z) { double v = z; return (float)gamma_log(&v); }