input: Write a python function to generate points for Monte Carlo simulation within the unit circle.
  assert len(monte_carlo_simulation(50)) == 50
  assert all(-1 <= x <= 1 and -1 <= y <= 1 for x, y in monte_carlo_simulation(100))
  assert isinstance(monte_carlo_simulation(10), list)
output: import random
  
  def monte_carlo_simulation(n):
      points = []
      for _ in range(n):
          while True:
              x = random.uniform(-1, 1)
              y = random.uniform(-1, 1)
              if x**2 + y**2 <= 1:
                  points.append((x, y))
                  break
      return points