if i == 1:
        return (a - logsf) / fact1  # b