  The anharmonic Bloch oscillations of a light beam in the array of optical
waveguides are considered. The coupling modes model (CMM) with the second order
interaction is used to describe the effect analytically. The formula obtained
predicts explicitly the path of the optical beam, in particular, the positions
of the turning points are found. A total agreement of this formula with the
numerical simulation is confirmed.
