for i in range(10000):
    # пытаюсь проверить 10000 случайных матриц 2Х2
    zero = [] # сюда записываю результаты проверки определителей
    i =  np.random.randint(0,10,[2,2])
    if int(np.linalg.det(i)) == 0: # сравниваю определитель одной случайной матрицы с 0
        zero.append(1) # если определитель случайной матрицы равен 0, то в переменную zero добавляю 1
    else:
        zero.append(0) # если определитель случайной матрицы отличен от 0, то в переменную zero добавляю 0
# пытаюсь оценить вероятность появления линейно-зависимых матриц среди 10000 случайных
p = sum(zero)/10000

  
>>> import numpy as np
>>> 
>>> zero = []
>>> for _ in range(10000):
...     det = np.random.randint(0, 10, [2, 2])
...     if int(round(np.linalg.det(det))) == 0:
...         zero.append(1)
... 
>>> p = sum(zero) / 10000
>>> print(p)
0.0557
>>>