I am using the following code. It's returning the block matrix (Z) raised to negative one (think about inputting 22/7 in a Casio fx-991ES PLUS).
import sympy as sp
from IPython.display import display
X = sp.Matrix([[1, 1, 1], [2, 2, 2], [3, 3, 3]])
i = sp.Matrix([[1], [1], [1]])
Z =...