“…The existence and uniqueness of the weak answers for this group of equations have been discussed in [6,11,38,44]. In recent decades, finding numerical solutions of Equation (1.1) using different ideas has attracted the attention of researchers such as relaxation scheme [18], adaptive multiresolution schemes [12,15,16], finite volume schemes [8,13,19,35], kinetic scheme [9], local discontinuous Galerkin finite element method [46], mixed finite element methods [7,37,39] and WENO schemes [3,25,33,40]. In recent decades, numerical solution of partial differential equations has been considered by many researchers [1,2,4,5,29,30,34,43].…”