function f= gh_log_fun( para,num,x);
mean=para(1,1);std=para(1,2);
y=1/sqrt(2*pi)/std*exp(-(x-mean).^2/2/std^2);
if num==1
f=log(y);
else
f=-sum(log(y));
end
注x是观测数据（最大似然法估计正态分布的均值和方差）
在命令窗口输入[para]=my_mle('gh_log_fun',[0.5;0;0;1;1],x)

para =

   6.9551e+03 结果只有一个参数估计量 我的正态分布明明有两个参数需要估计的