from random import * def pgcd(a,b) : while b > 0 : r = a % b a = b b = r return a # question 1 : def simulation_cesaro(n): f = 0 for i in range(n) : a = randint(1, 100) b = randint(1, 100) if pgcd(a,b) == 1 : f = f + 1 f = f / n return f # question 2 : def proba_cesaro() : num = 0 for a in range(1, 101) : for b in range(1, 101) : if pgcd(a,b) == 1 : num = num + 1 p = num / (100**2) return p