personal/time-to-botec
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix the base of the log for squiggle.bc

Browse Source
This commit is contained in:
NunoSempere 2023-11-02 23:26:06 +00:00
parent 249a1ff434
commit 8675d98784
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
bc/squiggle.bc
View File

@ -13,7 +13,7 @@ define sample_unit_uniform(){
define sample_unit_normal(){ define sample_unit_normal(){
u1=sample_unit_uniform() u1=sample_unit_uniform()
u2=sample_unit_uniform() u2=sample_unit_uniform()
z = sqrt(-2 * log(u1, 2)) * sin(2 * pi * u2) z = sqrt(-2 * l(u1)) * sin(2 * pi * u2)
return z return z
} }
@ -66,7 +66,7 @@ define sample_to(low, high){
we need but get a normal with 90% confidence interval [log(a), log(b)]. we need but get a normal with 90% confidence interval [log(a), log(b)].
Then see code for sample_normal_from_90_confidence_interval Then see code for sample_normal_from_90_confidence_interval
*/ */
loglow = log(low, 2) loglow = l(low)
loghigh = log(high, 2) loghigh = l(high)
return e(sample_normal_from_90_confidence_interval(loglow, loghigh)) return e(sample_normal_from_90_confidence_interval(loglow, loghigh))
} }