function D=contrast(C,s)
% contrast(C,s)
% adjust contrast of colourmap - moves all 
% colours towards a medium grey.
% C - colourmap (N by 3 matrix)
% s - contrast: 1 no change, 0 completely grey

D=s*C + (1-s)*0.5*ones(size(C));




