| # | Line 138 | Line 138 | un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/ | |
|---|---|---|
| 138 | (1+t*(1.432788+t*(.189269+t*.001308))) ; | |
| 139 | un1`private(p) : un2`private(sqrt(-2*log(p))) ; | |
| 140 | ||
| 141 | < | unif2norm(p) : if( .5-p, un1`private(p), -un1`private(1-p) ) ; |
| 141 | > | unif2norm(p) : if( .5-p, -un1`private(p), un1`private(1-p) ) ; |
| 142 | ||
| 143 | nrand(x) = unif2norm(rand(x)); | |
| 144 | ||
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |