 LIB "arr.lib";
ring R = 0,(x,y,z),dp;
arr A= ideal(x,y,xy+1); // centerless
arrCenter(A);
==> [1]:
==> 0
arr B= ideal(x,y,z); // center is a single point
arrCenter(B);
==> [1]:
==> 1
==> [2]:
==> _[1,1]=0
==> _[2,1]=0
==> _[3,1]=0
==> [3]:
==> _[1,1]=0
arr C= ideal(x,z,x+z); // center is a line
// here we get a wrong result because the matrix is simplified since A doesn't
// contain any "y" the matrix (Ab) will be 3x3 only.
arrCenter(C);
==> [1]:
==> 1
==> [2]:
==> _[1,1]=0
==> _[2,1]=0
==> _[3,1]=0
==> [3]:
==> _[1,1]=0
==> _[2,1]=1
==> _[3,1]=0
