  We provide a self-contained introduction to Gr\"obner bases of submodules of
$R[x_1, \ldots, x_n]^k$, where $R$ is a Euclidean domain, and explain how to
use these bases to solve linear systems over $R[x_1, \ldots, x_n]$.
