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