First-order approximation of the number-projected (NP) SO(2N) Tamm-Dancoff (TD) equation is developed to describe ground and excited states of superconducting fermion systems. We start from an NP Hartree-Bogoliubov (HB) wave function. The NP SO(2N) TD expansion is generated by quasi-particle pair excitations from the degenerate geminals in the number-projected HB wave function. The Schrödinger equation is cast into the NP SO(2N) TD equation by the variation principle. We approximate it up to first order. This approximate equation is reduced to a simpler form by the Schur function of group characters which has a close connection with the soliton theory on the group manifold.