Distribuţia Binomială: Modelare Statistică, Optimizare Numerică, cu ...

Lorentz JÄNTSCHI (principal investigator) & Sorana D. BOLBOACĂ (co-investigator)
define("EPS", 2.22e-16);
function logGamma(x){
//LogGammarelatedconstants
lg_d1=-0.5772156649015328605195174;
lg_d2=0.4227843350984671393993777;
lg_d4=1.791759469228055000094023;
lg_p1=array(4.945235359296727046734888,
201.8112620856775083915565,2290.838373831346393026739,
11319.67205903380828685045,28557.24635671635335736389,
38484.96228443793359990269,26377.48787624195437963534,
7225.813979700288197698961);
lg_p2=array(4.974607845568932035012064,
542.4138599891070494101986,15506.93864978364947665077,
184793.2904445632425417223,1088204.76946882876749847,
3338152.967987029735917223,5106661.678927352456275255,
3074109.054850539556250927);
lg_p4=array(14745.02166059939948905062,
2426813.369486704502836312,121475557.4045093227939592,
2663432449.630976949898078,29403789566.34553899906876,
170266573776.5398868392998,492612579337.743088758812,
560625185622.3951465078242);
lg_q1=array(67.48212550303777196073036,
1113.332393857199323513008,7738.757056935398733233834,
27639.87074403340708898585,54993.10206226157329794414,
61611.22180066002127833352,36351.27591501940507276287,
8785.536302431013170870835);
lg_q2=array(183.0328399370592604055942,
7765.049321445005871323047,133190.3827966074194402448,
1136705.821321969608938755,5267964.117437946917577538,
13467014.54311101692290052,17827365.30353274213975932,
9533095.591844353613395747);
lg_q4=array(2690.530175870899333379843,
639388.5654300092398984238,41355999.30241388052042842,
1120872109.61614794137657,14886137286.78813811542398,
101680358627.2438228077304,341747634550.7377132798597,
446315818741.9713286462081);
lg_c=array(-0.001910444077728,8.4171387781295e-4,
-5.952379913043012e-4,7.93650793500350248e-4,
-
0.002777777777777681622553,0.08333333333333333331554247,
0.0057083835261);
//Rough estimate of the fourth root of logGamma_xBig
lg_frtbig=2.25e76; pnt68=0.6796875;
if(x==this->logGammaCache_x){
return this->logGammaCache_res;
}
y=x;
if(y>0.0&&y