0001 function r=rand(g,varargin)
0002
0003
0004
0005
0006
0007
0008
0009 if nargin>1
0010 n=varargin{1};
0011 else
0012 n=1;
0013 end
0014
0015 if g.d==0
0016 r=repmat(g.m,1,n);
0017 else
0018 L=chol(g.S);
0019 A_temp=normrnd(0,1,g.dim,n);
0020 A=L*A_temp;
0021 r=repmat(g.m,1,n)+A;
0022 end