a = [1,9,6]
b = [0,1,0]
c = [7,1,9]
svk1 = 8
svk2 = 9
svk3 = 0

L1=-b[0]/a[0] # Значение коэффициента, на который нужно будет умножить каждый элемент первой строки
L2=-c[0]/a[0]

'''
В следующем цикле находим значение элементов 2 и 3 строки после умножения на коэффициент L1,L2 на элементы первой строки, затем вычитается каждый элемент 2 и 3 строк
'''

for j in range(3):
    b[j]=a[j]*L1+b[j]
    c[j]=a[j]*L2+c[j]

svk2 = svk1*L1+svk2 #Нахождение свободного коэффициента
svk3 = svk1*L2+svk3

print('Матрица после первого шага')
print(a)
print(b)
print(c)

print('Коэффициенты после первого шага')
print(svk1)
print(svk2)
print(svk3)


L1=-a[1]/b[1]
L2=-c[1]/b[1]

for j in range(3):
    a[j]=b[j]*L1+a[j]
    c[j]=b[j]*L2+c[j]

svk1 = svk2*L1+svk1
svk3 = svk2*L2+svk3


print('Матрица после второго шага')
print(a)
print(b)
print(c)

print('Коэффициенты после второго шага')
print(svk1)
print(svk2)
print(svk3)


L1=-a[2]/c[2]
L2=-b[2]/c[2]

for j in range(3):
    a[j]=c[j]*L1+a[j]
    b[j]=c[j]*L2+b[j]

svk1 = svk3*L1+svk1
svk2 = svk3*L2+svk2


print('Матрица после третьего шага')
print(a)
print(b)
print(c)

print('Коэффициенты после третьего шага')
print(svk1)
print(svk2)
print(svk3)

print(a[1])

x1 = round((svk1/a[0]),4)
x2 = round((svk2/b[1]),4)
x3 = round((svk3/c[2]),4)

print('Результаты')
print('X1= ', x1)
print('X2= ', x2)
print('X3= ', x3)