minimize
  meanabs_pen(matrix_regmatrixscenarios)
Constraint: = 1
  linear(matrix_sum)
Box: >= 0, <= 1
Solver: car, precision = 9