function y = expln(x)
  if (x > 0)
     y = exp(-x) - log(x);
  else
     y = 0;
  end
end
