from matplotlib import pyplot as plt
def solver_visualizer(A:np.ndarray,b:np.ndarray,x:np.ndarray,steps:int=100):
np.linspace(ax_min - 0.2 * rx,ax_max + 0.2 * rx,steps),
np.linspace(ay_min - 0.2 * ry,ay_max + 0.2 * ry,steps),
A @ np.array([X.reshape(-1),Y.reshape(-1)]) - b,
plt.figure(figsize=(10,6))
plt.plot(x[:,0],x[:,1],'-o')
A = np.array([[1,0.2],[-0.3,2]])
b = np.array([0.6,0.8]).reshape(2,1)
init_x = b / np.diag(A).reshape(2,1)
solver_visualizer(A,b,gauss_seidel_iteration(A,b,init_x,10))