% Function to evaluate the area under the normal distribution

% Passes back a single argument containing f(x) for the normal PDF function

function f = normal(x);
global MU SIGMA

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