function [data] = gaussian_distribution_in_2d ( gX, gY, cX, cY, width )

data = 0*(1:gX)'*(1:gY);

for iX = 0:(gX-1),
   for iY = 0:(gY-1),
	data(iX+1,iY+1) = exp(- ((iX-cX)^2 + (iY-cY)^2)/(width^2));
   end
end