[1] 授業単元：計算実習
[2] 問題文(含コード&リンク)：
乗算合同法を用いて擬似乱数を作成し、Box-Muller法により平均０、分散１の正規乱数系列に
変換した場合に（1×10^6）個の系列から以下の方法で確率密度関数を推定し、ガウス分布
p(x) = {1 / (2π)^(1/2) } × exp(-x^2 / 2)
と比較せよ。
〜確率密度関数の推定〜
確率密度関数p(x)とした場合に確率密度関数のヒストグラムからの推定値は以下で与えられる。
p(x_i) = {N(x_i - Δ/2 ; x_i + Δ/2)} / MΔ
ここで、N(x_1 ; x_2)は区間(x_1 , x_2] に含まれる個数を表し、Ｍはサンプル数、Δはヒストグラムの
刻み幅である。
[3] 環境
　[3.1] OS： Windows XP
　[3.2] コンパイラ名とバージョン： gcc 3.4
　[3.3] 言語： C
[4] 期限： 2007年5月28日
[5] その他の制限： 特にありません

よろしくお願いします。