We present a technique which permits the calculation of two-point functions of operators containing one heavy quark and an arbitrary number of other fixed mass quarks as analytic functions of the heavy-quark mass. It is based on the standard Jacobi linear solver used for the calculation of quark propagators. Masses and dispersion relations for heavy-light mesons are obtained on 16 3 × 48 lattices at β = 6.0 using an O( α)-improved “clover” fermion action and the standard Wilson action.