function f = normal(x);
global MU SIGMA

f = 1./(SIGMA .* sqrt(2*pi)) .* exp(-1/2 .* ((x-MU)./SIGMA).^2 );