A=matrix( c( 0, 0, 0.36, 0.5, 3.9, 5, 0, 0, 23.1, 58.7, 145, 217.9, 0.26, 0.17, 23.5, 56.1, 158.2, 231.3, 0, 0, 0.64, 0.81, 0.38, 0, 0, 0, 0, 0.03, 0.49, 0.28, 0, 0, 0, 0, 0.10, 0.72 ), nrow=6, byrow=TRUE, dimnames=list(c("CHseed", "CLseed", "small", "medium", "large", "xlarge"),c("CHseed", "CLseed", "small", "medium", "large", "xlarge"))) ev <- eigen(A) lmax <- which(Re(ev$values) == max(Re(ev$values))) lmax <- lmax[1] lambda <- Re(ev$values[lmax]) W <- ev$vectors w <- abs(Re(W[, lmax])) V <- try(Conj(solve(W)), silent=TRUE) v <- abs(Re(V[lmax, ])) s <- v %o% w e <- s * A/lambda