p=c(.4,.1,.2,.3) dp=diag(p) sig=dp-p%*%t(p) eigen(sig)