| 66 |
|
bound(a,x,b) : if(a-x, a, if(x-b, b, x)); |
| 67 |
|
Acos(x) : acos(bound(-1,x,1)); |
| 68 |
|
Asin(x) : asin(bound(-1,x,1)); |
| 69 |
< |
Exp(x) : if(-x-60, 0, exp(x)); |
| 69 |
> |
Exp(x) : if(-x-100, 0, exp(x)); |
| 70 |
|
Sqrt(x) : if(x, sqrt(x), 0); |
| 71 |
|
|
| 72 |
|
{ Useful constants } |
| 135 |
|
|
| 136 |
|
un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/ |
| 137 |
|
(1+t*(1.432788+t*(.189269+t*.001308))) ; |
| 138 |
< |
un1`private(p) : un2`private(sqrt(log(1/p/p))) ; |
| 138 |
> |
un1`private(p) : un2`private(sqrt(-2*log(p))) ; |
| 139 |
|
|
| 140 |
|
unif2norm(p) : if( .5-p, un1`private(p), -un1`private(1-p) ) ; |
| 141 |
|
|
