19 lines
309 B
Text
19 lines
309 B
Text
include ..\Utilitys.pmt
|
|
|
|
def RandomNormal
|
|
drop rand log -2 * sqrt 2 pi * rand * cos * 0.5 * 1 +
|
|
enddef
|
|
|
|
1000 var n
|
|
0 n repeat
|
|
|
|
getid RandomNormal map
|
|
|
|
dup
|
|
sum n / var mean
|
|
"Mean: " print mean print nl
|
|
|
|
0 swap n for
|
|
get mean - 2 power rot + swap
|
|
endfor
|
|
swap n / sqrt "Standard deviation: " print print
|