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','-')