Under the influence of gravity, light scalar fields can form bound compact objects called boson stars. We use the technique of matching asymptotic expansions to obtain the profile for boson stars where the constituent particles have self-interactions. We obtain parametric representations of these profiles as a function of the self-interactions, including the case of very strong self-interactions. We show that our methods agree with solutions obtained by purely numerical methods. Significant distortions are found as compared to the noninteracting case.

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.