MATLAB Answers

mutah
0

lauricella's hypergeometric function

Asked by mutah
on 29 Nov 2013
Latest activity Commented on by Roger Stafford 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);
FA = triplequad(Q,0,1,0,1,0,1)
end
any one help ??

  1 Comment

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?

Sign in to comment.

0 Answers