Estimating pi by the monte carlo method

pi_monte_carlo.jl

using Distributions

function estimate_pi(n) # the higher n is, the more accurate the estimate will be
  hits = 0
  for i=1:n
    x = rand(Uniform(-1,1))
    y = rand(Uniform(-1,1))
    if x*x + y*y <= 1
      hits += 1
    end
  end
  4*(hits/n)
end