#Monte Carlo method
import random

def estimate_pi(num_samples):
    inside_circle = 0

    for _ in range(num_samples):
        x = random.random()
        y = random.random()
        distance = x ** 2 + y ** 2
        if distance <= 1:
            inside_circle += 1

    return (inside_circle / num_samples) * 4

# Number of random samples for estimation
num_samples = 1000000
estimated_pi = estimate_pi(num_samples)
print(f"Estimated value of pi after {num_samples} samples: {estimated_pi}")