男女の確率問題(3)の解答例を示します。
(3)の解
(考察)
問題(3)のシミュレーションを行うプログラムを書き、実行した。
●プログラム(BG311.bas)
' << BG311.bas >> ' Tiny Basic ' 問題(3)のシミュレーション ' Dim A(2): ' 男女の状態を表す配列。 ' A(1)=0:1番目が男の場合。 ' A(1)=1:1番目が女の場合。 ' A(2)=0:2番目が男の場合。 ' A(2)=1:2番目が女の場合。 Dim B(3): ' 2人の子どもの状態を表す配列。 ' B(0):1番目が男、2番目が男の場合の数。 ' B(1):1番目が男、2番目が女の場合の数。 ' B(2):1番目が女、2番目が男の場合の数。 ' B(3):1番目が女、2番目が女の場合の数。 ' Do ' 実験回数NMAXを読み込む。 Read NMAX If NMAX <= 0 Then Exit Do ' ' 初期設定。 Randomize 123: ' 乱数の初期設定。 For K=0 to 3: B(K)=0: Next K COUNT=0: ' 『少なくとも1人が男の子』となる場合の数。 ' ' NMAX回実験を繰り返す。 For N=1 To NMAX ' 1番目の子供の判定。 ' 1番目の子供の男女の出生割合を、1:1と考える。 ' 1番目が男ならA(1)=0、女ならA(1)=1。 If Rnd < 0.5 Then A(1)=0 Else A(1)=1 ' 2番目の子供の判定。 ' 2番目の子供の男女の出生割合を、1:1と考える。 ' 2番目が男ならA(2)=0、女ならA(2)=1。 If Rnd < 0.5 Then A(2)=0 Else A(2)=1 ' ' 『少なくとも1人は男の子』という条件の下で考察する。 If (A(1) = 0) or (A(2) = 0) Then K=A(1)+2*A(2) COUNT=COUNT+1 ' 配列Bに記録する。 ' A(1)=0,A(2)=0 のとき、K=0 ' A(1)=1,A(2)=0 のとき、K=1 ' A(1)=0,A(2)=1 のとき、K=2 ' A(1)=1,A(2)=1 のとき、K=3 B(K)=B(K)+1 End If Next N ' ' 実験結果の表示。 Print"実験回数:";NMAX Print"2番目の子が男の子となる確率:"; Print Using"#.###";(B(0)+B(1))/COUNT Print Loop End ' ' データ。 Data 1000, 10000, 100000, 1000000, 0
実行結果