function grad = logpen(x,epsilon)

absx = abs(x);
grad = 1./(epsilon + absx);