Hi,

I would like to create M paths of 2 correlated brownian motion incrementals
where each path is of length N. I use mvrnorm to create all the increments,
i.e.

h      <-  1.0;
COV <-  matrix(c(1,0,0,1),nrow=2);
dW    <-  h * t(mvrnorm(n=N*M,mu=c(0,0),Sigma=COV));

The next step is that I'd like to wrap dW (2D matrix of size 2x(NM) into a
3D array where each slice is 2xN matrix and the entire array consists of M
slices. Just wondering if there is any way to do that? Or is there a better
way than the way I am describing? Thank you.

Roebert

	[[alternative HTML version deleted]]

