Example 3.  Consider the region [Graphics:Images/LinearProgrammingMod_gr_289.gif] in the plane defined by the inequalities:  

        [Graphics:Images/LinearProgrammingMod_gr_290.gif]  

3 (a).  Find the maximum of   [Graphics:Images/LinearProgrammingMod_gr_291.gif]   over the region [Graphics:Images/LinearProgrammingMod_gr_292.gif].  

Solution 3 (a).

Graph the domain defined by the constraints.

[Graphics:../Images/LinearProgrammingMod_gr_295.gif]


[Graphics:../Images/LinearProgrammingMod_gr_296.gif]

[Graphics:../Images/LinearProgrammingMod_gr_297.gif]

 

 

Enter the objective function.

[Graphics:../Images/LinearProgrammingMod_gr_298.gif]


[Graphics:../Images/LinearProgrammingMod_gr_299.gif]

 

 

Use the Mathematica subroutine LinearProgramming to find the solution (note carefully the "-" signs in front of the matrix and vector in the second and third arguments).

 

[Graphics:../Images/LinearProgrammingMod_gr_300.gif]
[Graphics:../Images/LinearProgrammingMod_gr_301.gif]

Use the simplex algorithm to find this solution.  

Enter the coefficients of the decision variables and right side of the tableau.  

 

[Graphics:../Images/LinearProgrammingMod_gr_302.gif]

Mathematica will fill in the columns for the slack variables and append a column of zeros which will be use to calculate the ratios [Graphics:../Images/LinearProgrammingMod_gr_303.gif] for determining the pivot rows.

[Graphics:../Images/LinearProgrammingMod_gr_304.gif]


[Graphics:../Images/LinearProgrammingMod_gr_305.gif]

 

 

Run the subroutine  SimplexMethod.

[Graphics:../Images/LinearProgrammingMod_gr_306.gif]


[Graphics:../Images/LinearProgrammingMod_gr_307.gif]

[Graphics:../Images/LinearProgrammingMod_gr_308.gif]

[Graphics:../Images/LinearProgrammingMod_gr_309.gif]

[Graphics:../Images/LinearProgrammingMod_gr_310.gif]

[Graphics:../Images/LinearProgrammingMod_gr_311.gif]

 

 

There are no negative coefficients in the bottom row, so we are done.

From column 1 and row 6 we see that [Graphics:../Images/LinearProgrammingMod_gr_312.gif].  From column 2 and row 5 we see that [Graphics:../Images/LinearProgrammingMod_gr_313.gif].  The final point  [Graphics:../Images/LinearProgrammingMod_gr_314.gif]  is the optimal feasible solution.  The simplex method has moved from the origin along the edges to the point [Graphics:../Images/LinearProgrammingMod_gr_315.gif].  The value of the objective function at the solution point is located in the bottom row of the column "[Graphics:../Images/LinearProgrammingMod_gr_316.gif]",  i.e.  [Graphics:../Images/LinearProgrammingMod_gr_317.gif].

[Graphics:../Images/LinearProgrammingMod_gr_318.gif]

[Graphics:../Images/LinearProgrammingMod_gr_319.gif]

 

 

We are done!

Aside.  The slack variables could be computed from (4) using the initial coefficients of the tableau:

    
[Graphics:../Images/LinearProgrammingMod_gr_320.gif]   for rows   [Graphics:../Images/LinearProgrammingMod_gr_321.gif] in the tableau.    

The bottom line of the tableau corresponds to the augmented objective function.

[Graphics:../Images/LinearProgrammingMod_gr_322.gif]


[Graphics:../Images/LinearProgrammingMod_gr_323.gif]

 

 

[Graphics:../Images/LinearProgrammingMod_gr_324.gif]


[Graphics:../Images/LinearProgrammingMod_gr_325.gif]

[Graphics:../Images/LinearProgrammingMod_gr_326.gif]

 

 

Graph the domain defined by the constraints and the level curves of the objective function.

[Graphics:../Images/LinearProgrammingMod_gr_327.gif]


[Graphics:../Images/LinearProgrammingMod_gr_328.gif]

[Graphics:../Images/LinearProgrammingMod_gr_329.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005