用MATLAB实现报童模型，求大神帮我看下程序，万分感谢！！

我在做“需求信息对零售商订货策略影响”的论文，用报童模型，用MATLAB编了出来，但是验算出来的和答案不同。求大神帮我看看，很急，万分谢谢。
function Untitled2
clear;clc
w1 = 35;
p = 60;
w2 = 45;
r = 10;
miu = 2000;
sigma = 600;
for w1 = 35:45  
    v1 = norminv((w2-w1)/(w2-r),miu,sigma); % 正态分布反函数
    v1 = normpdf(v1,miu,sigma); % 正态分布概率密度函数
    fun1 = @(x) (x-v1).*normpdf(x,0,1);  
    L1 = quadl(fun1,v1,10000);  % 无穷积分,上限用10000
    E1 = (p-w1)*miu-(w1-r)*v1*sigma-(p-r)*L1*sigma; % 期望
   
    v2 = norminv((p-w2)/(p-r),miu,sigma);
    v2 = normpdf(v2,miu,sigma);
    fun2 = @(x) (x-v2).*normpdf(x,0,1);
    L2 = quadl(fun2,v2,10000);
    E2 = (p-w2)*miu-(w2-r)*v2*sigma-(p-r)*L2*sigma;
    plot(w1,E1,'r*',w1,E2,'o')
    hold on
end

附上参考：http://www.