  Using the theory of analytic functions of several complex variables, we prove
that if an analytic function in several variables satisfies a system of
$q$-partial differential equations, then, it can be expanded in terms of the
product of the homogeneous Hahn polynomials. Some non-trivial applications of
this expansion theorem to $q$-series are discussed.
