pi_monte_carlo.py

import random
total=1000000
in_circle=0
for i in range(total):
  (x,y)= (random.random(), random.random())
  if x**2 + y**2 <= 1:
    in_circle+=1
  
4*(in_circle/total)