I'm fairly new with R and I can't get my head around a for loop. I've developed a big markov model based on a loop (now I know it's not the most efficient way but no time to change it now). the final outcome of the model is a pair of matrices for each value of one parameter j. A sample of what I'm trying to do is 


 for (j in 1:10){
 x<-3:13
 y<-rep(x[j],10)
 a<-matrix(y,5,2)
 y<-rep(2,10)
 b<-matrix(y,5,2)
 outcome<-array(0,dim=c(5,2,10))
 outcome[,,j]<-a*b
 
 }
 outcome

that is: saving in the third dimension of the array the output of every loop. But as I see it saves only the last loop in the last dimension of the array. The model is much more complex so that I can't really change the basic structure but only find a way of saving its outputs. Any thoughts?
hope it's clear and any help would be really appreciated!