function output = mvnrnd(centroid,variance,flag)
for i = 1:flag,
    output(i, :) = randn(1, length(centroid)) * chol(variance) + centroid;
end;
