integrate(-log(1 - y)/y, (y, 0, 1))