import math
import random


def euclidean_distance(x, y):
	return math.sqrt(x*x + y*y)

def montecarlo_simulation_for_calulcating_pii(num_iterations=100):

	desired_count = 0

	for i in range(0, num_iterations):
		x_coordinate = random.random()
		y_coordinate = random.random()

		if euclidean_distance(x_coordinate, y_coordinate) <= 1:
			desired_count += 1


	return float(desired_count) * 4/ float(num_iterations)


if __name__ == '__main__':
	print(montecarlo_simulation_for_calulcating_pii(10000))