53.sagews

M = matrix(2, 2, [1,2,pi,e])

M

show(M^(-1))

show(expand(M.charpoly()))

M = random_matrix(ZZ,200)
time M.det(proof=False)

x = random_matrix(ZZ,200,1)
time w = M \ x

show(w[0])

A = random_matrix(GF(2),2048,density=0.001)
time B = A.echelon_form()

B.visualize_structure(maxsize=2048)

# Linbox -- Clement Pernet

# M4RI -- Pernet, Albrecht

# IML --

# Lots of code by W. Stein, D. Harvey, R. Bradshaw, and M. Albrecht

# Client code

M = ModularSymbols(200,sign=1)

M.decomposition()

set_verbose(1)
ModularSymbols(37,sign=1).decomposition()