function output = normrnd(centroid,variance,size1,size2)
output = randn(size1, size2) .* sqrt(variance) + centroid;
