  We study primary submodules and primary decompositions from a differential
and computational point of view. Our main theoretical contribution is a general
structure theory and a representation theorem for primary submodules of an
arbitrary finitely generated module over a polynomial ring. We characterize
primary submodules in terms of differential operators and punctual Quot
schemes. Moreover, we introduce and implement an algorithm that computes a
minimal differential primary decomposition for a module.
