  Let $\zeta$ be a fixed nonzero element in a finite field $\mathbb F_q$ with
$q$ elements. In this article, we count the number of pairs $(A,B)$ of $n\times
n$ matrices over $\mathbb F_q$ satisfying $AB=\zeta BA$ by giving a generating
function. This generalizes a generating function of Feit and Fine that counts
pairs of commuting matrices. Our result can be also viewed as the point count
of the variety of modules over the quantum plane $xy=\zeta yx$, whose geometry
was described by Chen and Lu.
