function map = hot()

% generate hot colormap
%  map = hot();

cycle(1:24) = 0;
cycle(25:48) = 1/24:1/24:1; 
cycle(49:96) = 1;

r(1:64) = cycle(25:24+64);
g(1:64) = cycle(1:64);
b(1:48) = 0;
b(49:64) = 1/16:1/16:1;
map = [r' g' b'];
