Displaying 1 result from an estimated 1 matches for "integpsvbar".
2006 Sep 01
1
integration problem with gamma function
...vbar = NIntegrate[1/(UiA^(Ni/2)) 1/(UiiA^(Nii/2))
Gamma[Ni/2,UiA/(sH^2),UiA/(sL^2)]
Gamma[Nii/2,UiiA/(sH^2),UiiA/(sL^2)],{A,L,H},
MinRecursion->3];
PSVbar = psvbar/(4 Log[sH/sL]);
Print["p(s?v|D_1D_2I) = const. ",N[PSVbar,6]];
</--->
<--- translation to R --->
integpsvbar <- function(A)
{
# Mathematica: gamma[a,z0,z1] = gamma[a,z1] - gamma[a,z0]
UiA <- Ni*(Dsi-2*Di*A+A^2)/2
UiiA <- Nii*(Dsii-2*Dii*A+A^2)/2
1/UiA^(Ni/2) *
1/UiiA^(Nii/2) *
( pgamma(UiA/(sL^2), Ni/2) - pgamma(UiA/(sH^2), Ni/2) ) *
( pgamma(UiiA/(sL^2), Nii/2) - pgamma(UiiA/(sH^2), Nii/2)...