how to get the smith normal form of lambda matrix
Right now, I know sage matrix can get the smith_form and elementary_divisors on pure integer matrix. But I want to get the same stuff on lambda matrix. What I actually mean is:
http://www.solitaryroad.com/c153.html (www.solitaryroad.com/c153.html) (my karma is not enough to put a direct link)
I directly construct a matrix and the minus lambda times identity matrix and call smith_form is not working. So how to get the smith normal form of lambda matrix in sagemath? Or if sagemath can not do that, which computer algebra system can do the job?