demo.cpp

void demo(){
    double r = 0.5 ;
    double area = M_PI * r * r;
    double** samples = new double*[sample_count];
    for(int i=0;i<sample_count;i++){
        samples[i] = new double[2];
        samples[i][0] = 1.0*rand()/INT_MAX;
        samples[i][1] = 1.0*rand()/INT_MAX;
    }
    cout << "The exact area is " << area << endl ;
    mc(r, samples);
    verify(r, samples);
}