Create residual viscosity package +rv and generalize the ResidualViscosity class - Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update. - Create separate RungekuttaRV specifically using interior RV updates - Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class

function [update_data, line_handle] = updatableLine()
    figure_handle = gcf;
    axis_handle = gca;

    line_handle = line(0,0);

    function update(x,y)
        if ishandle(figure_handle) && ishandle(axis_handle)
            set(line_handle,'XData',x);
            set(line_handle,'YData',y);
            drawnow
        end
    end
    update_data = @update;
end