0

重複の可能性:
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;

私が犯している明らかな間違いはありますか?

4

0 に答える 0