Index: tol/bmath/matfun.cpp
===================================================================
--- tol/bmath/matfun.cpp	(revision bb4284bf2d937f0bac8162db4766f0bfef884dee)
+++ tol/bmath/matfun.cpp	(revision b8f55afafaa61e3fa51f432e800a65f731014274)
@@ -30,4 +30,6 @@
 #include <gsl/gsl_math.h>
 #include <gsl/gsl_sf_hyperg.h>
+
+#include <boost/math/common_factor_rt.hpp>
 
 #define LIMITFactorial	100
@@ -996,3 +998,22 @@
 };
 
-
+//--------------------------------------------------------------------
+BInt GCD(const BInt& a, const BInt& b)
+      
+/*! Returns the Greatest Common Divisor of two integers
+ */
+//--------------------------------------------------------------------
+{
+    return(boost::math::gcd(a, b));
+};
+
+//--------------------------------------------------------------------
+BInt LCM(const BInt& a, const BInt& b)
+      
+/*! Returns the Least Common Multiple of two integers
+ */
+//--------------------------------------------------------------------
+{
+    return(boost::math::lcm(a, b));
+};
+
Index: tol/bmath/mathobjects/tol_bpolyn.h
===================================================================
--- tol/bmath/mathobjects/tol_bpolyn.h	(revision bb4284bf2d937f0bac8162db4766f0bfef884dee)
+++ tol/bmath/mathobjects/tol_bpolyn.h	(revision b8f55afafaa61e3fa51f432e800a65f731014274)
@@ -24,5 +24,5 @@
 
 #include <tol/tol_bmonome.h>
-#include <boost/math/common_factor_rt.hpp>
+#include <tol/tol_bmatfun.h>
 #include <limits>
 
@@ -370,5 +370,5 @@
   {
     deg = (*this)(i).Degree();
-    if(deg) { gcd = boost::math::gcd(gcd,deg); }
+    if(deg) { gcd = GCD(gcd,deg); }
   }
   return(gcd);
Index: tol/bmath/tol_bmatfun.h
===================================================================
--- tol/bmath/tol_bmatfun.h	(revision bb4284bf2d937f0bac8162db4766f0bfef884dee)
+++ tol/bmath/tol_bmatfun.h	(revision b8f55afafaa61e3fa51f432e800a65f731014274)
@@ -69,4 +69,7 @@
 BDat InvFactorial(const BDat& x);
 
+BInt GCD(const BInt& m, const BInt& n);
+BInt LCM(const BInt& m, const BInt& n);
+
 #endif	// TOL_BMATFUN_H
 
