Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-pro t purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. We present an efficient computational method for calculating the binding energies of the bound states of ultracold alkali-metal dimers in the presence of magnetic fields. The method is based on propagation of coupled differential equations and does not use a basis set for the interatomic distance coordinate. It is much more efficient than the previous method based on a radial basis set and allows many more spin channels to be included. This is particularly important in the vicinity of avoided crossings between bound states. We characterize a number of different avoided crossings in Cs 2 and compare our converged calculations with experimental results. Small but significant discrepancies are observed in both crossing strengths and level positions, especially for levels with l symmetry ͑rotational angular momentum L =8͒. The discrepancies should allow the development of improved potential models in the future.