minimize
mip = 1
  linear(matrix_production_and_inventory_cost)
  +fxchg_pos(0.001, matrix_set_up_costs)  
Constraint: == vector_constraint_multiple_demand
  linearmulti(matrix_multiple_demand)  
Box: >= 0, <= point_upperbounds
Solver: cargrb