 LIB "dmodapp.lib";
ring r = 0,(x,y,Dx,Dy),dp;
def R = Weyl(); // Weyl algebra on the variables x,y,Dx,Dy
setring R;
poly F = x2y3;
ideal I = Dx*F, Dy*F;
// note, that I is not holonomic, since it's dimension is not 2
gkdim(I); // 3, while dim R = 4
==> 3
def W = SDLoc(I,F);
setring W; // = R[s], where s is a new variable
LD; // Groebner basis of sparametric presentation
==> LD[1]=3*x*Dx*s+2*y*Dy*s6*s^2+6*s
==> LD[2]=3*y^2*Dx*s+2*x*Dy*s
==> LD[3]=y^3*Dyx^2*Dy3*y^2*s+3*y^2
==> LD[4]=y^3*Dxx^2*Dx+2*x*s2*x
