Remove separate directory for gaussian_elimination_pivoting.py

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)
@@ -863,6 +862,7 @@
   * [Newtons Second Law Of Motion](physics/newtons_second_law_of_motion.py)
   * [Photoelectric Effect](physics/photoelectric_effect.py)
   * [Potential Energy](physics/potential_energy.py)
+  * [Rainfall Intensity](physics/rainfall_intensity.py)
   * [Reynolds Number](physics/reynolds_number.py)
   * [Rms Speed Of Molecule](physics/rms_speed_of_molecule.py)
   * [Shear Stress](physics/shear_stress.py)

linear_algebra/src/gaussian_elimination_pivoting/gaussian_elimination_pivoting.py → 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) = }")