d=[.001:.001:.009,.01:.01:.09,.1:.1:.9,1:9,10:10:90,100:100:900,1000:1000:10000];
nd=numel(d);
f1=.0001;
f2=1000000;
hold on
for i =1:nd
pointd(1,:)=[f1,d(i)*2*pi*f1];
pointd(2,:)=[f2,d(i)*2*pi*f2];
xd=logspace(log10(pointd(1,2)),log10(pointd(2,2)));
yd=logspace(log10(f2),log10(f1));
plot(xd,yd,'color','black')
end
g=[.001:.001:.009,.01:.01:.09,.1:.1:.9,1:9,10:10:90,100:100:900,1000:1000:10000];
ng=numel(g);
for i =1:nd
pointg(1,:)=[f1,g(i)*386/(2*pi*f1)];
pointg(2,:)=[f2,g(i)*386/(2*pi*f2)];
xg=logspace(log10(pointg(1,2)),log10(pointg(2,2)));
yg=logspace(log10(f2),log10(f1));
plot(xg,yg,'color','black')
end
set(gca,'xscale','log')
set(gca,'yscale','log')
axis([1 10000 1 1000])
grid on
set(gca,'GridLineStyle','-')