double betacf(const double, const double, const double);
double betai(const double, const double, const double);
double gammln(const double);
double betacf(const double, const double, const double);
double betai(const double, const double, const double);
double gammln(const double);