Procedure from library mprimdec.lib (see mprimdec_lib).
mprimdec.lib
annil(N); module N
ideal ann=std(quotient(N,freemodule(nrows(N)))); the annihilator of M/N in the basering
ann is a std basis in the basering
Example:
LIB "mprimdec.lib"; ring r=0,(x,y,z),dp; module N=x*gen(1), y*gen(2); ideal ann=annil(N); ann; → ann[1]=xy