Procedure from library mprimdec.lib (see mprimdec_lib).
Usage:
modDec (N[, i]); module N, int i
Return:
list l
a minimal primary decomposition of N
computed by an generalized version of the algorithm of Shimoyama/Yokoyama,
if i=1 is given, the factorizing Groebner basis algorithm is used internally.