I don't have a code, but the plot was probably made in Matlab with the function 'shadederrorbar'.

Thank you Stefan for your answer, that's right I will try it.

Hello!

I use function patch in MATLAB

for example:

dd = [cDp; cDp(end:-1:1)] / 1e3;

clf

patch(err2(jj), dd(jj), 'y', 'FaceColor', 0.7 * ones(1, 3), 'EdgeColor', 0.7 * ones(1, 3))

hold on

plot(nanmedian(aa, 2), cDp / 1e3, 'rs-', 'LineWidth', 2, 'MarkerFaceColor', [1, 0, 0])

Andrey, thank you very much for your code sharing.