/**
* 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);
}
