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)