 LIB "dmodapp.lib";
ring r = 0,(x,y,Dx,Dy),dp;
def R = Weyl(); setring R; // Weyl algebra in variables x,y,Dx,Dy
poly F = x2y3;
ideal I = (y^3  x^2)*Dx  2*x, (y^3  x^2)*Dy + 3*y^2; // I = Dx*F, Dy*F;
// moreover I is not holonomic, since its dimension is not 2 = 4/2
gkdim(I); // 3
==> 3
def W = SDLoc(I,F); setring W; // creates ideal LD in W = R[s]
def U = DLoc0(LD, x2y3); setring U; // compute in R
LD0; // Groebner basis of the presentation of localization
==> LD0[1]=3*x*Dx+2*y*Dy+12
==> LD0[2]=3*y^2*Dx+2*x*Dy
==> LD0[3]=y^3*Dyx^2*Dy+6*y^2
BS; // description of bfunction for localization
==> [1]:
==> _[1]=0
==> _[2]=1/6
==> _[3]=1/6
==> [2]:
==> 1,1,1
