A = 5000 ;
k = 10000 ;
m = 1 ;
c = 10 ;
w = linspace( 10 , 1000 , 1000 ) ;
B = A ./ abs( k - m*w .^ 2 + 1i*w*c ) ;
p = plot( w , B , 'LineWidth' , 1.5 ) ;
assi = gca ;
assi.XScale = 'log' ;
assi.YScale = 'log' ;
assi.XLabel.String = 'w' ;
assi.YLabel.String = 'B(w)' ;
assi.TickDir = 'out' ;
assi.XGrid = true ;
assi.YGrid = true ;
assi.XMinorGrid = true ;
assi.YMinorGrid = true ;