2013西北大学数学建模竞赛(陈思、李瑶、张瑜)(8)

end y(m,n)=1; y=y';

for i=1:citynum for j=1:citynum

mubiaob=c(i,j)*y(i,:); end end xuq1=0;

for i=1:citynum for j=1:citynum

xuq1=xuq1+s(i)*y(i,:)-q(i); end

xuqiu=max((xuq1),0)*M; end end shij1=0; shij2=0; for i=1:cir for j=1:cir for l=1:cir

shij1=shij1+t(i)-a(i); shij2=shij12+b(i)-t(i); end

shij3=max((shij1),0); shij4=max((shij2),0); shijian=M*shij3+M*shij4; end end

f=mubiao+xuqiu+shijian; f=1/f; end end }

{ int fsum=0; for i=1:inn fsum=fsum+f(i);