diff --git a/fermi b/fermi
index 1e308de..c43d74c 100755
Binary files a/fermi and b/fermi differ
diff --git a/fermi.go b/fermi.go
index 43b4874..a391888 100644
--- a/fermi.go
+++ b/fermi.go
@@ -244,10 +244,6 @@ func multiplyLogDists(l1 Lognormal, l2 Lognormal) Lognormal {
 
 }
 
-func multiplyBetaDists(beta1 Beta, beta2 Beta) Beta {
-	return Beta{a: beta1.a + beta2.a, b: beta1.b + beta2.b}
-}
-
 func multiplyLogDistAndScalar(l Lognormal, s Scalar) (Dist, error) {
 	if s == 0.0 {
 		return Scalar(0.0), nil
@@ -285,11 +281,6 @@ func multiplyDists(old_dist Dist, new_dist Dist) (Dist, error) {
 				return Scalar(float64(o) * float64(n)), nil
 			}
 		}
-	case Beta:
-		switch n := new_dist.(type) {
-		case Beta:
-			return multiplyBetaDists(o, n), nil
-		}
 	}
 	return operateDistsAsSamples(old_dist, new_dist, "*")
 }
diff --git a/simple/f_simple.go b/simple/f_simple.go
index 85838db..c712a71 100644
--- a/simple/f_simple.go
+++ b/simple/f_simple.go
@@ -173,10 +173,6 @@ func multiplyLogDists(l1 Lognormal, l2 Lognormal) Lognormal {
 
 }
 
-func multiplyBetaDists(beta1 Beta, beta2 Beta) Beta {
-	return Beta{a: beta1.a + beta2.a, b: beta1.b + beta2.b}
-}
-
 func multiplyDists(old_dist Dist, new_dist Dist) (Dist, error) {
 
 	switch o := old_dist.(type) {
@@ -205,16 +201,8 @@ func multiplyDists(old_dist Dist, new_dist Dist) (Dist, error) {
 				return operateDistsAsSamples(old_dist, new_dist, "*")
 			}
 		}
-	case Beta:
-		switch n := new_dist.(type) {
-		case Beta:
-			return multiplyBetaDists(o, n), nil
-		default:
-			return operateDistsAsSamples(old_dist, new_dist, "*")
-		}
-	default:
-		return operateDistsAsSamples(old_dist, new_dist, "*")
 	}
+	return operateDistsAsSamples(old_dist, new_dist, "*")
 }
 
 func divideDists(old_dist Dist, new_dist Dist) (Dist, error) {