Calculate π using Monte Carlo

mc.cpp

#include <iostream>

using namespace std;

int main()
{
    double x, y;
    double N = 100000;
    int c = 0;

    for(int i = 0; i < N; i++)
    {
        x = double(random()) / RAND_MAX;
        y = double(random()) / RAND_MAX;

        if(x * x + y * y <= 1.0)
            c++;
    }

    cout << 4.0 * c / N << endl;

    return 0;

}