For bosons in a harmonic oscillator-like trap we will use is a spherical (S) or an elliptical (E) harmonic trap in one, two and finally three dimensions, with the latter given by
Vext(r)={12mω2hor2(S)12m[ω2ho(x2+y2)+ω2zz2](E)where (S) stands for symmetric and
ˆH=N∑i(−ℏ22m▽2i+Vext(ri))+N∑i<jVint(ri,rj),as the two-body Hamiltonian of the system.