1 | initial version |

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)

2 | No.2 Revision |

```
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)~~show(x/1000)

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.