function y = lorentzian(x,params)
	x0 = params(1);
	delta = params(2);
	yS = params(3);
	yInf = params(4);
	y = yInf+(yS-yInf)*x0^2./(x0^2+1i*delta*x-x.^2);
end