Symbol Error Probability of DF Relay Selection over Arbitrary Nakagami-𝑚 Fading Channels

Algorithm 1

MATLAB Program for evaluating 𝐹₂(⋅).

Function  F2 = Appell(a,b1,b2,c1,c2,x,y);

f1 = gamma(c1).*gamma(c2);

f2 = gamma(b1).*gamma(b2);

f3 = gamma(c1-b1).*gamma(c2-b2);

f = f1./(f2.*f3);

Q = @(u,v)f.*u.̂‌(b1-1).*(v.̂‌(b2-1)).*...

   ((1-u).̂(c1-b1-1)).*...

   ((1-v).̂(c2-b2-1))./...

   ((1-u.*x-v.*y).̂a);

F2 = dblquad(Q,0,1,0,1);