Através do uso da matriz A, com os elementos [1,3],[4,5]; depois que a matriz é inserida, o comando expand(charpoly(A, x)) é usado, onde A é o nome da matriz, x a variável polinomial, obtendo como resultado x^2-6*x-7