diff --git a/GeoLib/Triangle.cpp b/GeoLib/Triangle.cpp
index c7b6c2574c9b8b30273ca00128f0561d18775027..6ddb2b2fad92b56278fea81dc21f0485c97f3ea0 100644
--- a/GeoLib/Triangle.cpp
+++ b/GeoLib/Triangle.cpp
@@ -33,26 +33,4 @@ bool Triangle::containsPoint(MathLib::Point3d const& q, double eps) const
     GeoLib::Point const& c(*(_pnts[_pnt_ids[2]]));
     return MathLib::isPointInTriangle(q, a, b, c, eps);
 }
-
-void getPlaneCoefficients(Triangle const& tri, double c[3])
-{
-    GeoLib::Point const& p0 (*(tri.getPoint(0)));
-    GeoLib::Point const& p1 (*(tri.getPoint(1)));
-    GeoLib::Point const& p2 (*(tri.getPoint(2)));
-    Eigen::Matrix3d mat;
-    mat(0,0) = p0[0];
-    mat(0,1) = p0[1];
-    mat(0,2) = 1.0;
-    mat(1,0) = p1[0];
-    mat(1,1) = p1[1];
-    mat(1,2) = 1.0;
-    mat(2,0) = p2[0];
-    mat(2,1) = p2[1];
-    mat(2,2) = 1.0;
-    Eigen::Vector3d y;
-    y << p0[2], p1[2], p2[2];
-
-    Eigen::Map<Eigen::Vector3d>(c,3) = mat.partialPivLu().solve(y);
-}
-
-} // end namespace GeoLib
+}  // namespace GeoLib
diff --git a/GeoLib/Triangle.h b/GeoLib/Triangle.h
index 23f6e5dbfa0db39c6973203f2af671d66171da1f..e9388cd12bff3bb3d9eda0306fb0564d2c125688 100644
--- a/GeoLib/Triangle.h
+++ b/GeoLib/Triangle.h
@@ -68,7 +68,4 @@ private:
     /// position of pointers to the geometric points
     std::array<std::size_t, 3> _pnt_ids;
 };
-
-void getPlaneCoefficients(Triangle const& tri, double c[3]);
-
-} // end namespace GeoLib
+}  // namespace GeoLib