minimize
  pr_pen(0,matrix_h)
Constraint: ==  0.421725668503024
  linear(matrix_orthogonal)
Solver: stages=30, init_point  = point_initial