| Hosted by CoCalc | Download
    m = matrix([[0, 1, 0],
                [0, 0, 1],
                [1, 0, 0]])
    v = vector([2, 3, 1])
    m \ v
(1, 2, 3) 
    sol = m \ v
    m * sol
(2, 3, 1)