Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics