pdf("rpd2_12.pdf") A <- matrix(c(10,3,3,8),byrow=T,ncol=2); eigenA <- eigen(A) lambda <- eigenA$val Z <- eigenA$vec A1 <- lambda[1]*Z[,1]%*%t(Z[,1]) A2 <- lambda[2]*Z[,2]%*%t(Z[,2]) A1A2 <- A1+A2 A A1 A2 A1A2 dev.off()