customColormap2.m

function customColormap = customColormap2(nColors)

% Define desired color values at specific points
colors = [ ...
    0.2 0.2 1;    % Light blue for negative values
    1 1 1;        % White around 0
    1 1 1;        % White between -0.2 and 0
    1 1 0.8;      % Yellow for low positive values
    1 0.8 0.4;    % Orange for mid-range positive values
    1 0 0;        % Red for high positive values
];

% Define values representing color transitions
values = [-0.8, -0.2, 0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6];

% Create empty colormap array
customColormap = zeros(nColors, 3);

% Define new values linearly spaced throughout the desired range
newValues = linspace(-0.8, 1.6, nColors);

% Interpolate each color component separately
for i = 1:3
    interpolatedValues = interp1(values, colors(:, i), newValues, 'linear');
    % Ensure RGB values are within [0, 1] range
    customColormap(:, i) = max(0, min(1, interpolatedValues));
end

end