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 prof(f)
    profile on
    try
        f();
        profile viewer
    catch e
        fprintf(2, '\n%s', getReport(e));
        profile clear
    end
end