srand((unsigned)time(NULL));
for(i = 0 ; i <= K ; i++){
r1 = rand()/(float)RAND_MAX;
r2 = rand()/(float)RAND_MAX;
u = sqrt(-2*log(r1));
v1 = u * cos(2.0 * PI * r2);
v2 = u * sin(2.0 * PI * r2);
fprintf(f1,"%lf\n",v1);
fprintf(f2,"%lf\n",v2);
��
评论0