07/11/09 16:16:37
>>143
//ボックス＝ミューラー法で100個の正規乱数を生成するテスト
#include<iostream>
#include<math.h>
#include<stdio.h>
#include<time.h>
int main()
{
int i;
double u1,u2,n1,n2;
srand(time(NULL));
for(i=0;i<50;i++){
u1=(double)(rand())/RAND_MAX;
u2=(double)(rand())/RAND_MAX;
n1=sqrt(-2*log(u1))*sin(2*M_PI*u2);
n2=sqrt(-2*log(u1))*cos(2*M_PI*u2);
std::cout <<n1<<" "<<n2<<std::endl;
}
return 0;
}