home *** CD-ROM | disk | FTP | other *** search
- 100 REM subroutine RNDN.BAS - to calculate random normal variables x=0 var=1
- 105 '
- 110 V1.RND = (2!*RND(1))-1!
- 120 V2.RND = (2!*RND(1))-1!
- 130 S.RND = (V1.RND*V1.RND) + (V2.RND*V2.RND)
- 140 IF S.RND > 1! THEN 110
- 150 X1.RND = V1.RND * SQR(-2!*LOG(S.RND)/S.RND)
- 160 X2.RND = V2.RND * SQR(-2!*LOG(S.RND)/S.RND)
- 170 ' x1.rnd and x2.rnd are random normals
- RND * SQR(-2!*LOG(S.R