钢管订购和运输程序实现的程序

T=[104 301 750 606 194 205 201 680 480 300 220 210 420 500];
S=[800 800 1000 2000 2000 2000 2000];
C=[170.7 160.3 140.2 98.6 38 20.5 3.1 21.2 64.2 92 96 106 121.2 128 142;
215.7 205.3 190.2 171.6 111 95.5 86 71.2 114.2 142 146 156 171.2 178 192;
230.7 220.3 200.2 181.6 121 105.5 96 86.2 48.2 82 86 96 111.2 118 132;
260.7 250.3 235.2 216.6 156 140.5 131 116.2 84.2 62 51 61 76.2 83 97;
255.7 245.3 225.2 206.6 146 130.5 121 111.2 79.2 57 33 51 71.2 73 87;
265.7 255.3 235.2 216.6 156 140.5 131 121.2 84.2 62 51 45 26.2 11 28;
275.7 265.3 245.2 226.6 166 150.5 141 131.2 99.2 76 66 56 38.2 26 2];
P=[160 155 155 160 155 150 160]
P=kron(P',ones(1,15));
D=[C+P;0.05*ones(2,15)]';
f=D(:);
H=[zeros(105,105) zeros(105,30);
  zeros(30,105) 0.1*eye(30,30)];
A1=kron(eye(7,9),ones(1,15));
b1= S';
temp=eye(15);
e1=[zeros(1,15);temp(1,:)];
el5=[temp(15,:);zeros(1,15)];
A2=[zeros(2,105) el5 e1]
b2=zeros(2,1);
A=[A1;A2];
b=[b1;b2];
Aeq1=kron([ones(1,7) -ones(1,2)],eye(15));
beq1=zeros(15,1);
Aeq2=[zeros(14,105) temp(1:14,:) temp(2:15,:)];
beq2=T';
Aeq=[Aeq1;Aeq2];
beq=[beq1;beq2];