+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Software Libraries (/forum-10.html)
+--- Forum: General Software Library (/forum-13.html)
+--- Thread: (71B) Buffon's Needle (/thread-13610.html)
(71B) Buffon's Needle - SlideRule - 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
RE: (71B) Buffon's Needle - Valentin Albillo - 09-08-2019 09:54 PM
(09-08-2019 04:00 PM)SlideRule Wrote: 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
You will get a syntax error entering line 30, it's just RND, without the parentheses, not RND().
Also, it should run somewhat faster if you include the NEXT J also in line 30, like this:
30 H=0 @ FOR J=1 TO 10 @ H=H+(RND<SIN(RND*PI)) @ NEXT J
(written and posted on the go using a tablet, still not at my laptop)
Regards.
V.
.
RE: (71B) Buffon's Needle - J-F Garnier - 09-09-2019 07:20 AM
(09-08-2019 09:54 PM)Valentin Albillo Wrote:(09-08-2019 04:00 PM)SlideRule Wrote: 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
You will get a syntax error entering line 30, it's just RND, without the parentheses, not RND().
Actually, RND(), with parentheses, is a valid syntax on the HP71, as for other functions without argument. Even PI can be entered as PI() since PI is a function from the HP71 system point of view.
J-F