Hi Jeanne,
As per my understanding, you want to shade the area representing the standard deviation with a different colour.
You can consider the following example to fill colour between two different plots:
x = 1 : 20;
% curve1 and curve2 representing the data and its deviation
curve1 = sin(x);
curve2 = sin(x) + 0.5;
plot(x, curve1, 'r', 'LineWidth', 2);
hold on;
plot(x, curve2, 'b', 'LineWidth', 2);
x2 = [x, fliplr(x)];
inBetween = [curve1, fliplr(curve2)];
fill(x2, inBetween, [1 0.7 1]);
For further information on the functions used in the code, kindly refer to the following links:
