Generate Simulation

Hallo everybody,

I'm new in r and I"ll appreciate some help!

I have a matrix of nrow=30 and ncoll=54,and I would like to generate 50
simulations with tha same size of the matrix!!!That is to say that I want
to generate 50 matrices -for my 50 simulations - with the same dimensions!
I took my 1st matrix according to the formula that I want to implement:
D<-mean_m + U_i*mat_DELTA
mean_m:vector
U_I:vector
mat_Delta:matrix(54,30)

Here is the loop:
#Define Simulations(IS)
?#Define Time Step(IT)

vec_IS<-c(1:50)
vec_IT<-c(1:30)
mat_delta<-matrix(nrow=54, ncol=30)
mat_DELTA<-matrix(nrow=54, ncol=30)
mat_Yr_m<-matrix(nrow=54, ncol=30)
D<-matrix(nrow=30, ncol=54)
DELTA<-(matrix(nrow=30,ncol=54))*50

for (i in 1:length(vec_IS)){
? ? ? ?for (j in 1:length(vec_IT)){
? ? ? ? ? ? ? ?Yr_m<-rnorm(54,m=0,sd=1)
? ? ? ? ? ? ? ?mat_Yr_m[,j]<-Yr_m
? ? ? ? ? ? ? ?delta_i<-lower_m%*%mat_Yr_m[,j]
? ? ? ? ? ? ? ?mat_delta[,j]<-as.vector(delta_i)
? ? ? ?}
? ? ? ?DELTA_1<-mat_delta[,1]
? ? ? ?DELTA_2<-mat_delta[,2]-a_1*mat_delta[,1]
? ? ? ? ? ? ? ?for (t in 3:length(vec_IT)){
? ? ? ? ? ? ? ? ? ? ? ?mat_DELTA[,t]<-mat_delta[,t] - a_1*mat_delta[,t-1]
- a_2*mat_delta[,t-2]
? ? ? ? ? ? ? ?}
? ? ? ? ? ? ? ?mat_DELTA[,1]<-DELTA_1
? ? ? ? ? ? ? ?mat_DELTA[,2]<-DELTA_2
? ? ? ? ? ? ? ? ? ? ? ?for (k in 1:length(nmesh)){
? ? ? ? ? ? ? ? ? ? ? ?mean_m<-as.numeric(vec_mean_col)
? ? ? ? ? ? ? ? ? ? ? ?DELTA<-U_i*mat_DELTA
? ? ? ? ? ? ? ? ? ? ? ?D<-mean_m + t(DELTA)
? ? ? ? ? ? ? ?}

? ? ? ?}

I want to implement this formula for 50 simulations!

any idea?

Thanks a lot!

? ? ? ?[[alternative HTML version deleted]]