DATA602_Discussion 3

install.packages(“pracma”)

library(pracma)

A <- matrix(c(1, 3, 2, 4), nrow = 2)
A

##      [,1] [,2]
## [1,]    1    2
## [2,]    3    4

A <- charpoly(A, info = TRUE)

## Error term: 4

A$cp

## [1]  1 -5 -2

So the characteristic polynomial is

\[\lambda^2\ - 5 \lambda\ - 2 = 0 \]