I am writing my own (orbital-free, for now) 3D density functional theory code from scratch. It currently works correctly for the simple harmonic oscillator potential $V_\text{ext}(\mathbf r) = \frac12\omega^2 |\mathbf r|^2$, yielding a ground state energy of $\frac32\omega$ (atomic units). I am t...