Contragenic functions are defined to be reduced-quaternion-valued harmonic functions which are orthogonal to all monogenic and antimonogenic functions in the \(L^2\) norm of a given domain. The parallelism between the spaces of contragenic functions in the interior and exterior of the unit sphere in \(\mathbb {R}^3\) is described in detail. Bergman reproducing kernels for the spaces of contragenic functions are given, mirroring the corresponding kernels for the spaces of vector parts of monogenic functions. Numerical examples are given showing the accuracy of truncations of the integral kernels. A striking duality is observed between the basic interior contragenic functions and the vector parts of exterior monogenic functions, and vice versa.
