09-08-2019, 04:00 PM

An extract from Buffon, [kiyoshiakima url]

"2 Simulation Programs

Buffon's Needle can be implemented easily on all but the smallest HP programmable …

2.2 BASIC

For the HP-71B the following BASIC program simulates the experiment. It prompts for the number of ten-throw trials, then prints the occurrences of each outcome and an estimate for π.

10 DESTROY ALL @ OPTION BASE 0 @ RADIANS @ DIM R(11) @ P=0

20 INPUT N @ FOR I=1 TO N

30 H=0 @ FOR J=1 TO 10 @ IF RND()<SIN(RND()*PI) THEN H=H+1

40 NEXT J @ R(H)=R(H)+1 @ P=P+H

50 NEXT I

60 FOR I=0 TO 10 @ PRINT I;R(I) @ NEXT I

70 PRINT 20*N/P"

BEST!

SlideRule

"2 Simulation Programs

Buffon's Needle can be implemented easily on all but the smallest HP programmable …

2.2 BASIC

For the HP-71B the following BASIC program simulates the experiment. It prompts for the number of ten-throw trials, then prints the occurrences of each outcome and an estimate for π.

10 DESTROY ALL @ OPTION BASE 0 @ RADIANS @ DIM R(11) @ P=0

20 INPUT N @ FOR I=1 TO N

30 H=0 @ FOR J=1 TO 10 @ IF RND()<SIN(RND()*PI) THEN H=H+1

40 NEXT J @ R(H)=R(H)+1 @ P=P+H

50 NEXT I

60 FOR I=0 TO 10 @ PRINT I;R(I) @ NEXT I

70 PRINT 20*N/P"

BEST!

SlideRule