id	summary	reporter	owner	description	type	status	milestone	component	version	severity	resolution	keywords	cc
5433	beta function: when b<epsilon return tgamma(b)?	Eric Butter <egbutter@…>	John Maddock	"In the beta.hpp header for the beta function (i.e., *not* the beta distribution) around line 131, you have: 

   if((c == a) && (b < tools::epsilon<T>()))
      return boost::math::tgamma(b, pol);
   else if((c == b) && (a < tools::epsilon<T>()))
      return boost::math::tgamma(a, pol);

Don't you actually want if b<epsilon, tgamma(a)?  beta(a,b) reduces to (approximately) gamma(a)b^-a for a>>b, which seems like the case for which you are trying to account?

Cheers,  

Eric "	Bugs	closed	To Be Determined	math	Boost 1.47.0	Problem	invalid