How to correctly express this expression（20）? And the code I wrote as follows,is there an error?.
pjam is the variable， other parameters are known.
rie_ub=[rie_ub1;rie_ub2]；r0 is a constant，hie_max=[ hie_max1 ; hie_max2 ]，
hke_min=hie_min=[ hie_min1 ; hie_min2 ]，heu_min=[ heu_min(1) ; heu_min(2) ]
rie_ub1(1,8)， hie_max1(1,8)， hie_min1(1，8)，pi=pii_f
/////////////////////////////////////////////////////////////////////////////////////////
C1=[];C2=[];A1=[];A2=[];
for i=1:8
C1=[C1 pii_f(i) * hie_max1(i)]; %(1,8)
C2=[C2 pii_f(i) * hie_max2(i)]; %(1,8)
end
for i=1:8
h1=hie_min1; h2=hie_min2;
h1(i)=0; h2(i)=0;
A1=[A1 pii_f * h1’+sigma]; %(1,8)
A2=[A2 pii_f * h2’+sigma]; %(1,8)
end
for i=1:8
X1(i)=A1(i) / C1(i)+heu_min(1) * pjam / C1(i);
X2(i)=A2(i) / C2(i)+heu_min(2) * pjam / C2(i);
rie_ub1(i)=(rel_entr(1,X1(i))+rel_entr(X1(i),1)) / log(2);
rie_ub2(i)=(rel_entr(1,X2(i))+rel_entr(X2(i),1)) / log(2);
end
rie_ub=[rie_ub1 ; rie_ub2];
////////////////////////////////////////////////////////////////////////////////////////////
subject to
for i=1:8
if a_f(i)==1
rie_ub(1,i)<=r0;
rie_ub(2,i)<=r0;
end
end
///////////////////////////////////////////////////////////////////////////////////////