from libsolve import *

m = Matrix([
    [1, 1, -1, 1],
    [-1, -2, 2, -2],
    [1, 2, -1, 3],
    [-1, -2, 1, -2]
])

i = Matrix([
    [1, 0, 0, 0],
    [0, 1, 0, 0],
    [0, 0, 1, 0],
    [0, 0, 0, 1]
])


for mat in [m, i]:
    mat.make_S(1, 0, 2)
    mat.make_S(2, 3, 1)
    mat.make_U(0, 1)
    mat.make_U(2, 3)
    mat.make_S(2, 1, 1)
    mat.make_S(1, 2, 1)
    mat.make_S(2, 3, 1)
    mat.make_S(1, 0, -1)
    mat.make_U(1, 2)
    mat.make_D(1, -1)
    mat.make_D(2, -1)


print(m)
print(i)