function [fval,grad] = Lorentzian(y,delta)

absy = abs(y);
fval = sum(log(1 + absy/delta^2));
grad = 1./(absy + delta^2); % grad at |y|
