C语言生成正态分布（标准分布）数

#include<stdio.h>
#include<math.h>
#include <stdlib.h>

int meandata=0;
int variance=1;
#define PI 3.141592654
#define naturale 2.7182818284

double gaussrand(double x)       //高斯分布函数
{

    double U, Z;
    U=-pow(x-meandata,2)/2*pow(variance,2);
//    printf("U:  %f\n",U);fsd
    Z=pow(naturale,U)*(1/(sqrt(2*PI)*variance));
 //   printf("Z:  %f\n",Z);
    return Z;
}


double GAUSSRAND(int x)
{
    double s1 = 0;
    float i;      // 循环变量
    for(i=x-10; i<=x;i=i+0.1)
    {
        s1 = s1+gaussrand(i)*0.1;

        printf("%f\n",s1*72*640);
    }
    return s1;
}


int main()
	{           
		double s= GAUSSRAND(0);
		printf("%f\n",s);
	}