r <- 4; d <- 2; p <- 3
Lamd <- matrix(sample(1:50-25, d*p), nrow=d)
ieg<- eigen(diag(p)+t(Lamd)%*%Lamd)
V <- ieg$vectors
Delta <-Lamd %*% V %*% diag(1/sqrt(ieg$values)) %*% t(V)
MomCum <- MomCumCFUSN(r,d,p,Delta)
Run the code above in your browser using DataLab