  We numerically construct the BPS and non-BPS wavefunctions of an
$\mathcal{N}=4$ quiver quantum mechanics with two Abelian nodes and a single
arrow. This model captures the dynamics of a pair of wrapped D-branes
interacting via a single light string mode. A dimensionless parameter $\nu$,
which is inversely proportional to the Fayet-Iliopoulos parameter, controls
whether the bulk of the wavefunctions are supported on the Higgs branch or the
Coulomb branch. We demonstrate how the BPS and excited states morph as $\nu$ is
tuned. We also numerically compute the energy gap between the ground state and
the first excited states as a function of $\nu$. An expression for the gap,
computed on the Coulomb branch, matches nicely with our numerics at large $\nu$
but deviates at small $\nu$ where the Higgs branch becomes the relevant
description of the physics. In the appendix, we provide the Schr\"{o}dinger
equations fully reduced via symmetries which, in principle, allow for the
numerical determination of the entire spectrum at any point in moduli space.
For the ground states, this numerical determination of the spectrum can be
thought of as the first \emph{in silico} check of various Witten index
calculations.
