In this paper we present the numerical solutions of neutral atoms with attractive interaction in a harmonic trap. The calculated result shows the bistability curve of the number of condensate atoms versus the energy eigenvalues. The maximum number of atoms in and attractive Bose-Einstein condensate deduced hereby is in agreement with the experimental on the whole.