Index: tol/bmath/bstat/stat.cpp
===================================================================
--- tol/bmath/bstat/stat.cpp	(revision 4aebed6aea2aa2a383725a38127dcfce40d5185f)
+++ tol/bmath/bstat/stat.cpp	(revision dbf50f65d993aff72dc60cc4b22f2533930dc264)
@@ -1098,8 +1098,8 @@
   int q2 = max(0,2*m-n+1)/2;
   int q3 = max(0,3*m-n+1)/2;
-  int q4 = max(0,4*m-n+1)/2;
-  double EQBP2 = (1/(n2-1)*(n+3)*(n+5)) * 
+  int q4 = max(0,4*m-n+2)/2;
+  double EQBP2 = (1/((n2-1)*(n+3)*(n+5))) * 
   (
-    (m/12*n2)*(
+    (m/(12*n2))*(
       24*n*(-2-2*n-n2+3*n3-2*n4+n5)
       +(35+114*n2+60*n3-261*n4+84*n5+12*n6)*m
@@ -1112,5 +1112,5 @@
     +4*n*q4*(q4+1)*(4*q4-12*m+3*n-4)
     +4*q3*(q3+1)*(2*q3*(q3+1)-(3*m-n+1)*(3*m-n+2)-1)
-    +(q2/6*n)*(
+    +(q2/(6*n))*(
       -(2*m-n+2)*(
         (24-84*n-150*n2+59*n3+16*n4)