LIB "tst.lib";
tst_init();
tst_ignore("CVS ID $Id$");
LIB "ncfactor.lib";

ring R = 0,(x1,d1,x2,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,1,0,0,
                  -1,0,0,0,
                  0,0,0,1,
                  0,0,-1,0;
 def r = nc_algebra(C,D);
 setring(r);

poly h = (x1*d2 + 5*x1 + 1)*(x2*d1 - 3*d2 +4);
list lsng = facWeyl(h);
print(lsng);
testNCfac(lsng, h);

tst_status();
tst_status(1); $