RosettaCodeData/Task/Monte-Carlo-methods/Python/monte-carlo-methods-3.py

import numpy as np

n = input('Number of samples: ')
print np.sum(np.random.rand(n)**2+np.random.rand(n)**2<1)/float(n)*4