popN1 <- matCory %*% popN0 popN1 popN2 <- matCory %*% popN1 popN2 popN3 <- matCory %*% popN2 popN3