model:
sets:
owall/1..4/;
nwall/1..8/;
year/1..5/;
cji(year,owall):c;
bik(year,nwall):b,m,x;
endsets
min=@sum(bik:m*x)+0.66*(Q^0.51)*L;
Q>=100;
@sum(owall(j):c(1,j))+@sum(nwall(k):b(1,k)*x(1,k))>=150;
@sum(owall(j):c(2,j))+@sum(nwall(k):b(1,k)*x(1,k))*0.9+@sum(nwall(k):b(2,k)*x(2,k))>=160;
@sum(owall(j):c(3,j))+@sum(nwall(k):b(1,k)*x(1,k))*0.9^2+@sum(nwall(k):b(2,k)*x(2,k))*0.9+@sum(nwall(k):b(3,k)*x(3,k))>=170;
@sum(owall(j):c(4,j))+@sum(nwall(k):b(1,k)*x(1,k))*0.9^3+@sum(nwall(k):b(2,k)*x(2,k))*0.9^2+@sum(nwall(k):b(3,k)*x(3,k))*0.9+Q>=180;
@sum(owall(j):c(5,j))+@sum(nwall(k):b(1,k)*x(1,k))*0.9^4+@sum(nwall(k):b(2,k)*x(2,k))*0.9^3+@sum(nwall(k):b(3,k)*x(3,k))*0.9^2+Q>=190;
0.22*Q^0.51*L+@sum(nwall(k):m(1,k)*x(1,k))