3) Approximately 120/1000(Depending on the randomization)

var('x')
x=0
for y in range (1000):
    F=matrix([[randint(0,1) for i in range(1,11)] for j in range(1,11)])
    if (F.echelon_form()==identity_matrix(10)):
        x=x+1
show(x/1000)

for y in range (1000):
    F=matrix([[randint(0,1) for i in range(1,11)] for j in range(1,11)])
    F.echelon_form()==identity_matrix(10)

I need help with this formula! We are supposed to find the number of outputted Trues in this statement. Either a percent answer or just the number true out of 1000 would be awesome!

I'm not very well versed in this programming language and any help with a description would be great!!