function y=normaldensity(x,m,s)

y=exp(-((x-m).^2)/(2*s*s));
y=y/sqrt(2*3.14152*s*s);