Allow for the residual normalization function to return a vector. Also change the default normalization from normalizing on norm(u-mean(u),inf) to norm(u/2), since this proved to be better for 1d burgers

function saveFigurePosition()
    defaultPosition = get(0,'defaultfigureposition');
    f = gcf;
    defaultPosition(3:4) = f.Position(3:4);
    set(0,'defaultfigureposition',defaultPosition);
end