Untitled

def jacobi(a,b,x,n,iterations):
    x_new=x.copy()
    
    for it in range(iterations):
        for i in range(n):
            sum1=0
            for j in range(n):
                if j!=i:
                    sum1+=a[i][j]*x[j]
            x_new[i]=(b[i]-sum1)/a[i][i]
        x=x_new.copy()
    return x
n=3
a=[[4,-1,1],
   [-1,5,-1],
   [2,-1,6]]
b=[7,-8,6]
x=[0,0,0]

iteration=10
solution=jacobi(a,b,x,n,iteration)
print("solution")
for i in range(n):
    print(f"x{i+1}={solution[i]}")