grid on
gh = get(gcf);
axes(gh.Children(3));
axis([10 100000 1e-2 1e4])
r1 = rectangle('position',[10 1e-2 145-1 2e3-1e-2]);
col = [0.85 0.85 0.85];
r1.FaceColor = col;
r1.EdgeColor = col;
uistack(r1,'down')
r2 = rectangle('position',[12566 1 1e5-12566 1e4-1]);
col = [0.85 0.85 0.85];
r2.FaceColor = col;
r2.EdgeColor = col;
uistack(r2,'down')

