RosettaCodeData/Task/Monte-Carlo-methods/HicEst/monte-carlo-methods.hicest
Ingy döt Net db842d013d A-M baby
2013-04-10 21:29:02 -07:00

12 lines
311 B
Text

FUNCTION Pi(samples)
inside = 0
DO i = 1, samples
inside = inside + ( (RAN(1)^2 + RAN(1)^2)^0.5 <= 1)
ENDDO
Pi = 4 * inside / samples
END
WRITE(ClipBoard) Pi(1E4) ! 3.1504
WRITE(ClipBoard) Pi(1E5) ! 3.14204
WRITE(ClipBoard) Pi(1E6) ! 3.141672
WRITE(ClipBoard) Pi(1E7) ! 3.1412856