  Recently we have derived an effective one-dimensional nonpolynomial
Schr\"odinger equation (NPSE) [Phys. Rev. A {\bf 65}, 043614 (2002)] that
accurately describes atomic Bose-Einstein condensates under transverse harmonic
confinement. In this paper we analyze the stability of the solitary-wave
solutions of NPSE by means of the Vakhitov-Kolokolov criterion. Stable
self-focusing solutions of NPSE are precisely the 3D bright solitons
experimentally found in Bose-Einstein condensates [Nature {\bf 417}, 150
(2002); Science {\bf 296}, 1290 (2002)]. These self-focusing solutions
generalize the ones obtained with the standard nonlinear Schr\"odinger equation
(NLSE) and take into account the formation of instabilities in the transverse
direction. We prove that neither cubic nor cubic-quintic approximations of NPSE
are able to reproduce the unstable branch of the solitary-wave solutions.
Moreover, we show that the analitically found critical point of NPSE is in
remarkable good agreement with recent numerical computations of 3D GPE. We also
discuss the formation of multiple solitons by a sudden change of the sign of
the nonlinear strength. This formation can be interpreted in terms of the
modulational instability of the time-dependent macroscopic wave function of the
Bose condensate. In particular, we derive an accurate analytical formula for
the number of multiple bright solitons.
