minimize
  pr_pen(0,matrix_L3)
Box: >= point_lb, <= point_ub
Solver: stages=30