maximize                                       
  logexp_sum(matrix_allscenarios)              
Constraint: <= 4
  cardn(0.000001, matrix_cardinality)
Solver: precision = 9     
                                               
Value:                                         
  logistic(matrix_allscenarios)