close all
clear
clc

% Colormap name     <-------------------------------------- PARAMETER HERE
cmapname = 'CmapOctaveGenerated.xml';

% Vector of values (equispaced on HSV map)     <----------- PARAMETER HERE
Yval = [2e10, 4e10, 6e10, 8e10, 1e11, 2e11, 4e11, 6e11];

% ==============   MAIN   ==============
nval = numel(Yval);
Hval = linspace(0.66667,0,nval); % from blue (H = 0.66667) to red (H = 0)

% HSV values
HSVval = [Hval', ones(nval,1), ones(nval,1)];
RGBval = hsv2rgb(HSVval);

% Plot colormap just for fun
figure
hold on
for ii = 1:nval
  plot(ii,Yval(ii),'o','color',RGBval(ii,:),'linewidth',3)
end

% Writing colormap
fid = fopen(cmapname,'w');

fprintf(fid, '<ColorMap name="%s" space="HSV">\n', cmapname(1:end-4)); % Not the extension!

for ii = 1:nval-1
  Ynow1 = Yval(ii);
  Ynow2 = Ynow1*(1+0.00001);

  fprintf(fid, '  <Point x="%f" o="1" r="%f" g="%f" b="%f"/>\n', Ynow1, RGBval(ii,1), RGBval(ii,2), RGBval(ii,3));
  fprintf(fid, '  <Point x="%f" o="1" r="%f" g="%f" b="%f"/>\n', Ynow2, RGBval(ii+1,1), RGBval(ii+1,2), RGBval(ii+1,3));
end

fprintf(fid, '  <Point x="%f" o="1" r="%f" g="%f" b="%f"/>\n', Yval(end), RGBval(end,1), RGBval(end,2), RGBval(end,3));

fwrite(fid, '</ColorMap>');
fclose(fid);