diff --git a/GeoLib/AnalyticalGeometry.cpp b/GeoLib/AnalyticalGeometry.cpp
index 138b14f1183119865f174699bf9ccff466f20423..eb19bfb4b165366e65ae4361bb85e29dc2b227c1 100644
--- a/GeoLib/AnalyticalGeometry.cpp
+++ b/GeoLib/AnalyticalGeometry.cpp
@@ -274,28 +274,6 @@ bool lineSegmentsIntersect(const GeoLib::Polyline* ply,
     return false;
 }
 
-void computeRotationMatrixToXZ(MathLib::Vector3 const& plane_normal, MathLib::DenseMatrix<double> & rot_mat)
-{
-    // *** some frequently used terms ***
-    // n_1^2 + n_2^2
-    const double h0(plane_normal[0] * plane_normal[0] + plane_normal[1] * plane_normal[1]);
-    // 1 / sqrt (n_1^2 + n_2^2)
-    const double h1(1.0 / sqrt(h0));
-    // 1 / sqrt (n_1^2 + n_2^2 + n_3^2)
-    const double h2(1.0 / sqrt(h0 + plane_normal[2] * plane_normal[2]));
-
-    // calc rotation matrix
-    rot_mat(0, 0) = plane_normal[1] * h1;
-    rot_mat(0, 1) = -plane_normal[0] * h1;
-    rot_mat(0, 2) = 0.0;
-    rot_mat(1, 0) = plane_normal[0] * h2;
-    rot_mat(1, 1) = plane_normal[1] * h2;
-    rot_mat(1, 2) = plane_normal[2] * h2;
-    rot_mat(2, 0) = plane_normal[0] * plane_normal[2] * h1 * h2;
-    rot_mat(2, 1) = plane_normal[1] * plane_normal[2] * h1 * h2;
-    rot_mat(2, 2) = -sqrt(h0) * h2;
-}
-
 void rotatePoints(MathLib::DenseMatrix<double> const& rot_mat, std::vector<GeoLib::Point*> &pnts)
 {
     rotatePoints(rot_mat, pnts.begin(), pnts.end());
diff --git a/GeoLib/AnalyticalGeometry.h b/GeoLib/AnalyticalGeometry.h
index 58a06eb54277170c71dc3cf553d9bedefc917b5a..8cc6c1644265a747683e34503b0050cc28463b94 100644
--- a/GeoLib/AnalyticalGeometry.h
+++ b/GeoLib/AnalyticalGeometry.h
@@ -106,14 +106,6 @@ void compute3DRotationMatrixToX(MathLib::Vector3  const& v, T_MATRIX & rot_mat);
 template <class T_MATRIX>
 void computeRotationMatrixToXY(MathLib::Vector3 const& n, T_MATRIX& rot_mat);
 
-/**
- * Method computes the rotation matrix that rotates the given vector parallel to the \f$y\f$ axis.
- * @param plane_normal the (3d) vector that is rotated parallel to the \f$y\f$ axis
- * @param rot_mat 3x3 rotation matrix
- */
-void computeRotationMatrixToXZ(MathLib::Vector3 const& plane_normal,
-                               MathLib::DenseMatrix<double> & rot_mat);
-
 /**
  * rotate points according to the rotation matrix
  * @param rot_mat 3x3 dimensional rotation matrix