Harmonic functions in three-dimensional Euclidean space conceal a subtle and unexpected structure. Unlike their counterparts in the complex plane or in quaternionic space, they cannot, in general, be decomposed even locally into the sum of a monogenic and an antimonogenic function. This failure gives rise to “contragenics”—harmonic functions lying in the \(L^2\) -orthogonal complement of the space of all such sums. This survey traces the development of the theory of contragenics in stages. It begins with the spheroidal harmonics, both internal and external, on spheroidal domains of arbitrary eccentricity. Building on these, we construct graded one-parameter orthogonal bases of spheroidal monogenics and ambigenics, emphasizing their algebraic structure and limiting behavior that reflect the geometry of the underlying spheroid. Contragenics then emerge, first within the interior of spheroidal domains, and later in both the interior and exterior of the unit sphere, where they furnish a richer \(L^2\) -orthogonal decomposition of harmonic functions in quaternionic analysis on three-dimensional, simply connected, unbounded Euclidean domains. Finally, we extend the framework to the exterior of spheroidal domains, where logarithmic terms intervene and further enrich the theory.
