f = 0:0.01:0.5;
plot([0.25 0.25],[1 0], '--', ...
[0.125 0.125],[1 0], ...
f,sin(pi*f) ./ (pi*f)), ...
label('',18,'f','SH',18)