Singular

LIB "ncfactor.lib";
ring R = 0,(x,s),dp;
def r = nc_algebra(1,s);
setring(r);
poly h = (s^2*x+x)*s;
facFirstShift(h);
==> [1]:
==>    [1]:
==>       1
==>    [2]:
==>       s
==>    [3]:
==>       s2+1
==>    [4]:
==>       x-1
==> [2]:
==>    [1]:
==>       1
==>    [2]:
==>       s2+1
==>    [3]:
==>       s
==>    [4]:
==>       x-1
==> [3]:
==>    [1]:
==>       1
==>    [2]:
==>       s2+1
==>    [3]:
==>       x
==>    [4]:
==>       s
==> [4]:
==>    [1]:
==>       1
==>    [2]:
==>       s
==>    [3]:
==>       xs2+s2+x-1