Аннотация:A new approach is considered for increasing the order of accuracy of the gridcharacteristic method in the region of coefficient jumps. The approach is based on piecewise
polynomial interpolation for schemes of the second and third orders of accuracy for the case
where the interface between the media is consistent with a finite-difference grid. The method is
intended for numerical simulation of the propagation of dynamic wave disturbances in heterogeneous media. Systems of hyperbolic equations with variable coefficients are used to describe
the considered physical processes. The description of the numerical method and the results of
its testing are given.