ring r = 0,(x,y,z),dp
module MD =
[1,x2+y,x*y-y*z,0,0],[0,x*z-y,z-x,0,0],[0,0,0,x2+y,x*y-y*z],[1,0,0,x*z-y,z-x];
module MK = syz(MD);
matrix MM = MK;
print(MM);
module N =[x2+y,x*z-y],[x*y-y*z,z-x];
quotient(N,freemodule(nrows(MK)));
from sympy import QQ
from sympy.abc import x
QQ.old_poly
I'm CCing Tom Bachmann, as he's the one who wrote this code, and is
probably the only one who can answer here.
Aaron Meurer
2014-01-31 Gerhard Martin :
> ring r = 0,(x,y,z),dp
>
> module MD =
> [1,x2+y,x*y-y*z,0,0],[0,x*z-y,z-x,0,0],[0,0,0,x2+y,x*y-y*z],[1,0,0,x*z-y,z-x];
>
> module MK = syz(MD);