Calculate the average distance between 2 random points inside a 1x1 square using a Monte Carlo simulation

random_points_in_square.py

from random import random  # the random function yelds values between 0 and 1


def distance(x1,y1,x2,y2):
    return ((x1-x2)**2+(y1-y2)**2)**(1/2)

def results(n):
    for i in range(n):
        yield distance(random(),random(),random(),random())

N = 10**6
print(sum(results(N))/N)  # result is a very unintuitive ~0.521. An analitical solution is presented at youtube.com/watch?v=i4VqXRRXi68