diff --git a/Eigen/src/Core/SpecialFunctions.h b/Eigen/src/Core/SpecialFunctions.h
index 8cf26f4d1..bd022946c 100644
--- a/Eigen/src/Core/SpecialFunctions.h
+++ b/Eigen/src/Core/SpecialFunctions.h
@@ -13,6 +13,29 @@
 namespace Eigen {
 namespace internal {
 
+//  Parts of this code are based on the Cephes Math Library.
+//
+//  Cephes Math Library Release 2.8:  June, 2000
+//  Copyright 1984, 1987, 1992, 2000 by Stephen L. Moshier
+//
+//  Permission has been kindly provided by the original author
+//  to incorporate the Cephes software into the Eigen codebase:
+//
+//    From: Stephen Moshier
+//    To: Eugene Brevdo
+//    Subject: Re: Permission to wrap several cephes functions in Eigen
+//
+//    Hello Eugene,
+//
+//    Thank you for writing.
+//
+//    If your licensing is similar to BSD, the formal way that has been
+//    handled is simply to add a statement to the effect that you are incorporating
+//    the Cephes software by permission of the author.
+//
+//    Good luck with your project,
+//    Steve
+
 namespace cephes {
 
 /* polevl (modified for Eigen)
@@ -178,11 +201,6 @@ struct digamma_impl<float> {
    *     message         condition      value returned
    * psi singularity    x integer <=0      INFINITY
    */
-  /*
-    Cephes Math Library Release 2.2:  June, 1992
-    Copyright 1984, 1987, 1992 by Stephen L. Moshier
-    Direct inquiries to 30 Frost Street, Cambridge, MA 02140
-  */
   EIGEN_DEVICE_FUNC
   static float run(float xx) {
     float p, q, nz, x, s, w, y, z;
@@ -297,10 +315,6 @@ struct digamma_impl<double> {
      * psi singularity    x integer <=0      INFINITY
      */
 
-    /*
-     * Cephes Math Library Release 2.8:  June, 2000
-     * Copyright 1984, 1987, 1992, 2000 by Stephen L. Moshier
-     */
     double p, q, nz, s, w, y, z;
     bool negative;
 
diff --git a/Eigen/src/plugins/ArrayCwiseUnaryOps.h b/Eigen/src/plugins/ArrayCwiseUnaryOps.h
index e818ac588..2ce7414a1 100644
--- a/Eigen/src/plugins/ArrayCwiseUnaryOps.h
+++ b/Eigen/src/plugins/ArrayCwiseUnaryOps.h
@@ -320,9 +320,6 @@ lgamma() const
 }
 
 /** \returns an expression of the coefficient-wise digamma (psi, derivative of lgamma).
- *
- * Example: \include Cwise_digamma.cpp
- * Output: \verbinclude Cwise_digamma.out
  *
  * \sa cos(), sin(), tan()
  */