pdf("rpd2_3.pdf") B <- matrix(c(1,3,2,4,5,6,8,7,9),byrow=T,ncol=3) detB <- det(B) invB <- solve(B) B detB invB BinvB <- B%*% invB BinvB dev.off()