モンテカルロ法だけど、
図は視覚的なものほしいからけずらないとして
これで１００００くらいやると、すごく時間かかる
改良してちょ

10 cls :r=0
20 input"k=";k:cls
30 line(0,0)-(200,200),3,b
40 circle(0,0),200,3
50 for i=1 to k
60 x=rnd*1:y=rnd*1
70 n=x*x+y*y:c=3
80 if n<=1 then 90 else 100
90 r=r+1:c=5
100 circle(x*200,y*200),1,c
110 next i
120 pi=4*r/k
130 locate 20,20:print "r=";r,"pi=";pi
140 end