We present a method for tabulating all cubic function fields over $\mathbb{F}_{q}(t)$ whose discriminant $D$ has either odd degree or even degree and the leading coefficient of $-3D$ is a non-square in $\mathbb_{q}^*$, up to a given bound $X$ on $|D| = q^{\deg(D)}$. Our method is based on a generalization of Belabas' method for tabulating cubic number fields. The main theoretical ingredient is a generalization of a theorem of Davenport and Heilbronn to cubic function fields, along with a reduction theory for binary cubic forms that provides an efficient way to compute equivalence classes of binary cubic forms. The algorithm requires $O(q^4 X^{1+\epsilon})$ field operations when $D$ has odd degree, and $O(q^5 X^{1+\epsilon})$ when $D$ has even degree. It performs quite well in practice. The algorithm, examples and numerical data for $q=5,7,11,13$ are included.