diff --git a/unsupported/Eigen/NonLinear b/unsupported/Eigen/NonLinear
index e98a907a5..334b6d65d 100644
--- a/unsupported/Eigen/NonLinear
+++ b/unsupported/Eigen/NonLinear
@@ -45,7 +45,6 @@ namespace Eigen {
 #include "src/NonLinear/r1mpyq.h"
 #include "src/NonLinear/rwupdt.h"
 #include "src/NonLinear/qrfac.h"
-#include "src/NonLinear/fdjac2.h"
 #include "src/NonLinear/fdjac1.h"
 #include "src/NonLinear/qform.h"
 #include "src/NonLinear/lmpar.h"
diff --git a/unsupported/Eigen/src/NonLinear/fdjac2.h b/unsupported/Eigen/src/NonLinear/fdjac2.h
deleted file mode 100644
index cb360d84c..000000000
--- a/unsupported/Eigen/src/NonLinear/fdjac2.h
+++ /dev/null
@@ -1,35 +0,0 @@
-
-template<typename FunctorType, typename Scalar>
-int ei_fdjac2(
-        const FunctorType &Functor,
-        Matrix< Scalar, Dynamic, 1 >  &x,
-        Matrix< Scalar, Dynamic, 1 >  &fvec,
-        Matrix< Scalar, Dynamic, Dynamic > &fjac,
-        Scalar epsfcn)
-{
-    /* Local variables */
-    Scalar h, temp;
-    int iflag;
-
-    /* Function Body */
-    const Scalar epsmch = epsilon<Scalar>();
-    const int n = x.size();
-    const Scalar eps = ei_sqrt((std::max(epsfcn,epsmch)));
-    Matrix< Scalar, Dynamic, 1 >  wa(fvec.size());
-
-    for (int j = 0; j < n; ++j) {
-        temp = x[j];
-        h = eps * ei_abs(temp);
-        if (h == 0.) {
-            h = eps;
-        }
-        x[j] = temp + h;
-        iflag = Functor(x, wa);
-        if (iflag < 0)
-            return iflag;
-        x[j] = temp;
-        fjac.col(j) = (wa-fvec)/h;
-    }
-    return iflag;
-}
-