## lauricella's hypergeometric function

### mutah (view profile)

on 29 Nov 2013
Latest activity Commented on by Roger Stafford

### Roger Stafford (view profile)

on 29 Nov 2013
I have this code but still found errors
function FA = Lauricella(a,b1,b2,b3,c1,x,y,z)
mu=0.5;
eta(1)=0.6;
h1=(2+eta(1)^-1+eta(1))/4;
H1=(eta(1)^-1-eta(1))/4;
SNR_avg=1;
a=1;b1=1;b2=mu;b3=mu;c1=2*mu+1.5;x=0.5;y=4*(h1-H1)*mu+SNR_avg/(4*(h1-H1)*mu+2*SNR_avg);
z=4*(h1+H1)*mu+SNR_avg/(4*(h1+H1)*mu+2*SNR_avg);
f1 = gamma(c1);
f2 = gamma(b1).*gamma(b2).*gamma(b3);
f3 = gamma(c1-b1).*gamma(c1-b2).*gamma(c1-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).^ (c1-b2-1)).*...
((1-w).^ (c1-b3-1))./...
((1-u.*x-v.*y-w.*z).^ a);
end
any one help ??

Roger Stafford

### Roger Stafford (view profile)

on 29 Nov 2013
I admit I am stumped on your expression for FA. It does not seem to match up with the one I see in Wikipedia at:
http://en.wikipedia.org/wiki/Lauricella_hypergeometric_series
In particular, FA in Wikipedia contains input arguments c2 and c3 which you do not have. Your input arguments appear to match those of FD. Can you please explain that?
Also please explain the nature of the error you are getting where you say "but still found errors". Is there an error message from matlab, or is it that the results are not as you expect?