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

Lorentz JÄNTSCHI (principal investigator) & Sorana D. BOLBOACĂ (co-investigator)
Calculul funcţiei erf(·) implică uneori şi calculul funcţiei complementare erfc(·), care este
dată de formula:
erfc(x) = ∫ e
32(157)
∞
2 2
−t
dt
π z
Calculul aproximativ al celor două funcţii se bazează pe polinoame de interpolare pe diferite
intervale pentru variabila transmisă ca argument. Astfel calculul este destul de laborios. Algoritmul
implementării celor două funcţii este prezentat mai jos:
function error(x) { //erf(·)
e_efx=1.28379167095512586316e-01;
ePp=array(1.28379167095512558561e-01,
-3.25042107247001499370e-01,-2.84817495755985104766e-02,
-5.77027029648944159157e-03,-2.37630166566501626084e-
05);
eQq=array(3.97917223959155352819e-01,
6.50222499887672944485e-02,5.08130628187576562776e-03,
1.32494738004321644526e-04,-3.96022827877536812320e-06);
//Coefficients for approximation to erf in[0.84375,1.25]
ePa=array(-2.36211856075265944077e-03,
4.14856118683748331666e-01,-3.72207876035701323847e-01,
3.18346619901161753674e-01,-1.10894694282396677476e-01,
3.54783043256182359371e-02,-2.16637559486879084300e-03);
eQa=array(1.06420880400844228286e-01,
5.40397917702171048937e-01,7.18286544141962662868e-02,
1.26171219808761642112e-01,1.36370839120290507362e-02,
1.19844998467991074170e-02);
e_erx = 8.45062911510467529297e-01;
abs_x = (x >= 0.0 ? x : -x);
if(abs_x