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