f=.3:.1:3.4;
for n=1:length(f)
   c(n)=pi*f(n)*125*10^(-3)/sin(pi*f(n)*125*10^(-3));
   d(n)=20*log10(c(n));
end
plot(f,d)