Моделирование процессов в физике частиц (7 сем) / prog5

#include <iostream>
#include <math.h>
#include <time.h>
#include "hist.hpp"

#ifndef M_PI
double M_PI = 3.14159265;
#endif

using namespace std;

// This function generates basic random
double gamma(){
        return (rand()+0.5)/(RAND_MAX+1.0);
}
// simulate exp(-3x) + 1/3/sqrt(2*M_PI)*exp((x-2)^2/8)
int main(){
        double ksi;  // random value
        srand(time(NULL)); // initialize the random generator
        int num_bin = 20;
        Hist *hist = new Hist(-1, 7, num_bin);
        for( int i=0; i<4000000; i++ ){
// *************************************************************

// *************************************************************
            hist->fill(ksi); // fill histogram
        }
        for(int i=0; i<hist->get_num_bin(); i++){
                cout << hist->get_bin(i) << "\t";   // get middle of bin boundary
                cout << hist->get_bin_value(i) << "\n"; // get number of counts in the bin
        }
        hist->draw(60);
        delete hist;  // destroy histogram

}