蒙特卡洛求圆周率PI

直接上代码：

#include <iostream>
#include <ctime>
using namespace std;

void main()
{
	srand(unsigned(time(NULL)));
	int total = 10000000;
	float k, j;
	int num = 0;
	for (int i = 0; i < total; i++)
	{
		k = 1.0*rand() / RAND_MAX;
		j = 1.0*rand() / RAND_MAX;
		if (k*k + j*j <= 1)
			num++;
	}
	cout << "Pi = " <<4.0* num / total << endl;
	system("pause");
}

蒙特卡洛探讨待续：