| LIB "ncfactor.lib";
ring R = 0,(x1,x2,d1,d2),dp;
matrix C[4][4] = 1,1,1,1,
1,1,1,1,
1,1,1,1,
1,1,1,1;
matrix D[4][4] = 0,0,1,0,
0,0,0,1,
-1,0,0,0,
0,-1,0,0;
def r = nc_algebra(C,D);
setring(r);
poly h = (d1+1)^2*(d1 + x1*d2);
facWeyl(h);
==> [1]:
==> [1]:
==> 1
==> [2]:
==> d1+1
==> [3]:
==> d1+1
==> [4]:
==> x1*d2+d1
==> [2]:
==> [1]:
==> 1
==> [2]:
==> x1*d1*d2+d1^2+x1*d2+d1+2*d2
==> [3]:
==> d1+1
|