Call a function and provide output in a matrix

I have a function in R which I call

RS1 = t(cbind(Data[,18], Data[,20]))
RS2 = t(cbind(Data[,19], Data[,21]))
p = t(Data[23:24])

rand_x <- function (p, x) {
  n.goods <- dim (p)[1]
  n.obs   <- dim (p)[2]
  xRC = NaN*matrix(1, n.goods, n.obs)
  for(i in 1:n.obs) {
    xRC[1,i] <- RS1[1,i] + RS1[2,i]
    xRC[2,i] <- RS2[1,i] + RS2[2,i]
  }
  result <- xRC
  return(result)
}

This function by having these two inputs generates a vector (2x50) with some random numbers. I want to call this function rand_x 1000 times and derive 1000 matrices and then bind the results in a final matrix. I have tried to create a loop to sort this problem but I am still struggling. Any help will be much appreciated.