重複の可能性:
C#でのBox-Mueller乱数ジェネレーターの実装
Box-Muller変換法を使用して、正規分布の乱数を生成しています。
do {
rand1 = 2*(rand()/(double)RAND_MAX)-1 ;
rand2 = 2*(rand()/(double)RAND_MAX)-1;
distSqr = rand1 * rand1 + rand2 * rand2;
} while (distSqr >= 1 );
mp1 = sqrt((-2*log(distSqr))/distSqr);
res = rand1 * mp1;
norm_dist=res*var+mean;
私が犯している明らかな間違いはありますか?