% plot reference solution 
%   (must first be computed on finer grid in subdirectory qref)

plotqref