% X = generate_gauss(muvec,covmat,numpoints);
%
% numpoints 俺薦 父怒狼 漂隆 氦磐 單撈磐 籠欽闌 積己. 
% 阿阿狼 單撈磐綽 瘤瀝茄 啊快矯救 盒器甫 蝶弗促. 
%
% 澇仿 牢磊:
% muvec: d-瞞盔狼 乞閉 氦磐
% covmat: d*d 傍盒魂 青紡
%
% 免仿 牢磊:
% X綽 青撈 d瞞盔狼 漂隆氦磐牢 單撈磐 青紡撈促. 

function X = generate_gauss(muvec,covmat,numpoints);

d = length(muvec);
if (size(covmat) == [d,d]) == 0
   error('Incompatible dimensions!');
end

X = inv(whiten(covmat)') * randn(d,numpoints) + muvec'*ones(1,numpoints);