lauricella's hypergeometric function

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);

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

end

any one help ??

[Roger Stafford]

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?