Hello Onur,

Your can use the Climate Data Toolbox (https://github.com/chadagreene/CDT) with its funcionalities described here: Article The Climate Data Toolbox for MATLAB

The shadow range represents the uncertainty between models, and you can use quartiles to define upper and lower limits, for example:

Slow=quantile(ysample,0.125); Shigh=quantile(ysample,0.875);

h1=patch([yr' yr(end:-1:1)'],[Shigh Slow(end:-1:1)],[0.8 0.8 0.8],'linestyle','none');

Ilow=quantile(ys,0.125); Ihigh=quantile(ys,0.875);

h3=patch([yr' yr(end:-1:1)'],[Ihigh Ilow(end:-1:1)],[1 0.7 0.7],'linestyle','none','facealpha',0.8);

Climate Data Toolbox in Matlab can do this.