A probe-based harmonic balance method to simulate coupled oscillators