报童问题(3)

D=var(x,1) s=std(x,1) D1=var(x) s1=std(x)

[muhat,sigmahat,muci,sigmaci]= normfit(x)

附录4：求解问题二中报纸买进量及最大收益的程序 Myfun1.m

function f=myfun1(x,n);

f=0.2.*n.*0.2044.*exp(-1.*((x-190.1)./54.98).^2);

myfun2.m

function f=myfun2(x,n);

f=(0.2.*x-0.1.*(n-x)).*0.2044.*exp(-1.*((x-190.1)./54.98).^2);

主程序 i=1;

for n=100:1:300

c(i)=quad(@(x)myfun2(x,n),0,n)+quad(@(x)myfun1(x,n),n,3000); i=i+1; end

[d,i]=max(c)

11