Text this: Nonlinear forced response of non-uniform beams carrying point masses using the harmonic balance method with an arc-length continuation scheme