**parameter**

numb(m)

G(m,n)

PDGmax(m)

;

**binary variable**

zz(m,t)

;

**variable**

Pm(m,MEG,t)

;

e18(m,t)$(PDGmax(m) ne 0 or

**Pm(m,MEG,t) ne 0**).. zz(m,t)=e=1;

e19(m,t)$(PDGmax(m) eq 0 and

**Pm(m,MEG,t) eq 0**).. sum(n$(ord(m) ne ord(n) and G(m,n) ne 0),delta(m,n,t))/numb(m)=l=zz(m,t);

e20(m,t)$(PDGmax(m) eq 0 and

**Pm(m,MEG,t) eq 0**) .. sum(n$(ord(m) ne ord(n) and G(m,n) ne 0),delta(m,n,t))=g=zz(m,t);