We define U(n) gauge theory on fuzzy S2 N × S2 N as a multi-matrix model, which reduces to ordinary Yang-Mills theory on S2×S2 in the commutative limit N → ∞. The model can be used as a regularization of gauge theory on noncommutative Rθ in a particular scaling limit, which is studied in detail. We also find topologically non-trivial U(1) solutions, which reduce to the known “fluxon” solutions ...