function map = phasemap(n)
% PHASEMAP(N) circular colormap with N entries
% See also redblue, pixcmap, phasemap, clmap rscolmap, zebra, cmap5,
% makecmap
%
% John Wilkin

if nargin == 0
  n = size(colormap,1);
end

s = linspace(0,2*pi,n);
b = 0.5*(1+sin(s));
r = 0.5*(1+sin(s+pi/3));
g = 0.5*(1+sin(s+2*pi/3));
map = [r; g; b]';
