![](https://csdnimg.cn/release/download_crawler_static/85912575/bg1.jpg)
Min=100000000;
for n=5:1:50
for m=200:1:600
A=0; B=0; L=fix(m/n);
for i=1:L
A=A+(normcdf((i*n),600,196.6292)-normcdf(((i-1)*n),600,196.6292))*(10*i+100*n+3000);
B=B+(i-0.5)*n* (normcdf((i*n),600,196.6292)-normcdf(((i-1)*n),600,196.6292));
end
Q=(A+(normcdf(m,600,196.6292)-normcdf((L*n),600,196.6292))*(10*L+100*(m-n*L)+3000)+(
1-normcdf(m,600,196.6292))*(10*L+1000))/(B+(normcdf(m,600,196.6292)-normcdf((L*n),600,1
96.6292))*(m+n*L)/2+(m*(1-normcdf(m,600,196.6292))))+5/19;
if Q<Min
Min=Q;
J=n;
T=m;
end
end
end