請問這個問題的lingo編程是什麽啊，謝謝啦！

sets:

quarters/1..4/:x,q,a,p,b;

endsets

data:

q=1000,800,800,1000;

a=10,15,25,20;

p=0.00065,0.00065,0.00065,0.00065;

b=25,35,30,10;

c=0.15;

enddata

min=@sum(quarters:x*q*p)+c*(3*(x(1)-a(1))+2*(x(2)-a(2))+(x(3)-a(3)))+0.00645372-0.0668;

x(1)>=a(1);

x(1)+x(2)-b(1)>=a(2);

x(1)+x(2)+x(3)-b(1)-b(2)>=a(3);

x(1)+x(2)+x(3)+x(4)-b(1)-b(2)-b(3)>=a(4);

x(1)+x(2)+x(3)+x(4)=a(1)+a(2)+a(3)+a(4);

end