 LIB "ncModslimgb.lib";
ring r = 0,(x,y),dp;
poly P = y^4+x^3+x*y^3; // a (3,4)Reiffen curve
def A = Sannfs(P); setring A; // computed Dmodule data from P
ideal bs = LD, imap(r,P); // preparing the computation of the BernsteinSato polynomial
ideal I1 = ncmodslimgb(bs,0,2); // no final verification, use 2 cores
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
==> // ** going to redefine the algebra structure
I1[1]; // the BernsteinSato polynomial of P, univariate in s
==> s^7+7*s^6+499/24*s^5+815/24*s^4+227563/6912*s^3+43627/2304*s^2+4461779/74\
6496*s+595595/746496
ideal I2 = ncmodslimgb(bs); // do the final verification, use 1 core (default)
I2[1]; // the BernsteinSato polynomial of P, univariate in s
==> s^7+7*s^6+499/24*s^5+815/24*s^4+227563/6912*s^3+43627/2304*s^2+4461779/74\
6496*s+595595/746496
