text: We will consider the fundamentals of the MUSCL scheme by considering the following simple first-order, scalar, 1D system, which is assumed to have a wave propagating in the positive direction,
  Where  represents a state variable and  represents a flux variable.
  The basic scheme of Godunov uses piecewise constant approximations for each cell, and results in a first-order upwind discretisation of the above problem with cell centres indexed as .  A semi-discrete scheme can be defined as follows,
  This basic scheme is not able to handle shocks or sharp discontinuities as they tend to become smeared. An example of this effect is shown in the diagram opposite, which illustrates a 1D advective equation with a step wave propagating to the right. The simulation was carried out with a mesh of 200 cells and used a 4th order Runge–Kutta time integrator (RK4).
  To provide higher resolution of discontinuities, Godunovs scheme can be extended to use piecewise linear approximations of each cell, which results in a central difference scheme that is second-order' accurate in space. The piecewise linear approximations are obtained from
  Thus, evaluating fluxes at the cell edges we get the following semi-discrete scheme
  where  and  are the piecewise approximate values of cell edge variables, i.e.,
  Although the above second-order scheme provides greater accuracy for smooth solutions, it is not a total variation diminishing (TVD) scheme and introduces spurious oscillations into the solution where discontinuities or shocks are present. An example of this effect is shown in the diagram opposite, which illustrates a 1D advective equation , with a step wave propagating to the right. This loss of accuracy is to be expected due to Godunov's theorem. The simulation was carried out with a mesh of 200 cells and used RK4 for time integration.
  MUSCL based numerical schemes extend the idea of using a linear piecewise approximation to each cell by using slope limited left and right extrapolated states. This results in the following high resolution, TVD discretisation scheme,
  Which, alternatively, can be written in the more succinct form,
  The numerical fluxes  correspond to a nonlinear combination of first and second-order approximations to the continuous flux function.
  The symbols  and  represent scheme dependent functions (of the limited extrapolated cell edge variables), i.e.,
  where, using downwind slopes:
  and
  The function  is a limiter function that limits the slope of the piecewise approximations to ensure the solution is TVD, thereby avoiding the spurious oscillations that would otherwise occur around discontinuities or shocks - see Flux limiter section. The limiter is equal to zero when  and is equal to unity when . Thus, the accuracy of a TVD discretization degrades to first order at local extrema, but tends to second order over smooth parts of the domain.
  The algorithm is straight forward to implement. Once a suitable scheme for  has been chosen, such as the Kurganov and Tadmor scheme (see below), the solution can proceed using standard numerical integration techniques.