litery=CharacterRange["a","z"]

n=3

A=Table[litery[[ii]]<>litery[[jj]],{ii,1,n},{jj,1,n}]

wielomiancharakterystyczny = CharacteristicPolynomial[A,λ]

wyznacznik=Det[A]

wielomiancharakterystyczny-wyznacznik /. {λ->A,Power->MatrixPower}

% + IdentityMatrix[n] wyznacznik

Simplify[%]