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

Algorithm 2

MATLAB Program for evaluating 𝐹_𝐴⁽³⁾(⋅).

Function  FA = Lauricella  (a,b1,b2,b3,c1,c2,c3,x,y,z);

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

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

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

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

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

   (v.̂‌(b2-1)).*(w.̂‌(b3-1)).*...

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

   ((1-v).̂‌(c2-b2-1)).*...

   ((1-w).̂‌(c3-b3-1))./...

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

FA = triplequad(Q,0,1,0,1,0,1);