lauricella's hypergeometric function

6 次查看（过去 30 天）

显示更早的评论

mutah 2013-11-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/108066-lauricella-s-hypergeometric-function

评论： Roger Stafford 2013-11-29

在 MATLAB Online 中打开

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 个评论
显示 -1更早的评论隐藏 -1更早的评论

Roger Stafford 2013-11-29

在 MATLAB Online 中打开

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?

请先登录，再进行评论。

请先登录，再回答此问题。