completeMatrix = matrix(c(1:20), nrow=10, ncol=2)
for (i in 1:3) {
  newMatrix <- matrix(c(1:20), nrow=10, ncol=2)
  completeMatrix <- abind(completeMatrix, newMatrix , along = 3)
}
newMatrix <- matrix(c(1:18), nrow=9, ncol=2)
completeMatrix <- abind(completeMatrix, newMatrix, along = 3)