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 (b).  Find the maximum of   [Graphics:Images/LinearProgrammingMod_gr_293.gif]   over the region [Graphics:Images/LinearProgrammingMod_gr_294.gif].  

Solution 3 (b).

Graph the domain defined by the constraints.

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


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

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

 

 

Enter the objective function.

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


[Graphics:../Images/LinearProgrammingMod_gr_334.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_335.gif]
[Graphics:../Images/LinearProgrammingMod_gr_336.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_337.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_338.gif] for determining the pivot rows.

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


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

 

 

Run the subroutine  SimplexMethod.

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


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

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

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

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

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

 

 

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

From column 1 and row 2 we see that [Graphics:../Images/LinearProgrammingMod_gr_347.gif].  From column 2 and row 1 we see that [Graphics:../Images/LinearProgrammingMod_gr_348.gif].  The final point  [Graphics:../Images/LinearProgrammingMod_gr_349.gif]  is the optimal feasible solution.  The simplex method has moved from the origin along the edges to the point [Graphics:../Images/LinearProgrammingMod_gr_350.gif].  The value of the objective function at the solution point is located in the bottom row of the column "[Graphics:../Images/LinearProgrammingMod_gr_351.gif]",  i.e.  [Graphics:../Images/LinearProgrammingMod_gr_352.gif].

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

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

 

 

We are done!

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

    
[Graphics:../Images/LinearProgrammingMod_gr_355.gif]   for rows   [Graphics:../Images/LinearProgrammingMod_gr_356.gif] in the tableau.    

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

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


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

 

 

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


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

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

 

 

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

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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005