The Gross-Pitaevskii equation admits soliton solutions such as: $$\psi(x)=\psi_0 sech(x/\xi),$$ where $\xi$ is the healing length defined by: $\xi=\frac{\hbar}{\sqrt{m \mu}}$, with $\mu$ being the chemical potential. To obtain this solution, one of the boundary conditions is given by assuming that for lengths much greater than the healing length, the condensate is homogeneus with density $\psi_0=\sqrt{\frac{\mu}{g}}$.
I don't see how to obtain this solution by solving the non-linear GPE: $$\mu \psi(x) = -\frac{\hbar^2}{2m} \partial^2_x \psi(x) + g|\psi(x)|^2 \psi(x).$$