% Function to evaluate the first moment of the normal exponential


function f = normals(x);
global MU SIGMA

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