function [x,y] = lorentz(mu,gamma)

% lorentz(mu,gamma) models the Resonance or Lorentz Profile
% see, Bevington and Robinson, 2nd Edition, Section 2.4
% mu is the resonance frequency
% gamma is the line width

x=[mu-5*gamma:.1:mu+5*gamma];
y=(gamma/(2*pi))./(((x-mu).^2)+(gamma/2)^2);
plot(x,y);
