Following the general theory developed in paper I [ Chin. Phys. Lett. 26 (2009) 120506 ] we present here the numerical results for an N spinless boson system in a 1D harmonic trap with repulsive delta function interparticle interaction as N → ∞, with √ N / g = fixed.