function y=mynorm(x，mu，sigma) y=1/sqrt(2*pi)/sigma*exp(-(x-mu).^2/2/sigma.^2);end保存成mynorm.m调用时 用 mynorm(x，mu，sigma)传入适当的参数就可以了
