function l = log_(a)

l = log(a + (a == 0)*0.0001);

end