We present a new model for Yang-Mills theory on the fuzzy sphere in which the configuration space of gauge fields is given by a coadjoint orbit. In the classical limit it reduces to ordinary Yang-Mills theory on the sphere. We find all classical solutions of the gauge theory and use nonabelian localization techniques to write the partition function entirely as a sum over local contributions from critical points of the action, which are evaluated explicitly. The partition function of ordinary Yang-Mills theory on the sphere is recovered in the classical limit as a sum over instantons. We also apply abelian localization techniques and the geometry of symmetric spaces to derive an explicit combinatorial expression for the partition function, and compare the two approaches. These extend the standard techniques for solving gauge theory on the sphere to the fuzzy case in a rigorous framework. © 2007 Springer-Verlag.