//SINGULAR Example 2.8.7
ring R = 0,(x,y,z),(c,dp);
ideal I = x2y2-xyz2;
qring A = std(I);
poly f = xy*(y-z)*(y-1);
module N = [x,xz,y2],[y,yz,z2],[x2,xy,y2],[xy,xz,yz];
ideal ann= quotient(N,freemodule(3)); //annihilator of
// Coker(N)
ring Rt = 0,(t,x,y,z),dp;
ideal I = imap(R,I);
ideal ann= imap(A,ann),I;
poly f = imap(A,f);
ideal J = ann,1-t*f;
eliminate(J,t);
setring A;
size(quotient(N,f));