As a final note: some functions are defined for Domains, others for Polyhedra. A domain is simply a list of polyhedrons. Every polyhedron structure has a "next" pointer used to make a list of polyhedrons.
For instance, the union of two disjoint polyhedra is a domain consisting of two polyhedra. If you want the enclosing convex domain, you have to call DomainConvex(Pol1, 200) explicity. Note that the inclusion function does not work on domains, only on simple polyhedrons...
The output produced by the program is shown in Fig. .