   echo = 1
? mathnf([0,2])

[2]

? mathnf([0,2],1)
[Mat(2), [1, 0; 0, 1]]
? mathnf([0,x])

[x]

? mathnf([0,x],1)
[Mat(x), [1, 0; 0, 1]]
? mathnf([x;1])

[x]

[1]

? mathnf([x;1],1)
[[x; 1], Mat(1)]
? mathnf([x,x^2+1;x^3+x+1,x+2]*Mod(1,5))

[Mod(1, 5)*x^5 + Mod(2, 5)*x^3 + Mod(4, 5)*x + Mod(1, 5) Mod(4, 5)*x^4 + Mod
(2, 5)*x^3 + Mod(4, 5)*x^2 + Mod(3, 5)*x]

[0 1]

? v=[116085838,181081878,314252913,10346840];
? [H,U]=mathnf(v,1);[v*U,norml2(U)]
[[0, 0, 0, 1], 2833319]
? [H,U]=mathnf(v,5);[v*U,norml2(U)]
[[0, 0, 0, 1], 765585180708864230567243002686057927228240493]
? M=[0,0,0,0,0,0,0,0,0,13;0,0,0,0,0,0,0,0,23,6;0,0,0,0,0,0,0,23,-4,-7;0,0,0,0,0,0,17,-3,5,-5;0,0,0,0,0,56,16,-16,-15,-17;0,0,0,0,57,24,-16,-25,2,-21;0,0,0,114,9,56,51,-52,25,-55;0,0,113,-31,-11,24,0,28,34,-16;0,50,3,2,16,-6,-2,7,-19,-21;118,43,51,23,37,-52,18,38,51,28];mathnf(M)

[787850171872400 32189386376004 356588299060422 742392731867995 282253457851
430 665185047494955 664535243562463 744564809133574 113975061998590 52745901
3372200]

[0 12 6 11 5 3 7 6 6 0]

[0 0 3 1 2 1 1 0 0 0]

[0 0 0 1 0 0 0 0 0 0]

[0 0 0 0 1 0 0 0 0 0]

[0 0 0 0 0 1 0 0 0 0]

[0 0 0 0 0 0 1 0 0 0]

[0 0 0 0 0 0 0 1 0 0]

[0 0 0 0 0 0 0 0 1 0]

[0 0 0 0 0 0 0 0 0 1]

? mathnf([])
[;]
? mathnf([],1)
[[;], [;]]
? mathnf([;])
[;]
? mathnf([;],1)
[[;], [;]]
? mathnfmodid(matrix(0,2),[])
[;]
? mathnfmodid([0,7;-1,0;-1,-1],[6,2,2])

[2 1 1]

[0 1 0]

[0 0 1]

? mathnfmodid([;],[2])

[2]

? matsolvemod([;],[]~,[]~,1)
[[]~, [;]]
? matsolvemod([;],[],[]~,1)
  ***   at top-level: matsolvemod([;],[],[]~,1)
  ***                 ^-------------------------
  *** matsolvemod: incorrect type in matsolvemod (D) (t_VEC).
? matsolvemod([;],[]~,[],1)
  ***   at top-level: matsolvemod([;],[]~,[],1)
  ***                 ^-------------------------
  *** matsolvemod: incorrect type in matsolvemod (Y) (t_VEC).
? matsolvemod([;],2,[1]~,1)
0
? matsolvemod([;],[2]~,[1]~,1)
0
? matsolvemod([;],[2]~,1,1)
0
? matsolvemod([;],1,1,1)
[[]~, [;]]
? matsolvemod([1,2;3,4],1,2,1)
[[0, 0]~, [1, 0; 0, 1]]
? matsolvemod([;],[2]~,[2,3]~,1)
  ***   at top-level: matsolvemod([;],[2]~,[2,3]~,1)
  ***                 ^------------------------------
  *** matsolvemod: inconsistent dimensions in matsolvemod [2].
? print("Total time spent: ",gettime);
Total time spent: 1
