  We construct the CP^n model on fuzzy sphere. The Bogomolny bound is saturated
by (anti-)self-dual solitons and the general solutions of BPS equation are
constructed. The dimension of moduli space describing the BPS solution on fuzzy
sphere is exactly the same as that of the commutative sphere or the
(noncommutative) plane. We show that in the soliton backgrounds, the number of
zero modes of Dirac operator on fuzzy sphere, Atiyah-Singer index, is exactly
given by the topological charge of the background solitons.
