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
A-M baby 2013-04-10 21:29:02 -07:00			`import numpy as np`

			`n = input('Number of samples: ')`
Just another update 2015-02-20 00:35:01 -05:00			`print np.sum(np.random.rand(n)2+np.random.rand(n)2<1)/float(n)*4`