The goal of this article is to develop a robust and accurate numerical method for solving hyperbolic conservation laws in three dimensions. The basic equations are the three-dimensional Euler equations describing the motion of an inviscid gas. The mathematical description of the model is represented by the system of equations of continuity, motion and energy (three dimensional nonstationary partial differential equations).