pdf("rpd2_11.pdf")

A <- matrix(c(10,3,3,8),byrow=T,ncol=2);

eigenA <- eigen(A)

lambda <- eigenA$val

Z <- eigenA$vec

lambda1i <- lambda[1]*diag(1,nrow=2,ncol=2)

lambda2i <- lambda[2]*diag(1,nrow=2,ncol=2)

detA_lambda1 <- det(A-lambda1i)

detA_lambda2 <- det(A-lambda2i)

A

lambda

Z

detA_lambda1

detA_lambda2

dev.off()