c=[1.080;1.095;1.110;1.125;1.110;1.125;1.140;1.100;1.115;1.130];
A=[1 1 1 1 0 0 0 0 0 0;
0 0 0 0 1 1 1 0 0 0;
0 0 0 0 0 0 0 1 1 0;
0 0 0 0 0 0 0 0 0 1];
b=[25;35;30;10];
aeq=[1 0 0 0 0 0 0 0 0 0
0 1 0 0 1 0 0 0 0 0
0 0 1 0 0 1 0 1 0 0
0 0 0 1 0 0 1 0 1 1];
beq=[10;15;25;20];
vlb=[0;0;0;0;0;0;0;0;0;0];
[x,fval]=linprog(c,A,b,aeq,beq,vlb);
P=zeros(1,4);
P(1)=x(1)+x(2)+x(3)+x(4);P(2)=x(5)+x(6)+x(7);P(3)=x(8)+x(9);P(4)=x(10);
disp('各月生产的引擎数量为')
disp(P)
disp('总费用为')
C=1.08*x(1)+1.095*x(2)+1.110*x(3)+1.125*x(4)+1.110*x(5)+1.125*x(6)+1.14*x(7)+1.10*x(8)+1.115*x(9)+1.13*x(10);
disp(C)