円周率の算出

モンテカルロ法を使った π の算出

興味を持ってもらうため
JIS Full BASIC準拠の十進BASICで作ってみました。

100 ! 'モンテカルロ法によるπの算出
110 ! '
120 SET WINDOW -1, 4, -1, 4
130 RANDOMIZE
140 INPUT PROMPT "データ数 ":N
150 LET COUNT = 0
160 SET POINT STYLE 3
170 FOR I=1 TO N STEP 1
180    LET X = RND
190    LET Y = RND
200    IF X^2+Y^2 <1 THEN
210       SET POINT COLOR 4
220       LET COUNT = COUNT+1
230    ELSE
240       SET POINT COLOR 2
250    END IF
260    PLOT POINTS : X*3, Y*3
270 NEXT I
280 LET PAI = COUNT*4/N
290 PRINT "πの近似値は　";PAI
300 END