今天做实验需要用蒙特卡罗法产生满足正态分布的随机数，下面是我用C语言编写的程序：(转载请注明出处)

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#include<time.h>
double random(void)
{
 double r;
 r=rand()/32767.00;
 return r;
}
int sta(double mu,double sigma)
{
 double r1,r2;
 double pi=3.1415926;
 r1=random();
 r2=random();
 return (int)(sqrt(-2*log(r1))*cos(2*pi*r2)*sigma+mu);
}
int main()
{
 int i,n;
 double mu,sigma;
 FILE *fp;
 fp=fopen("data.dat","w");
 srand((unsigned)time(NULL));
 printf("please input mu:");
 scanf("%lf",&mu);
 printf("please input sigma:");
 scanf("%lf",&sigma);
 printf("please input n:");
 scanf("%d",&n);
 printf("Random number:\n");
 for(i=0;i<n;i++)
 {
  printf("%d\n",sta(mu,sigma));
  fprintf(fp,"%d\n",sta(mu,sigma));
 }
 fclose(fp);
 return 0;
}