//SINGULAR Procedure4.6.4 proc prepareSat(ideal i) { int k; poly p=leadcoef(i[1]); for(k=2;k<=size(i);k++) { p=p*leadcoef(i[k]); } return(p); }