x=xlsread('data.xlsx',1,'B2:B10000');    %load data
[mu,sigma] = normfit(x);     %Looking for data mean and standard deviation
b1=1.1;   %left edge of color area
b2=2.2;    %right edge of color area
z=x-b1;       
[~,f1]=min(abs(z));  %Find the coordinates of the left boundary
z=x-b2;       
[~,f2]=min(abs(z));  %Find the coordinates of the right boundary

d=pdf('norm',x,mu,sigma);      %Get probability density
figure
plot(x,d,'.')                  %Draw probability density curve
hold on
fill([x(f1:f2);x(f2);x(f1)],[d(f1:f2);0;0],[0.3,0.4,1])   %Color filling