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