This paper presents an efficient numerical scheme for calculating the periodic motion of a harmonically forced piecewise linear oscillator very accurately. The scheme is based on the shooting technique with the traditional numerical Poincare mapping and its Jacobian replaced by the piecewise analytic ones. Thus, the scheme gets rid of the requirement of the current schemes for an assumed order of the oscillator trajectory passing through different linear regions. The numerical examples in the paper demonstrate that the new scheme, compared with the current schemes, enables one to cope with more complicated dynamics of harmonically forced piecewise linear oscillators.
