Use "log1p(x)" instead of "log(1 + x)".
diff --git a/commons-statistics-distribution/src/main/java/org/apache/commons/statistics/distribution/TDistribution.java b/commons-statistics-distribution/src/main/java/org/apache/commons/statistics/distribution/TDistribution.java
index a579e15..e5baad3 100644
--- a/commons-statistics-distribution/src/main/java/org/apache/commons/statistics/distribution/TDistribution.java
+++ b/commons-statistics-distribution/src/main/java/org/apache/commons/statistics/distribution/TDistribution.java
@@ -77,7 +77,7 @@
     @Override
     public double logDensity(double x) {
         final double nPlus1Over2 = dofOver2 + 0.5;
-        return factor - nPlus1Over2 * Math.log(1 + x * x / degreesOfFreedom);
+        return factor - nPlus1Over2 * Math.log1p(x * x / degreesOfFreedom);
     }
 
     /** {@inheritDoc} */