R = QQ[x,y]
I = ideal(x^3-1,y-x)
cv = isCellular (I,ReturnCellVars=>true)
pd = cellularBinomialPrimaryDecomposition (I,CellVariables=>cv)
mingens \ pd