The Wikipedia article (http://en.wikipedia.org/wiki/Mathieu_function) is quite extensive, I believe, and in general there is NO way to explicitly calculate either eigenvalues or eigenfunctions of the Mathieu equation. 

The Wikipedia article also explicitly states that the solutions generally cannot be expressed in terms of hypergeometric functions. 

You can do perturbation theory around $q=\infty$ and recursively calculate asymptotic expansions of eigenfunctions and eigenvalues to any order (semiclassical analysis, basically). However, you seem not to be interested in this. (A reference would be a quite old Paper, unfortunately in german: J. Meixner, "Asymptotische Entwicklung der Eigenwerte und Eigenfunktionen der Differentialgleichungen der Sphäroid-Funktionen und der Mathieuschen Funktionen", Zeitschrift füt angewandte Mathematik und Mechanik, 1948)

-------------
Edit:

What I didn't mention if you are interested in approximations: Apart from the semiclassical expansion, the operator family
$$H(q) = - \frac{d^2}{d x^2} - q \cos(x)$$
is a self-adjoint family of unbounded operators, where each operator has the same domain (the Sobolev space $H^2(S^1)$ on the circle). By general perturbation theory for these operators (the standard reference is Kato, "Perturbation Theory for Linear Operators"), the eigenvalues are holomorphic in some domain including the real axis, and the eigenspaces are finite-dimensional as well and parametrized by the parameter $q$ in a holomorphic way.

I believe that you can (theoretically) calculate all coefficients of the corresponding Taylor series at $q=0$, both of eigenvalues and eigenfunctions.

As opposed to the other case above, this expansion will actually converge in a neighborhood near zero (the semiclassical expansion is only asymptotic). If I recall correctly, the eigenvalues $a_n(q)$ and $b_n(q)$ are holomorphic on some Riemann surface covering $\mathbb{C}$ which sadly has a very complicated branch structure (there are no real branch points, however, as no splitting of eigenvalues occurs on the real axis in the self-adjoint case). The problem for your purpose is that these branch points will limit the radius of convergence of the Taylor series at zero. If you have explicit values of $q$ in mind, it may or may not be enough. The locus of these branchpoints should probably be known, but I do not have a reference.