monte-carlo-pi

monte-carlo-pi.py

from math import sqrt                                                                                                                                                             
import random                                                                   
                                                                                
points, inside, radius = 50_000, 0, 1                                           
                                                                                
for point in range(points):                                                     
    # Generate random x and y coordinate                                        
    rand_x = random.random()                                                    
    rand_y = random.random()                                                    
                                                                                
    # For the point (x,y) in a quadrant of a circle the                         
    # lengths x and y become the legs of a right triangle                       
    # whose hypotenuse is less than 1 if inside the circle.                     
    hypotenuse = sqrt(rand_x  ** 2 + rand_y ** 2)                               
    if hypotenuse < radius:                                                     
        inside += 1                                                             
                                                                                
approx = 4 * inside / points                                                    
print("Pi Approx:\t{}".format(approx))