function g = lorentzian(delta,t2);
% Computes the super-lorentzian integral given:
% delta -- the offset freq in Hz
% t2 -- T2 of bound pool protons in sec 

theta = 0 : pi/2000 :pi/2;

g = pi/2000*sqrt(2/pi)*sum(t2./abs(3.*cos(theta).^2-1).*exp(-2*((2*pi*delta*t2./(3*cos(theta).^2-1)).^2)).*sin(theta));
