Harmonic balance approach to the periodic solutions of the (an)harmonic relativistic oscillator

The first-order harmonic balance method via the first Fourier coefficient is used to construct two approximate frequency-amplitude relations for the relativistic oscillator for which the nonlinearity (anharmonicity) is a relativistic effect due to the time line dilation along the world line. Making a change of variable, a new nonlinear differential equation is obtained and two procedures are used to approximately solve this differential equation. In the first the differential equation is rewritten in a form that does not contain a square-root expression, while in the second the differential equation is solved directly. The approximate frequency obtained using the second procedure is more accurate than the frequency obtained with the first due to the fact that, in the second procedure, application of the harmonic balance method produces an infinite set of harmonics, while in the first procedure only two harmonics are produced. Both approximate frequencies are valid for the complete range of oscillation amplitudes, and excellent agreement of the approximate frequencies with the exact one are demonstrated and discussed. The discrepancy between the first-order approximate frequency obtained by means of the second procedure and the exact frequency never exceeds 1.6%. We also obtained the approximate frequency by applying the second-order harmonic balance method and in this case the relative error is as low 0.31% for all the range of values of amplitude of oscillation A.