def __cdf(self, params, x, backscatter_model):
        mode = params[0]
        k = params[1]
        offset = params[2]
        if backscatter_model == RadarHistogram.BACKSCATTER_MODEL_GAUSSIAN:
            return 0.5 * (1 + scipy.special.erf((x - mode) / (k * math.sqrt(2))))
        theta = (mode - offset) / (k - 1)
        return scipy.special.gammainc(k, (offset - x) / theta)

    # find x where __cdf(params, offset, x) = percentile
dataset/ETHPy150Open nasa/CrisisMappingToolkit/cmt/radar/histogram.py/RadarHistogram.__cdf/swapped operands of `-`