dmod.lib
LIB "dmod.lib"; ring X = 0,(x,y,z,t),dp; poly q = arrange(3); factorize(q,1); ==> _[1]=x ==> _[2]=x+z ==> _[3]=y ==> _[4]=y+z ==> _[5]=x+y ==> _[6]=x+y+z ==> _[7]=z