TR5 gauss

def gauss():
    matrix = [[1, 2, 8],
              [7, 1, 6],
              [1, 7, 2]]

    ans = [-5,
           -11,
           11]
    matrix1, matrix2, matrix3 = [[0, 0, 0], [0, 0, 0], [0, 0, 0]], [[0, 0, 0], [0, 0, 0], [0, 0, 0]], [[0, 0, 0], [0, 0, 0], [0, 0, 0]]
    ans1, ans2, ans3 = [0, 0, 0], [0, 0, 0], [0, 0, 0]
    matrix1[0][0] = matrix[0][0] / matrix[0][0]
    matrix1[0][1] = matrix[0][1] / matrix[0][0]
    matrix1[0][2] = matrix[0][2] / matrix[0][0]
    ans1[0] = ans[0] / matrix[0][0]

    matrix1[1][0] = matrix[1][0] - matrix[1][0] * (matrix1[0][0] / matrix[0][0])
    matrix1[1][1] = matrix[1][1] - matrix[1][0] * (matrix1[0][1] / matrix[0][0])
    matrix1[1][2] = matrix[1][2] - matrix[1][0] * (matrix1[0][2] / matrix[0][0])

    ans1[1] = ans[1] - matrix[1][0] * (ans1[0] / matrix[0][0])

    matrix1[2][0] = matrix[2][0] - matrix[2][0] * (matrix1[0][0] / matrix[0][0])
    matrix1[2][1] = matrix[2][1] - matrix[2][0] * (matrix1[0][1] / matrix[0][0])
    matrix1[2][2] = matrix[2][2] - matrix[2][0] * (matrix1[0][2] / matrix[0][0])

    ans1[2] = ans[2] - matrix[2][0] * (ans1[0] / matrix[0][0])


    matrix2
    matrix2[1][1] = matrix1[1][1] / matrix1[1][1]
    matrix2[1][2] = matrix1[1][2] / matrix1[1][1]
    ans2[1] = ans1[1] / matrix1[1][1]

    matrix2[2][1] = matrix1[2][1] - matrix1[2][1] * (matrix1[1][1] / matrix1[1][1])
    matrix2[2][2] = matrix1[2][2] - matrix1[2][1] * (matrix1[1][2] / matrix1[1][1])
    ans2[2] = ans1[2] - matrix1[2][1] * (ans1[1] / matrix1[1][1])

    matrix3[2][2] = matrix2[2][2] / matrix2[2][2]
    ans3[2] = ans2[2] / matrix2[2][2]

    x3 = ans3[2]
    x2 = ans2[1] - matrix2[1][2] * x3
    x1 = ans1[0] - matrix1[0][1] * x2 - matrix1[0][2] * x3
    print(x1, x2, x3)