Question: Sample of random matrix. Normal law.

X is a real random matrix of dimension n; its entries are i.i.d. variables.

They follow the standard normal law N(0,1).

I want to write a big sample. For that, I use the following two lines:

roll := RandomVariable(Normal(0, 1))

X := Matrix(n, proc (i, j) options operator, arrow; (Sample(roll))(1)[1] end proc)

Is there a method that does the calculations faster ??

Thanks in advance.

Please Wait...