[飞逝1991]

a=0;
while(1.1-a)>1
c=[-0.27 -0.19 -0.185 -0.185];
Aeq=[1.01 1.02 1.045 1.065];
beq=[0.2];
A=[0.00125 0 0 0;0 0.00075 0 0;0 0 0.00275 0;0 0 0 0.0013];
b=[a;a;a;a];
vlb=[0,0,0,0];
vub=[];
[x,val]=linprog(c,A,b,Aeq,beq,vlb,vub);
a;
x=x';
Q=-val;
plot(a,Q,'.'),axis([0 0.2 0 6]), hold on
a=a+0.01;
end
xlabel('a'),ylabel('Q')