Add RK6 with residual viscosity update and reduce computational effort of spatial scheme - Add RK6 with residual updates - Change the D2 operator for upwind schemes to one less computationally expensive.

function benchmark(func, N)
    default_arg('N',100);

    tic
    profile on

    for i = 1:N
        func();
    end

    profile viewer
end