#include #include #include #include using namespace std; #define X_MIN -5 #define X_MAX 5 #define RAND_NUM 10000 double gaus(double x){ // return exp(-x*x/2.)/sqrt(2*M_PI); return exp(-x*x/2.);//normalize gaus(0)=1 } int main(int argc, char *argv[] ) { int i,bin; double x,r; srand(time(NULL)); i=0; while(RAND_NUM>i){ x=(X_MAX-X_MIN)*((double)rand()/RAND_MAX-0.5); if((double)rand()/RAND_MAX