    The hyperbolic conservation laws may develop discontinuities
in its solution even if the initial conditions are smooth. Therefore
classical numerical methods which depend on Taylor's expansion fails,
so the spurious oscillations occur in the solution. To resolve this,
the Total Variation Diminishing (TVD) schemes are constructed by Harten
in [1]}, [2]} based on the principle that, the total variation
to the approximation of numerical solution must be non-increasing
in time. But the TVD schemes are at most first order accurate near
smooth extrema [3]}.
