File : hark-data.inc used by : hark-monop.gms [html] hark-oligop.gms [html]
Sets
n nodes / one, two, three, four, five, six /
L(n) regions / one, two, three /
NL(n) 'complement of L'
;
NL(n) = not(L(n));
Alias (n,i,j);
Alias (L,Q,QQ,Qi);
Set Q1(Q);
Table coefs(L,*) demand and supply data
alpha beta rho eta
one 1.0 .5 19 .2
two 2.0 .4 27 .01
three 1.5 .3 30 .3 ;
Table pairs(i,j,*) transport data
kappa nu
one.four 1 .5
one.five 2 .2
two.six 3 .3
three.six 1 .4
four.one 2 .3
four.five 1 .1
four.six 1 .1
five.one 3 .5
five.four 2 .2
five.six 1 1.0
six.two 2 .25
six.three 2 .2
six.four 1 .9
six.five 3 .8 ;
parameter
alpha(L), beta(L), rho(L), eta(L)
kappa(i,j), nu(i,j);
alpha(L) = coefs(L,'alpha');
beta(L) = coefs(L,'beta');
rho(L) = coefs(L,'rho');
eta(L) = coefs(L,'eta');
kappa(i,j) = pairs(i,j,'kappa');
nu(i,j) = pairs(i,j,'nu');
Set arc(i,j);
arc(i,j) = yes$kappa(i,j);
$macro repMonop(style) rep1(i,j,"style") = t.l(i,j); rep2("supply",L,"style") = s.l(L); rep2("demand",L,"style") = d.l(L); rep2("price ",L,"style") = price.l(L); tab6("totProfit","style") = obj.l; tab6("totSupply","style") = sum{L, s.l(L)}; tab6("theta1","style") = price.l('one'); tab6("theta2","style") = price.l('two'); tab6("theta3","style") = price.l('three');