lauricella's hypergeometric function

10 vues (au cours des 30 derniers jours)
mutah
mutah le 29 Nov 2013
Commenté : Roger Stafford le 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 commentaire
Roger Stafford
Roger Stafford le 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?

Connectez-vous pour commenter.

Réponses (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by