Singular

LIB "nctools.lib";
ring r = 0,(x,y,z),dp;
matrix C[3][3]=0,1,2,0,0,-1,0,0,0;
print(C);
==> 0,1,2, 
==> 0,0,-1,
==> 0,0,0  
matrix D[3][3]=0,1,2y,0,0,-2x+y+1;
print(D);
==> 0,1,2y,     
==> 0,0,-2x+y+1,
==> 0,0,0       
def S=nc_algebra(C,D);setring S; S;
==> // coefficients: QQ
==> // number of vars : 3
==> //        block   1 : ordering dp
==> //                  : names    x y z
==> //        block   2 : ordering C
==> // noncommutative relations:
==> //    yx=xy+1
==> //    zx=2xz+2y
==> //    zy=-yz-2x+y+1
def l=ncRelations(S);
print (l[1]);
==> 0,1,2, 
==> 0,0,-1,
==> 0,0,0  
print (l[2]);
==> 0,1,2y,     
==> 0,0,-2x+y+1,
==> 0,0,0