function output = gaussDistribution(r,z)

cr = 0;
cz = 0.15;

dr = abs(r-cr);
dz = abs(z-cz);
mu = 0;
 tao = 10000;
%tao = 1000000;
dis = sqrt(dr.^2 + dz.^2);
output = 0.*dis;
for index_i = 1:length(dis)

    %if dr(index) > 0.004 || dz(index)> 0.004
     %   output(index_i) = 0;
    %else
    output(index_i) =tao*1e9./sqrt(2*pi).*exp(-1*tao.*(dis(index_i) - mu).^2/2);
    %end 
end