model: sets: red/r1..r3/:r,ra;
time/1..3/:t,ta,k,e,ea,w,wa,m;
!inttime/i1..i3/:k;
!emt/et1..et3/:e,ea;
!wt/w1..w3/:w,wa;
!loopinter/lo1..lo3/:m;
endsets data:
r=0.5,0.5,0.5;
ra=0.5,0.5,0.5;
t=3,2,3;
ta=3,2,3;
k=0.1,0.1,0.1;
e=0.1,0.1,0.1;
enddata
max=b;
@for(time(i):w(i)+b<=1-r(i));
@for(time(i):wa(i)+ba<=1-ra(i));
@for(time(i):w(i)+a(i)-w(i+1)-wa(i+1)+t(i)+ta(i)+k(i)-k(i+1)=-0.5*(r(i)+ra(i))+0.5*(r(i+1)+ra(i+1))+e(i+1)+ea(i)+m(i));
ba=b;
end


time/1..3/:t,ta,k,e,ea,w,wa,m;
!inttime/i1..i3/:k;
!emt/et1..et3/:e,ea;
!wt/w1..w3/:w,wa;
!loopinter/lo1..lo3/:m;
endsets data:
r=0.5,0.5,0.5;
ra=0.5,0.5,0.5;
t=3,2,3;
ta=3,2,3;
k=0.1,0.1,0.1;
e=0.1,0.1,0.1;
enddata
max=b;
@for(time(i):w(i)+b<=1-r(i));
@for(time(i):wa(i)+ba<=1-ra(i));
@for(time(i):w(i)+a(i)-w(i+1)-wa(i+1)+t(i)+ta(i)+k(i)-k(i+1)=-0.5*(r(i)+ra(i))+0.5*(r(i+1)+ra(i+1))+e(i+1)+ea(i)+m(i));
ba=b;
end

