import sympy as sp sp.init_printing()
eps = sp.Symbol("epsilon") lam = sp.Symbol("lambda")
m = sp.Matrix([[1, eps], [eps, 1]]) m
m.charpoly(lam)
Observe the occurrence of $(1-\epsilon^2)$ above.