#!/usr/bin/python3 import random stats = [0]*10 N = 100000 i = 0 while i < N: a = random.randrange(1,100000) b = random.randrange(1,100000) c = random.randrange(1,100000) result = a*b*c firstdigit = int(str(result)[0]) stats[firstdigit] += 1 i += 1 for d, n in enumerate(stats): print('{0}% of results start with digit {1}'.format(round(n/N*100),d))