normaliz.lib
LIB "normaliz.lib"; ring R=0,(x,y,z),dp; intmat M[3][2]=3,1, 3,2, 1,3; normaliz(M,1); ==> 1,1, ==> 2,1, ==> 1,2, ==> 3,1, ==> 1,3