  A self-consistent mean-field theory for bosons for T>0 is used to reconcile
predictions of collapse with recent observations of Bose-Einstein condensation
of 7Li. Eigenfunctions of a (non-separable) Hamiltonian that includes the
anisotropic external trap field and atom-atom interactions are obtained by an
iteration process. A sum over the Bose distribution, and the ``alternating
direction implicit'' algorithm are used. Near Tc, the ensemble exhibits a
localized condensate composed of atoms in the few lowest states. For lower T,
numerical instability indicates collapse to a more dense phase.
