Add support for artificial viscosity to the 1d burgers scheme. - Add support for artificial viscosity by parametrizing the discretization operator and boundary closure on the viscosity. - Add a time stepper which evaluates and updates the residual viscosity of the solution.

function r = rowVector(v)
    if  size(v,1) == 1
        r = v
    elseif size(v,2) == 1
        r = v';
    else
        error('v is not a matrix');
    end
end