def zipf2min(params, x, data):
    """ model Zipf-Mandelbrot power law f = C/(r+b)^a """
    c = params['c']
    b = params['b']
    a = params['a']
    model = c/(x + b)**a
    # print(c, b, a, model)
    return model - data