math_22.diff

diff --git a/src/main/java/org/apache/commons/math3/distribution/FDistribution.java b/src/main/java/org/apache/commons/math3/distribution/FDistribution.java
index 8b0993c4d..9c57c8224 100644
--- a/src/main/java/org/apache/commons/math3/distribution/FDistribution.java
+++ b/src/main/java/org/apache/commons/math3/distribution/FDistribution.java
@@ -126,6 +126,33 @@ public class FDistribution extends AbstractRealDistribution {
      * @since 2.1
      */
     public double density(double x) {
+        if (Boolean.parseBoolean(System.getProperty("defects4j.instrumentation.enabled"))) { // defects4j.instrumentation
+            double returnValue = density_original(x); // defects4j.instrumentation
+            double upperBound = getSupportUpperBound(); // defects4j.instrumentation
+            double lowerBound = getSupportLowerBound(); // defects4j.instrumentation
+            if (x == upperBound) { // defects4j.instrumentation
+                boolean upperBoundRule = !isSupportUpperBoundInclusive() || !Double.isInfinite(returnValue) && !Double.isNaN(returnValue); // defects4j.instrumentation
+                if (!upperBoundRule) { // defects4j.instrumentation
+                    throw new RuntimeException("[Defects4J_BugReport_Violation]"); // defects4j.instrumentation
+                } else { // defects4j.instrumentation
+                    return returnValue; // defects4j.instrumentation
+                } // defects4j.instrumentation
+            } else if (x == lowerBound) { // defects4j.instrumentation
+                boolean lowerBoundRule = !isSupportLowerBoundInclusive() || !Double.isInfinite(returnValue) && !Double.isNaN(returnValue); // defects4j.instrumentation
+                if (!lowerBoundRule) { // defects4j.instrumentation
+                    throw new RuntimeException("[Defects4J_BugReport_Violation]"); // defects4j.instrumentation
+                } else { // defects4j.instrumentation
+                    return returnValue; // defects4j.instrumentation
+                } // defects4j.instrumentation
+            } else { // defects4j.instrumentation
+                return returnValue; // defects4j.instrumentation
+            } // defects4j.instrumentation
+        } else { // defects4j.instrumentation
+            return density_original(x); // defects4j.instrumentation
+        } // defects4j.instrumentation
+    } // defects4j.instrumentation
+
+    public double density_original(double x) { // defects4j.instrumentation
         final double nhalf = numeratorDegreesOfFreedom / 2;
         final double mhalf = denominatorDegreesOfFreedom / 2;
         final double logx = FastMath.log(x);
diff --git a/src/main/java/org/apache/commons/math3/distribution/UniformRealDistribution.java b/src/main/java/org/apache/commons/math3/distribution/UniformRealDistribution.java
index 5d32f6ebf..ae6afafd8 100644
--- a/src/main/java/org/apache/commons/math3/distribution/UniformRealDistribution.java
+++ b/src/main/java/org/apache/commons/math3/distribution/UniformRealDistribution.java
@@ -106,6 +106,33 @@ public class UniformRealDistribution extends AbstractRealDistribution {

     /** {@inheritDoc} */
     public double density(double x) {
+        if (Boolean.parseBoolean(System.getProperty("defects4j.instrumentation.enabled"))) { // defects4j.instrumentation
+            double returnValue = density_original(x); // defects4j.instrumentation
+            double upperBound = getSupportUpperBound(); // defects4j.instrumentation
+            double lowerBound = getSupportLowerBound(); // defects4j.instrumentation
+            if (x == upperBound) { // defects4j.instrumentation
+                boolean upperBoundRule = !isSupportUpperBoundInclusive() || !Double.isInfinite(returnValue) && !Double.isNaN(returnValue); // defects4j.instrumentation
+                if (!upperBoundRule) { // defects4j.instrumentation
+                    throw new RuntimeException("[Defects4J_BugReport_Violation]"); // defects4j.instrumentation
+                } else { // defects4j.instrumentation
+                    return returnValue; // defects4j.instrumentation
+                } // defects4j.instrumentation
+            } else if (x == lowerBound) { // defects4j.instrumentation
+                boolean lowerBoundRule = !isSupportLowerBoundInclusive() || !Double.isInfinite(returnValue) && !Double.isNaN(returnValue); // defects4j.instrumentation
+                if (!lowerBoundRule) { // defects4j.instrumentation
+                    throw new RuntimeException("[Defects4J_BugReport_Violation]"); // defects4j.instrumentation
+                } else { // defects4j.instrumentation
+                    return returnValue; // defects4j.instrumentation
+                } // defects4j.instrumentation
+            } else { // defects4j.instrumentation
+                return returnValue; // defects4j.instrumentation
+            } // defects4j.instrumentation
+        } else { // defects4j.instrumentation
+            return density_original(x); // defects4j.instrumentation
+        } // defects4j.instrumentation
+    } // defects4j.instrumentation
+
+    public double density_original(double x) { // defects4j.instrumentation
         if (x < lower || x > upper) {
             return 0.0;
         }