minpoly([-5,0,-1,1;0,-5,-1,-1;1,1,-5,0;-1,1,0,-5])
minpoly([-5,0,-1,1;0,-5,-1,-1;1,1,-5,0;-1,1,0,-5]/3)
minpoly([-5,0,-1,1;0,-5,-1,-1;1,1,-5,0;-1,1,0,-5]/a)
minpoly([-5,0,-1,1;0,-5,-1,-1;1,1,-5,0;-1,1,0,-5]*a)
matsnf([1,2;3,4],2)
matsnf([-X-5,-1,-1,0;0,X^2+10*X+26,-1,-X-5;1,-X-5,-X-5,1;-1,0,0,1],2)
minpoly([a+4/51*d,29/51*d,55/17*d,-9/17*d;-11/51*d,a+35/51*d,-11/17*d,12/17*d;23/51*d,1/51*d,a+6/17*d,-5/17*d;11/17*d,16/17*d,33/17*d,a-19/17*d])
{
  for(i=1,15,
      print(matfrobenius(matrix(i,i,k,j,k>=j),2)));
  for(i=1,10,
     print(matfrobenius(matpascal(i),2)));
  for(i=1,8,
     print(matfrobenius(mathilbert(i),2)));
}
lindep([1,2,3], -3)
matsnf(matdiagonal([x,0,x]), 2)
matsnf(matdiagonal([1+x,0,x]), 2)
\\Bug #1208
matsnf([0;1;2], 4)
matsnf([0;1;2], 5)
matsnf(Mat([0,1,2]), 4)
matsnf(Mat([0,1,2]), 5)
