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.

% Creates a shape from a number of curves. A shape is a cell array of curves.
function s = shape(varargin);
    s = varargin;
end