| LIB "mondromy.lib";
ring R=0,(x,y),ds;
poly f=x2y2+x6+y6;
matrix M=monodromyB(f);
print(M);
==> 7/6,0, 0,0, 0, 0,0, 0,-1/2,0, 0, 0, 0,
==> 0, 7/6,0,0, 0, 0,-1/2,0,0, 0, 0, 0, 0,
==> 0, 0, 1,0, 0, 0,0, 0,0, 0, 0, 0, 0,
==> 0, 0, 0,4/3,0, 0,0, 0,0, 0, 0, 0, 0,
==> 0, 0, 0,0, 4/3,0,0, 0,0, 0, 0, 0, 0,
==> 0, 0, 0,0, 0, 1,0, 0,0, 0, 0, 0, 0,
==> 0, 0, 0,0, 0, 0,5/6, 0,0, 0, 0, 0, 0,
==> 0, 0, 0,0, 0, 0,0, 1,0, 0, 0, 0, 0,
==> 0, 0, 0,0, 0, 0,0, 0,5/6, 0, 0, 0, 0,
==> 0, 0, 0,0, 0, 0,0, 0,0, 2/3,0, 0, 0,
==> 0, 0, 0,0, 0, 0,0, 0,0, 0, 2/3,0, 0,
==> 0, 0, 0,0, 0, 0,0, 0,0, 0, 0, 1, -1/3,
==> 0, 0, 0,0, 0, 0,0, 0,0, 0, 0, 3/4,0
|