  In the present work, we improve a numerical method, developed to solve the
Gross-Pitaevkii nonlinear Schroedinger equation. A particular scaling is used
in the equation, which permits to evaluate the wave-function normalization
after the numerical solution. We have a two point boundary value problem, where
the second point is taken at infinity. The differential equation is solved
using the shooting method and Runge-Kutta integration method, requiring that
the asymptotic constants, for the function and its derivative, are equal for
large distances. In order to obtain fast convergence, the secant method is
used.
