PRINT
PRINT "１〜１０のカードが１万枚ずつある。"
PRINT "カードを任意に３０枚取ったとき、"
PRINT "１〜１０のカードがすべてそろっている確率は？"
PRINT
PRINT "１〜１０のカードが１万枚ずつある。"
PRINT "カードを任意に１枚ずつ取っていったとき、"
PRINT "１〜１０のカードがすべてそろうためには、"
PRINT "平均何枚のカードを取らなければならないか？"
PRINT
FUNCTION RND10(q)     ! 1〜10の整数をランダムに出す。重複あり。
   LET X=INT(RND*10)+1
   LET RND10=X
END FUNCTION

FUNCTION HEIKIN(A(),N)
   LET  s=A(1)
   FOR i=2 TO N
      LET  s=s+A(i)
   NEXT i
   LET  HEIKIN=s/N
END FUNCTION

RANDOMIZE
PRINT
PRINT "１回目から２０回目までのシミュレーションの結果："
PRINT
PRINT "   ";
DIM Q(100),R(100),T(1000)
MAT R=ZER
LET U=0
FOR c=1 TO 1000
   FOR i=1 TO 100
      LET Q(i)=RND10(i)
      LET x=Q(i)
      IF c<21 THEN
         PRINT x;
      END IF
      LET R(x)=1
      IF R(1)*R(2)*R(3)*R(4)*R(5)*R(6)*R(7)*R(8)*R(9)*R(10)>0 THEN
         LET T(c)=i
         IF T(c)<31 THEN
            LET U=U+1
         END IF
         EXIT FOR
      END IF
   NEXT i
   MAT Q=ZER
   MAT R=ZER
   IF c<21 THEN
      PRINT
      PRINT
      PRINT "   ";
   END IF
NEXT c
PRINT
PRINT "   全部そろうまでに集めたカードの枚数："
PRINT
PRINT "       ";
FOR i=1 TO 20
   PRINT T(i);
NEXT i
PRINT
PRINT
LET Z=HEIKIN(T,1000)
PRINT
PRINT
PRINT "１０００回のシミュレーションの結果："
PRINT
PRINT "   カードを任意に３０枚取ったとき、１〜１０のカードがすべてそろっている確率は、"
PRINT "        ";U/10;"％ です。"
PRINT
PRINT "   平均";Z;"枚 のカードを集めると、１〜１０すべてのカードがそろいます。"
PRINT
END