diff --git "a/\352\263\265\355\225\231\354\210\230\355\225\231 1/gauss-seidel.py" "b/\352\263\265\355\225\231\354\210\230\355\225\231 1/gauss-seidel.py" index 21a5127..2ec0d17 100644 --- "a/\352\263\265\355\225\231\354\210\230\355\225\231 1/gauss-seidel.py" +++ "b/\352\263\265\355\225\231\354\210\230\355\225\231 1/gauss-seidel.py" @@ -13,14 +13,15 @@ ] B = [ - 50, - 50, - 25, - 25, + 50., + 50., + 25., + 25., ] # Initial values -u11, u21, u12, u22 = 0, 0, 0, 0 +u11, u21, u12, u22 = 0., 0., 0., 0. +print(f'U = [{u11} {u21} {u12} {u22}]') for _ in range(20): v11 = A[0][0]*u11 + A[0][1]*u21 + A[0][2]*u12 + A[0][3]*u22 + B[0] diff --git "a/\352\263\265\355\225\231\354\210\230\355\225\231 1/jacobi.py" "b/\352\263\265\355\225\231\354\210\230\355\225\231 1/jacobi.py" index 8c94748..5c51607 100644 --- "a/\352\263\265\355\225\231\354\210\230\355\225\231 1/jacobi.py" +++ "b/\352\263\265\355\225\231\354\210\230\355\225\231 1/jacobi.py" @@ -15,15 +15,17 @@ ]) B = np.array([ - 50, - 50, - 25, - 25, + 50., + 50., + 25., + 25., ]) # Initial values -U = np.array([0, 0, 0, 0]) +U = np.array([0., 0., 0., 0.]) +print(f'U = {U}') for _ in range(20): - U = A@U + B + V = A@U + B + U = A@V + B print(f'U = {U}')