sage days 11 -- linear algebra

1293 days ago by pub

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()