Matrix has det, inv, rref. Print them for a 2×2.

Goal

Invert [[1,2],[3,4]].

import sympy as sp
M = sp.Matrix([[1, 2], [3, 4]])
print(M.det())
print(M.inv())

import sympy as sp
M = sp.Matrix([[1, 2], [3, 4]])
print(M.rref())

import sympy as sp
print(sp.Matrix([[1, 0], [0, 1]]).eigenvals())

import sympy as sp
print(sp.Matrix([1, 2, 3]).T * sp.Matrix([4, 5, 6]))