dmod.lib
LIB "dmod.lib"; ring R = 0,(x,y),dp; poly F = x*y*(x+y); def A = annfsBM(F,0); setring A; LD; ==> LD[1]=x*Dx+y*Dy+3 ==> LD[2]=x*y*Dy+y^2*Dy+x+2*y ==> LD[3]=y^2*Dx*Dy-y^2*Dy^2+2*y*Dx-4*y*Dy-2 isHolonomic(LD); ==> 1 ideal I = std(LD[1]); I; ==> I[1]=x*Dx+y*Dy+3 isHolonomic(I); ==> 0