def service_rule (model,i):
return sum(model.arrival [o,d,veh] * model.lam[i,o,d,veh] for o,d,veh in model.q) <= model.y_pc[i]
model.service = Constraint(model.Nodes, rule = service_rule)
def minimum_charger_rule (model, i):
return model.y_pc[i]*0.98 <=sum(model.arrival [o,d,veh] * model.lam[i,o,d,veh] for o,d,veh in model.q)
model.min_charge = Constraint(model.Nodes, rule = minimum_charger_rule)