Solves the linear system Ax = b using Gauss-Seidel method.

GaussSeidel.py

import numpy as np

def Gauss_Seidel(A, b, N):
    """Solves the linear system Ax = b using Gauss-Seidel method."""
    x = np.zeros(len(b))
    e_inverse = np.linalg.inv(np.tril(A))
    
    for i in range(N):
        x += np.dot(e_inverse, b - np.dot(A, x))
                
    return x