def eoq(I,F,h,d,w,W,a0,aK,K):
    """eoq --  multi-item capacitated economic ordering quantity model
    Parameters:
        - I: set of items
        - F[i]: ordering cost for item i
        - h[i]: holding cost for item i
        - d[i]: demand for item i
        - w[i]: unit weight for item i
        - W: capacity (limit on order quantity)
        - a0: lower bound on the cycle time (x axis)
        - aK: upper bound on the cycle time (x axis)
        - K: number of linear pieces to use in the approximation
    Returns a model, ready to be solved.
    """

    # construct points for piecewise-linear relation, store in a,b
    a,b = {},{}
    delta = float(aK-a0)/K
    for i in I:
        for k in range(K):
            T = a0 + delta*k
            a[i,k] = T                          # abscissa: cycle time
            b[i,k] = F[i]/T + h[i]*d[i]*T/2.    # ordinate: (convex) cost for this cycle time

    model = Model("multi-item, capacitated EOQ")

    x,c,w_ = {},{},{}
    for i in I:
        x[i] = model.addVar(vtype="C", name="x(%s)"%i)  # cycle time for item i
        c[i] = model.addVar(vtype="C", name="c(%s)"%i)  # total cost for item i
        for k in range(K):
            w_[i,k] = model.addVar(ub=1, vtype="C", name="w(%s,%s)"%(i,k)) #todo ??

    for i in I:
        model.addCons(quicksum(w_[i,k] for k in range(K)) == 1)
        model.addCons(quicksum(a[i,k]*w_[i,k] for k in range(K)) == x[i])
        model.addCons(quicksum(b[i,k]*w_[i,k] for k in range(K)) == c[i])

    model.addCons(quicksum(w[i]*d[i]*x[i] for i in I) <= W)

    model.setObjective(quicksum(c[i] for i in I), "minimize")

    model.data = x,w
    return model
def eoq ( I , F , h , d , w , W , a0 , aK , K ) : a , b = { } , { } delta = float ( aK - a0 ) / K for i in I : for k in range ( K ) : T = a0 + delta * k a [ i , k ] = T b [ i , k ] = F [ i ] / T + h [ i ] * d [ i ] * T / 2. model = Model ( "multi-item, capacitated EOQ" ) x , c , w_ = { } , { } , { } for i in I : x [ i ] = model . addVar ( vtype = "C" , name = "x(%s)" % i ) c [ i ] = model . addVar ( vtype = "C" , name = "c(%s)" % i ) for k in range ( K ) : w_ [ i , k ] = model . addVar ( ub = 1 , vtype = "C" , name = "w(%s,%s)" % ( i , k ) ) for i in I : model . addCons ( quicksum ( w_ [ i , k ] for k in range ( K ) ) == 1 ) model . addCons ( quicksum ( a [ i , k ] * w_ [ i , k ] for k in range ( K ) ) == x [ i ] ) model . addCons ( quicksum ( b [ i , k ] * w_ [ i , k ] for k in range ( K ) ) == c [ i ] ) model . addCons ( quicksum ( w [ i ] * d [ i ] * x [ i ] for i in I ) <= W ) model . setObjective ( quicksum ( c [ i ] for i in I ) , "minimize" ) model . data = x , w return model
eoq --  multi-item capacitated economic ordering quantity model
    Parameters:
        - I: set of items
        - F[i]: ordering cost for item i
        - h[i]: holding cost for item i
        - d[i]: demand for item i
        - w[i]: unit weight for item i
        - W: capacity (limit on order quantity)
        - a0: lower bound on the cycle time (x axis)
        - aK: upper bound on the cycle time (x axis)
        - K: number of linear pieces to use in the approximation
    Returns a model, ready to be solved.