Задача 2 - игральный кубик

Пусть у нас есть кубик, который с равной долей вероятности выпадает на числа 1,2,3,4,5,6. Такой кубик в Python выполняет функция randint из модуля random; можно записывать случайно сгенерированное значение в переменную:

import random
a = random.randint(1, 6)  

На примере такого “кубика” посмотрим на закон больших чисел. Так, значение вероятности $P(\{1,2\})$ выпадения чисел 1 или 2 в теории равно $\frac{1}{3},$ а на практике стремится к $\frac{1}{3}$ с увеличением числа бросков. На $N$-ом броске будем считать “практической” вероятностью частное кол-ва нужных нам событий (т.е. выпадения 1 или 2) на общее число бросков кубика $N.$

Для числа $\epsilon >0$ (назовём погрешностью) найдите номер броска $N^{\ast },$ для которого отклонение практической и теоретической вероятностей меньше $\epsilon :$ $|P(\{1,2\})-\frac{1}{3}|<\epsilon .$

Далее, для следующего массива погрешностей $x=[{\epsilon }_1,...{\epsilon }_M]$ найдите массив необходимого числа бросков $y=[N_1^{\ast },...N_M^{\ast }],$ изобразите на графике:

import matplotlib.pyplot as plt
 
x = [1e-1, 1e-2, 1e-3, 1e-4, 1e-5]
 
y = ... # ваш код
 
plt.plot(x,y)
plt.show()