Re: Plotting feasible region of a linear programme

*To*: mathgroup at smc.vnet.net*Subject*: [mg112822] Re: Plotting feasible region of a linear programme*From*: Yaroslav Bulatov <yaroslavvb at gmail.com>*Date*: Sat, 2 Oct 2010 05:43:55 -0400 (EDT)*References*: <i81j79$a07$1@smc.vnet.net>

On Sep 30, 1:54 am, Yaroslav Bulatov <yarosla... at gmail.com> wrote: > What is a good way to visualize feasible region of a 3-variable linear > programme in Mathematica? I tried RegionPlot3D but you can see the > result is lacking in aesthetics > > randomCons := Module[{}, > hadamard = KroneckerProduct @@ Table[{{1, 1}, {1, -1}}, {3}]; > invHad = Inverse[hadamard]; > vs = Range[8]; > m = mm /@ vs; > sectionAnchors = Subsets[vs, {1, 7}]; > randomSection := > Mean[hadamard[[#]] & /@ #] & /@ > Prepend[RandomChoice[sectionAnchors, 3], vs]; {p0, p1, p2, p3} = > randomSection; > section = > Thread[m -> > p0 + {x, y, z}.Orthogonalize[{p1 - p0, p2 - p0, p3 - p0}]]; > And @@ Thread[invHad.m >= 0 /. section] > ]; > Table[RegionPlot3D @@ {randomCons, {x, -3, 3}, {y, -3, 3}, {z, -3, > 3}}, {10}] I got a suggestion from Maxim Rytin to use a legacy package InequalityPlot, and it works much better than RegionPlot3D for this problem. I've put comparison in parallel stackoverflow thread -- http://stackoverflow.com/questions/3815496/plotting-linear-inequalities-in-mathematica. Two issues remain -- 1) it creates more polygons than necessary, so can't use EdgeForm to outline edges of the region 2) sometimes faces are missed