OPTION ARITHMETIC decimal_HIGH
FUNCTION CT(r)
   LET w=1
   FOR i=1 TO r
      LET w=w*i
   NEXT i
   LET CT=w
END FUNCTION

PRINT
PRINT "１人目が１人目と誕生日が同じである確率："
PRINT "          365/365"
PRINT "２人目が１人目と誕生日が異なる確率："
PRINT "          365/365 × 364/365"
PRINT "３人目が１人目と２人目と誕生日が異なる確率："
PRINT "          365/365 × 364/365 × 363/365"
PRINT
PRINT " n 人みんなの誕生日が異なる確率："
PRINT "        （ 365！ ÷ (365-n)！）÷ 365^n"
PRINT
PRINT "　何人集まれば、誕生日が一致する人がいる"
PRINT "確率が５０％以上になるか？"
PRINT
FOR n=20 TO 36
   LET s=CT(365)/(365^n*CT(365-n))
   PRINT n;"人 → 確率：";ROUND(1-s,3)
NEXT n
PRINT
PRINT "よって、２３人以上です。"
PRINT
END