import numpy as np
a = np.array([[1, -5],
			  [7, -1]], float)
b = np.array([-4, 6], float)

(n,) = np.shape(b)
x = np.full(n, 0, float)
#xnew = np.full(n, 0, float)
#xnew = np.empty(n, float)
xnew = x.copy()
iterlimit = 100
epsilon = 0.001

#for j in range(n): #memberi nilai 0 x1, x2, xn pada iterasi ke 0
#	x[j] = 0

for iteration in range(iterlimit):
	for i in range(n):
		s = 0
		for j in range(n):
			if j != i:
				s += (-a[i, j]) * x[j]
		xnew[i] = 1 / a[i, i] * (s + b[i])

	print("Iteration ke - %d" % (iteration))
	print(x)
	print("")

	if(abs(xnew - x) <= epsilon).all():
		iteration = iteration + 1
		print("Iteration ke - %d" % iteration)
		x = np.copy(xnew)
		print(x)
		print("")
		break
	else:
		x = np.copy(xnew)

print("The solution: ")
print(x)
print("Number of iteration: %d" % (iteration))