In this paper, we attempt to display effective numerical simulations of shock wave propagating represented by the Burgers equations known as a significant mathematical model for turbulence. A high order hybrid approximation based on seventh order weighted essentially non-oscillatory finite difference together with the sixth order finite difference scheme implemented for spatial discretization is presented and applied without any transformation or linearization to the Burgers equation and its modified form. Then, the produced system of first order ordinary differential equations is solved by the MacCormack method. The efficiency, accuracy and applicability of the proposed technique are analyzed by considering three test problems for several values of viscosity that can be caused by the steep shock behavior. The performance of the method is measured by some error norms. The results are in good agreement with the results reported previously, and moreover, the suggested approximation relatively comes to the forefront in terms of its low cost and easy implementation.