戻る
■ 円周率の近似値を求める手法 2018. 3.26

高機能電卓の情報
https://akatuki-724.blogspot.jp/2018/03/blog-post.html?showComment=1521968112039&m=1
の記事により、プチコンで組みました。


D=0.0000001
P=ATAN(1)*4
FOR J=1 TO 10000
 I=FLOOR(J*P+.5)
 IF ABS(I/J-P) <= D THEN
  S=I:T=J
  WHILE S MOD T != 0
   R=S MOD T
   S=T:T=R
  WEND
  IF R==1 THEN PRINT I,J,I/J ELSE PRINT "..." ENDIF
 ENDIF
NEXT J
END