from numpy import *
from Jacobi import *
a = array([[5, -3, 1], [2, 4, -1], [2, -3, 8]], float)
b = array([[5], [6], [4]], float)
x = array([[1], [1], [1]], float)
x = jacobi(a, b, x); print(x)
x = jacobi(a, b, x); print(x)
x = jacobi(a, b, x); print(x)
x = jacobi(a, b, x); print(x)
x = jacobi(a, b, x); print(x)
x = jacobi(a, b, x); print(x)
x = jacobi(a, b, x); print(x)