小弟想要跑三层循环的ARL，以下是我的程式码
tic
for i=1:1000
mvnrnd([0 0],[1 0; 0 1]);
end
toc
tic
for i=1:1000
[x1 x2]=BivGamRND(1, 4, 1, 4, 1, 0);
Q=[x1 x2];
end
toc
结果为
Elapsed time is 0.056313 seconds.
Elapsed time is 50.921110 seconds.
第一种产生二元常态与第二种产生二元Gamma差了1000倍左右...
这只是第一层而已，第二层j要重复1000次，第三层k要run 100次....
想请问大大如何让我的二元gamma产生数据快一点呢?
谢谢