function mr = mom(x)
n = length(x); 
mu = mean(x); 
mr = (1/n)*sum( (x-mu).^2 );