diff --git a/DIRECTORY.md b/DIRECTORY.md
index 54bb8f148..11de569a2 100644
--- a/DIRECTORY.md
+++ b/DIRECTORY.md
@@ -540,8 +540,7 @@
   * [Lu Decomposition](linear_algebra/lu_decomposition.py)
   * Src
     * [Conjugate Gradient](linear_algebra/src/conjugate_gradient.py)
-    * Gaussian Elimination Pivoting
-      * [Gaussian Elimination Pivoting](linear_algebra/src/gaussian_elimination_pivoting/gaussian_elimination_pivoting.py)
+    * [Gaussian Elimination Pivoting](linear_algebra/src/gaussian_elimination_pivoting.py)
     * [Lib](linear_algebra/src/lib.py)
     * [Polynom For Points](linear_algebra/src/polynom_for_points.py)
     * [Power Iteration](linear_algebra/src/power_iteration.py)
diff --git a/linear_algebra/src/gaussian_elimination_pivoting/gaussian_elimination_pivoting.py b/linear_algebra/src/gaussian_elimination_pivoting.py
similarity index 83%
rename from linear_algebra/src/gaussian_elimination_pivoting/gaussian_elimination_pivoting.py
rename to linear_algebra/src/gaussian_elimination_pivoting.py
index 2a86350e9..ecaacce19 100644
--- a/linear_algebra/src/gaussian_elimination_pivoting/gaussian_elimination_pivoting.py
+++ b/linear_algebra/src/gaussian_elimination_pivoting.py
@@ -1,15 +1,5 @@
 import numpy as np
 
-matrix = np.array(
-    [
-        [5.0, -5.0, -3.0, 4.0, -11.0],
-        [1.0, -4.0, 6.0, -4.0, -10.0],
-        [-2.0, -5.0, 4.0, -5.0, -12.0],
-        [-3.0, -3.0, 5.0, -5.0, 8.0],
-    ],
-    dtype=float,
-)
-
 
 def solve_linear_system(matrix: np.ndarray) -> np.ndarray:
     """
@@ -87,15 +77,18 @@ def solve_linear_system(matrix: np.ndarray) -> np.ndarray:
 
 if __name__ == "__main__":
     from doctest import testmod
-    from pathlib import Path
 
     testmod()
-    file_path = Path(__file__).parent / "matrix.txt"
-    try:
-        matrix = np.loadtxt(file_path)
-    except FileNotFoundError:
-        print(f"Error: {file_path} not found.  Using default matrix instead.")
 
-    # Example usage:
-    print(f"Matrix:\n{matrix}")
-    print(f"{solve_linear_system(matrix) = }")
+    example_matrix = np.array(
+        [
+            [5.0, -5.0, -3.0, 4.0, -11.0],
+            [1.0, -4.0, 6.0, -4.0, -10.0],
+            [-2.0, -5.0, 4.0, -5.0, -12.0],
+            [-3.0, -3.0, 5.0, -5.0, 8.0],
+        ],
+        dtype=float,
+    )
+
+    print(f"Matrix:\n{example_matrix}")
+    print(f"{solve_linear_system(example_matrix) = }")
diff --git a/linear_algebra/src/gaussian_elimination_pivoting/__init__.py b/linear_algebra/src/gaussian_elimination_pivoting/__init__.py
deleted file mode 100644
index e69de29bb..000000000
diff --git a/linear_algebra/src/gaussian_elimination_pivoting/matrix.txt b/linear_algebra/src/gaussian_elimination_pivoting/matrix.txt
deleted file mode 100644
index dd895ad85..000000000
--- a/linear_algebra/src/gaussian_elimination_pivoting/matrix.txt
+++ /dev/null
@@ -1,4 +0,0 @@
-5.0 -5.0 -3.0 4.0 -11.0
-1.0 -4.0 6.0 -4.0 -10.0
--2.0 -5.0 4.0 -5.0 -12.0
--3.0 -3.0 5.0 -5.0 8.0
\ No newline at end of file