if i == 0:
        return b * fact1 + logsf  # a