for {alpha_op,step}={0.795000, start}, {0.785000, 0.795000}, {0.775000, 0.785000}, {0.765000, 0.775000}, {0.755000, 0.765000}, {0.745000, 0.755000}, {0.735415, 0.745000}

Problem: problem_alpha_op, type = minimize
  10*pm_pen_g(0,matrix_new_1)
  pm_pen(0,matrix_new_0)
Constraint: <= alpha_op 
  pm_pen(-1,L(matrix_0)-L(matrix_1))
Solver: init_point = point_problem_step

end for