minimize
mip=1
  pr_pen(0, matrix_H)
Constraint: = 1
  linearmulti(matrix_Aeq)
Box: >= 0
Solver: CARGRB,init_point=0.0001,timelimit=100