  We present initial results from Monte Carlo simulations of NRQCD-light,
static-light, and NRQCD-NRQCD mesons, using a variational technique (MOST), as
part of our ongoing calculation of the $f_{B}$ decay constant. The basis states
for the variational calculation are quark-antiquark operators separated by all
possible relative distances not equivalent under the cubic group (for example,
for a $20^{3}$ lattice there are 286 operators). The efficacy of the method is
demonstrated by the good plateaus obtained for the ground state and the clean
extraction of the wave functions of the ground and first radially excited
state.
