 LIB "purityfiltration.lib";
ring D = 0,(x,y,z),dp;
matrix R[6][4]=
0,2*x,z2*yx,1,
0,z2*x,2*y3*x,1,
z,6*x,2*y5*x,1,
0,yx,yx,0,
y,x,yx,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
