Abstract
Using the rational splitting of a cosine operator-function, the fourth order accuracy decomposition scheme is constructed for hyperbolic equation when the principal operator is self-adjoint positively defined and is represented as a sum of two summands. Stability of the constructed scheme is shown and the error of an approximate solution is estimated.
