function [y] = LorentzFunc(v, x)
background = v(1);
amplitude = v(2);
center = v(3);
width = v(4);

y = background - abs(amplitude).*width.^2./((x-center).^2+width.^2);

end