//SINGULAR Example 2.8.8




ring R   = 0,(x,y,z,u),dp;
ideal I  = x2y2-xyz2;
qring A  = std(I);
module B = [x,xz,y2],[y,yz,z2],[x2,xy,y2];
ideal ann= quotient(B,freemodule(3));  
ideal fit= minor(B,3);                 
ann;                            //annihilator of Coker(B)

fit;                            //Fitting ideal of Coker(B)

size(reduce(fit,std(ann)));
size(reduce(ann,std(fit)));
size(reduce(ann*ann,std(fit)));