model:
sets:
liu/1..5/:r,p,x,q;
wen/1..100/:shu,ju;
endsets
max=obj;
[email=obj=@sum(liu:(r-p)*x]obj=@sum(liu:(r-p)*x[/email]);
@for(liu:q*x<=a*m);
@sum(liu:(1+p)*x)=m;
data:
r=0.05 0.25 0.23 0.28 0.21;
q=0 0.025 0.015 0.055 0.025;
p=0 0.01 0.02 0.045 0.065;
a0=0.001;a=0;
m=1;
enddata
calc:
@SET('TERSEO', 2);
@for(wen(i):
   @solve();
   shu(i)=obj;
   ju(i)=a;
   a=a+a0;   
   );
@ole('f:liu.xls',"y")=shu;
@ole('f:liu.xls',"x")=ju;
endcalc
end