Towards an Analytic Construction of the Wavefunction of Boson Stars
  Light scalar fields can form gravitationally bound compact objects called
boson stars. The profile of boson stars in the Newtonian limit is described by
the Gross-Pitaevskii-Poisson equations. We present a semi-analytic solution to
these equations and construct the profile of boson stars formed by a
non-interacting scalar field. Our solution is stable with respect to numerical
errors and has accuracy better than $10^{-6}$ over the entire range.
