mprimdec.lib
LIB "mprimdec.lib"; ring r=0,(x,y,z),lp; ideal i=x+1,y-1,z; i=std(i); ideal primId=primTest(i,z); primId; ==> primId[1]=z ==> primId[2]=y-1 ==> primId[3]=x+1 i=x,z2,yz,y2; i=std(i); primId=primTest(i); primId; ==> primId[1]=x ==> primId[2]=y ==> primId[3]=z