 LIB "matrix.lib";
ring r=(0,a,b),(A,B,C),dp;
matrix m[8][6]=
0, 2*C, 0, 0, 0, 0,
0, 4*C,a*A, 0, 0, 0,
b*B, A, 0, 0, 0, 0,
A, B, 0, 0, 0, 0,
4*C, 0, B, 2, 0, 0,
2*A, B, 0, 0, 0, 0,
0, 3*B, 0, 0, 2b, 0,
0, AB, 0, 2*A,A, 2a;"";
==>
print(colred(m));"";
==> 0,0,0,0, 2*C, 0,
==> 0,0,0,0, 4*C,(a)*A,
==> 0,0,0,(b)*B,A, 0,
==> 0,0,0,A, B, 0,
==> 1,0,0,0, 0, 0,
==> 0,0,0,2*A, B, 0,
==> 0,1,0,0, 0, 0,
==> 0,0,1,0, 0, 0
==>
list L=colred(m,1);
print(L[1]);
==> 0,0,0,0, 2*C, 0,
==> 0,0,0,0, 4*C,(a)*A,
==> 0,0,0,(b)*B,A, 0,
==> 0,0,0,A, B, 0,
==> 1,0,0,0, 0, 0,
==> 0,0,0,2*A, B, 0,
==> 0,1,0,0, 0, 0,
==> 0,0,1,0, 0, 0
print(L[2]);
==> 0, 0, 0, 1, 0, 0,
==> 0, 0, 0, 0, 1, 0,
==> 0, 0, 0, 0, 0, 1,
==> 1/2, 0, 0, 2*C, 0, 1/2*B,
==> 0, 1/(2b), 0, 0, 3/(2b)*B, 0,
==> 1/(2a)*A,1/(4ab)*A,1/(2a),2/(a)*AC,(2b+3)/(4ab)*AB,1/(2a)*AB
