import copy

def det3(m):
    return m[0][0]*m[1][1]*m[2][2] + m[0][1]*m[1][2]*m[2][0] + m[0][2]*m[1][0]*m[2][1] - m[0][1]*m[1][0]*m[2][2] - m[0][0]*m[1][2]*m[2][1] - m[0][2]*m[1][1]*m[2][0]

def cramer3(equation):
    base_matrix = [x[:3] for x in equation]
    W = det3(base_matrix)
    print(W)
    if W==0:
        print("Nie można rozwiązać metodą Cramera")
    else:
        x_matrix = copy.deepcopy(base_matrix)
        y_matrix = copy.deepcopy(base_matrix)
        z_matrix = copy.deepcopy(base_matrix)
        for i in range(3):
            x_matrix[i][0] = equation[i][3]
            y_matrix[i][1] = equation[i][3]
            z_matrix[i][2] = equation[i][3]
        print(f'X:{det3(x_matrix)/W}, Y:{det3(y_matrix)/W}, Z:{det3(z_matrix)/W}')

equation_matrix = [
    [2,8,-1,5],
    [3,2,5,4],
    [9,-6,4,1]
]
cramer3(equation_matrix)