The homotopy analysis method and Galerkin spectral method are applied to find the analytical solutions for the Gross–Pitaevskii equations, a set of nonlinear Schr？dinger equation used in simulation of spin-1 Bose–Einstein condensates trapped in a harmonic potential. We investigate the one-dimensional case and get the approximate analytical solutions successfully. Comparisons between the analytical solutions and the numerical solutions have been made. The results indicate that they are in agreement well with each other when the atomic interaction is weakly. We also find a class of exact solutions for the stationary states of the spin-1 system with harmonic potential for a special case.
