interiorPoint -- computes a point in the relative interior of the Polyhedron

• Usage:

  p = interiorPoint P

• Inputs:
  • P, a convex polyhedron
• Outputs:
  • p, a matrix, over QQ with only one column representing a point

i1 : P = cyclicPolytope(3,5)

o1 = {ambient dimension => 3           }
      dimension of lineality space => 0
      dimension of polyhedron => 3
      number of facets => 6
      number of rays => 0
      number of vertices => 5

o1 : Polyhedron

i2 : interiorPoint P

o2 = | 2  |
     | 6  |
     | 20 |

              3       1
o2 : Matrix QQ  <-- QQ