| LIB "purityfiltration.lib";
ring D = 0,(x,y,z),dp;
matrix R[6][4]=
0,-2*x,z-2*y-x,-1,
0,z-2*x,2*y-3*x,1,
z,-6*x,-2*y-5*x,-1,
0,y-x,y-x,0,
y,-x,-y-x,0,
x,-x,-2*x,0;
is_pure(transpose(R));
==> 0
// R is not pure, so we do the purity filtration
list T=purityFiltration(transpose(R));
// all Elements of T[2] are either zero or pure
purelist(T[2]);
==> [1]:
==> 1
==> [2]:
==> 1
==> [3]:
==> 1
==> [4]:
==> 1
|