def METHOD_NAME(X):
    t = [X[i] for i in range(nbdates)]
    a = X[22]
    c = X[23]
    b = math.exp(c)
    y = [0.0] * nbdates
    for i in range(nbdates):
        y[i] = a * y0 / (b * y0 + (a - b * y0) * math.exp(-a * (t[i] - t0)))
    z = [yi / populationFactor for yi in y]  # Convert into millions
    return z