15ME81 or Module 2
15ME81 or Module 2
15ME81 or Module 2
Module 2
1.1 Introduction
General Linear Programming Problem (GLPP)
Subject to constraints
a11x1 + a12x2 + …..........+a1nxn (≤ or ≥) b1
a21x1 + a22x2 + ………..+a2nxn (≤ or ≥) b2
.
.
.
am1 x1 + am2x2 + ……….+amnxn (≤ or ≥) bm
and
x1 ≥ 0, x2 ≥ 0,…, xn ≥ 0
Where constraints may be in the form of any inequality (≤ or ≥) or even in the form of an
equation (=) and finally satisfy the non-negativity restrictions.
1
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Step 3 – The right side element of each constraint should be made non-negative
2x1 +x2 – s2 = -15
-2x1 - x2 + s2 = 15 (That is multiplying throughout by -1)
Step 4 – All variables must have non-negative values.
For example: x1 +x2 ≤ 3
x1 > 0, x2 is unrestricted in sign
Then x2 is written as x2 = x2 – ׳x2 ׳׳where x2 ׳, x20 ≤ ׳׳
Therefore the inequality takes the form of equation as x1 + (x2 – ׳x2 )׳׳+ s1 = 3
Using the above steps, we can write the GLPP in the form of SLPP.
Example 1
Maximize Z = 3x1 + x2
Subject to
2 x1 + x2 ≤ 2
3 x1 + 4 x2 ≥ 12
and x1 ≥ 0, x2 ≥ 0
SLPP
Maximize Z = 3x1 + x2
Subject to
2 x1 + x2 + s1 = 2
3 x1 + 4 x2 – s2 = 12
x1 ≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0
Example 2
Minimize Z = 4x1 + 2 x2
Subject to
3x1 + x2 ≥ 2
x1 + x2 ≥ 21
x1 + 2x2 ≥ 30
and x1 ≥ 0, x2 ≥ 0
SLPP
׳
Maximize Z4 – = x1 – 2 x2
Subject to
3x1 + x2 – s1 = 2
x1 + x2 – s2 = 21
x1 + 2x2 – s3 = 30
x1 ≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0, s3 ≥ 0
Example 3
Minimize Z = x1 + 2 x2 + 3x3
Subject to
2
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Basic solution
Is a solution obtained by setting any ‘n’ variable equal to zero and solving remaining ‘m’
variables. Such ‘m’ variables are called Basic variables and ‘n’ variables are called Non-basic
variables.
It was developed by G. Danztig in 1947. The simplex method provides an algorithm (a rule of
procedure usually involving repetitive application of a prescribed operation) which is based
on the fundamental theorem of linear programming.
The Simplex algorithm is an iterative procedure for solving LP problems in a finite number of
steps. It consists of
Having a trial basic feasible solution to constraint-equations
Testing whether it is an optimal solution
3
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Improving the first trial solution by a set of rules and repeating the process till an
optimal solution is obtained
Advantages
Simple to solve the problems
The solution of LPP of more than two variables can be obtained.
1.5 Computational Procedure of Simplex Method
Consider an example
Maximize Z = 3x1 + 2x2
Subject to
x1 + x2 ≤ 4
x1 – x2 ≤ 2
and x1 ≥ 0, x2 ≥ 0
Solution
Cj → 3 2 0 0
Basic CB XB X1 X2 S1 S2 Min ratio
Variables XB /Xk
s1 0 4 1 1 1 0
s2 0 2 1 -1 0 1
Z= CB XB Δj
Step 4 – Calculation of Z and Δj and test the basic feasible solution for optimality by the
rules given.
Z= CB XB
= 0 *4 + 0 * 2 = 0
4
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Δj = Zj – Cj
= CB X j – Cj
Δ1 = CB X1 – Cj = 0 * 1 + 0 * 1 – 3 = -3
Δ2 = CB X2 – Cj = 0 * 1 + 0 * -1 – 2 = -2
Δ3 = CB X3 – Cj = 0 * 1 + 0 * 0 – 0 = 0
Δ4 = CB X4 – Cj = 0 * 0 + 0 * 1 – 0 = 0
Procedure to test the basic feasible solution for optimality by the rules given
Rule 1 – If all Δj ≥ 0, the solution under the test will be optimal. Alternate optimal solution
will exist if any non-basic Δj is also zero.
Rule 2 – If atleast one Δj is negative, the solution is not optimal and then proceeds to improve
the solution in the next step.
Rule 3 – If corresponding to any negative Δj, all elements of the column Xj are negative or
zero, then the solution under test will be unbounded.
In this problem it is observed that Δ1 and Δ2 are negative. Hence proceed to improve this solution
Step 5 – To improve the basic feasible solution, the vector entering the basis matrix and
the vector to be removed from the basis matrix are determined.
Incoming vector
The incoming vector Xk is always selected corresponding to the most negative value of
Δj. It is indicated by (↑).
Outgoing vector
The outgoing vector is selected corresponding to the least positive value of minimum
ratio. It is indicated by (→).
Step 6 – Mark the key element or pivot element by ‘ ‘.The element at the intersection
of outgoing vector and incoming vector is the pivot element.
Cj → 3 2 0 0
Basic CB XB X1 X2 S1 S2 Min ratio
Variables (Xk) XB /Xk
s1 0 4 1 1 1 0 4/1=4
s2 0 2 1 -1 0 1 2 / 1 = 2 → outgoing
↑incoming
Z= CB XB = 0 Δ1= -3 Δ2= -2 Δ3=0 Δ4=0
If the number in the marked position is other than unity, divide all the elements of
that row by the key element.
Then subtract appropriate multiples of this new row from the remaining rows, so as
to obtain zeroes in the remaining position of the column Xk.
5
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Step 7 – Now repeat step 4 through step 6 until an optimal solution is obtained.
Example 1
Maximize Z = 80x1 + 55x2
Subject to
4x1 + 2x2 ≤ 40
2x1 + 4x2 ≤ 32
and x1 ≥ 0, x2 ≥ 0
Solution
SLPP
6
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Cj → 80 55 0 0
Basic CB XB X1 X2 S1 S2 Min ratio
Variables XB /Xk
s1 0 40 4 2 1 0 40 / 4 = 10→ outgoing
s2 0 32 2 4 0 1 32 / 2 = 16
↑incoming
Z= CB XB = 0 Δ1= -80 Δ2= -55 Δ3=0 Δ4=0
(R1=R1 / 4)
x1 80 10 1 1/2 1/4 0 10/1/2 = 20
(R2=R2– 2R1)
s2 0 12 0 3 -1/2 1 12/3 = 4→ outgoing
↑incoming
Z = 800 Δ1=0 Δ2= -15 Δ3=40 Δ4=0
(R1=R1– 1/2R2)
x1 80 8 1 0 1/3 -1/6
(R2=R2 / 3)
x2 55 4 0 1 -1/6 1/3
Example 2
Maximize Z = 5x1 + 3x2
Subject to
3x1 + 5x2 ≤ 15
5x1 + 2x2 ≤ 10
and x1 ≥ 0, x2 ≥ 0
Solution
7
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
SLPP
Maximize Z = 5x1 + 3x2 + 0s1 + 0s2
Subject to
3x1 + 5x2+ s1= 15
5x1 + 2x2 + s2= 10
x1 ≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0
Cj → 5 3 0 0
Basic CB XB X1 X2 S1 S2 Min ratio
Variables XB /Xk
s1 0 15 3 5 1 0 15 / 3 = 5
s2 0 10 5 2 0 1 10 / 5 = 2 → outgoing
↑incoming
Z= CB XB = 0 Δ1= -5 Δ2= -3 Δ3=0 Δ4=0
(R1=R1– 3R2)
s1 0 9 0 19/5 1 -3/5 9/19/5 = 45/19 →
(R2=R2 /5)
x1 5 2 1 2/5 0 1/5 2/2/5 = 5
↑
Z = 10 Δ1=0 Δ2= -1 Δ3=0 Δ4=1
(R1=R1 / 19/5)
x2 3 45/19 0 1 5/19 -3/19
Example 3
Maximize Z = 5x1 + 7x2
Subject to
x1 + x2 ≤ 4
3x1 – 8x2 ≤ 24
10x1 + 7x2 ≤ 35
8
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
and x1 ≥ 0, x2 ≥ 0
Solution
SLPP
Maximize Z = 5x1 + 7x2 + 0s1 + 0s2 + 0s3
Subject to
x1 + x2 + s1= 4
3x1 – 8x2 + s2= 24
10x1 + 7x2 + s3= 35
x1 ≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0, s3 ≥ 0
Cj → 5 7 0 0 0
Basic CB XB X1 X2 S1 S2 S3 Min ratio
Variables XB /Xk
s1 0 4 1 1 1 0 0 4 /1 = 4→outgoing
s2 0 24 3 -8 0 1 0 –
s3 0 35 10 7 0 0 1 35 / 7 = 5
↑incoming
Z= CB XB = 0 -5 -7 0 0 0 ←Δj
x2 7 4 1 1 1 0 0
(R2 = R2 + 8R1)
s2 0 56 11 0 8 1 0
(R3 = R3 – 7R1)
s3 0 7 3 0 -7 0 1
Z = 28 2 0 7 0 0 ←Δj
Example 4
Maximize Z = 2x – 3y + z
Subject to
3x + 6y + z ≤ 6
4x + 2y + z ≤ 4
x–y+z≤3
and x ≥ 0, y ≥ 0, z ≥ 0
Solution
9
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
SLPP
Maximize Z = 2x – 3y + z + 0s1 + 0s2 + 0s3
Subject to
3x + 6y + z + s1= 6
4x + 2y + z + s2= 4
x – y + z + s3= 3
x ≥ 0, y ≥ 0, z ≥ 0 s1 ≥ 0, s2 ≥ 0, s3 ≥ 0
Cj → 2 -3 1 0 0 0
Basic CB XB X Y Z S1 S2 S3 Min ratio
Variables XB /Xk
s1 0 6 3 6 1 1 0 0 6/3=2
s2 0 4 4 2 1 0 1 0 4 / 4 =1→ outgoing
s3 0 3 1 -1 1 0 0 1 3/1=3
↑incoming
Z=0 -2 3 -1 0 0 0 ←Δj
↑incoming
Z=2 0 4 1/2 0 1/2 0 ←Δj
8/3
10
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Example 5
11
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Solution
SLPP
Maximize Z = 3x1 + 5x2 + 0s1 + 0s2 + 0s3
Subject to
3x1 + 2x2 + s1= 18
x1 + s2= 4
x2 + s3= 6
x1 ≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0, s3 ≥ 0
Cj → 3 5 0 0 0
Z = 36 0 0 1 0 3 ←Δj
12
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Example 6
Minimize Z = x1 – 3x2 + 2x3
Subject to
3x1 – x2 + 3x3 ≤ 7
-2x1 + 4x2 ≤ 12
-4x1 + 3x2 + 8x3 ≤ 10
and x1 ≥ 0, x2 ≥ 0, x3 ≥ 0
Solution
SLPP
׳
Min (-Z) = Max Z = -x1 + 3x2 - 2x3 + 0s1 + 0s2 + 0s3
Subject to
3x1 – x2 + 3x3 + s1 = 7
-2x1 + 4x2 + s2 = 12
-4x1 + 3x2 + 8x3 + s3 = 10
x1 ≥ 0, x2 ≥ 0, x3 ≥ 0 s1 ≥ 0, s2 ≥ 0, s3 ≥ 0
Cj → -1 3 -2 0 0
Basic CB XB X1 X2 X3 S1 S2
Variables
s1 0 7 3 -1 3 1 0 0 -
s2 0 12 -2 4 0 0 1 0 3→
s3 0 10 -4 3 8 0 0 1 10/3
↑
Z' = 0 1 -3 2 0 0 0 ←Δj
(R1 = R1 + R2)
s1 0 10 5/2 0 3 1 1/4 0 4→
(R2 = R2 / 4)
x2 3 3 -1/2 1 0 0 1/4 0 -
(R3 = R3 – 3R2)
s3 0 1 -5/2 0 8 0 -3/4 1 -
↑
Z' = 9 -5/2 0 0 0 3/4 0 ←Δj
(R1 = R1 / 5/2)
x1 -1 4 1 0 6/5 2/5 1/10 0
(R2 = R2 + 1/2 R1)
x2 3 5 0 1 3/5 1/5 3/10 0
(R3 = R3 + 5/2R1)
s3 0 11 0 1 11 1 -1/2 1
13
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Example 7
Max Z = 2x + 5y
x + y ≤ 600
0 ≤ x ≤ 400
0 ≤ y ≤ 300
Solution
SLPP
Max Z = 2x + 5y + 0s1 + 0s2 + 0s3
x + y + s1 = 600
x + s2 = 400
y + s3 = 300
x1 ≥ 0, y ≥ 0, s1 ≥ 0, s2 ≥ 0, s3 ≥ 0
Cj → 2 5 0 0 0
Basic Min ratio
CB XB X Y S1 S2 S3
Variables XB /Xk
s1 0 600 1 1 1 0 0 600 / 1 = 600
s2 0 400 1 0 0 1 0 -
s3 0 300 0 1 0 0 1 300 /1 = 300→
↑
Z=0 -2 -5 0 0 0 ←Δj
(R1 = R1 – R3)
s1 0 300 1 0 1 0 -1 300 /1 = 300→
s2 0 400 1 0 0 1 0 400 / 1 = 400
y 5 300 0 1 0 0 1 -
↑
Z = 1500 -2 0 0 0 5 ←Δj
x 2 300 1 0 1 0 -1
(R2 = R2 – R1)
s2 0 100 0 0 -1 1 1
y 5 300 0 1 0 0 1
14
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Z = 2100 0 0 2 0 3 ←Δj
Step 2 – Add non-negative artificial variable to the left side of each of the equations
corresponding to the constraints of the type ‘≥’ or ‘=’.
When artificial variables are added, it causes violation of the corresponding constraints. This
difficulty is removed by introducing a condition which ensures that artificial variables will
be zero in the final solution (provided the solution of the problem exists).
On the other hand, if the problem does not have a solution, at least one of the artificial
variables will appear in the final solution with positive value. This is achieved by assigning a
very large price (per unit penalty) to these variables in the objective function. Such large
price will be designated by –M for maximization problems (+M for minimizing problem),
where M > 0.
Step 3 – In the last, use the artificial variables for the starting solution and proceed with the
usual simplex routine until the optimal solution is obtained.
2.2 Worked
Examples
Example 1
Max Z = -2x1 - x2
Subject to
3x1 + x2 = 3
4x1 + 3x2 ≥ 6
x1 + 2x2 ≤ 4
and x1 ≥ 0, x2 ≥ 0
Solution
SLPP
Max Z = -2x1 - x2 + 0s1 + 0s2 - M a1 - M a2
Subject to
3x1 + x2 + a1= 3
4x1 + 3x2 – s1 + a2 = 6
15
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
x1 + 2x2 + s2 = 4
x1 , x2 , s1, s2, a1, a2 ≥ 0
Cj → -2 -1 0 0 -M -M
Basic Min ratio
CB XB X1 X2 S1 S2 A1 A2
Variables XB /Xk
a1 -M 3 3 1 0 0 1 0 3 /3 = 1→
a2 -M 6 4 3 -1 0 0 1 6 / 4 =1.5
s2 0 4 1 2 0 1 0 0 4/1=4
↑
Z = -9M 2 – 7M 1 – 4M M 0 0 0 ←Δj
x1 -2 1 1 1/3 0 0 X 0 1/1/3 =3
a2 -M 2 0 5/3 -1 0 X 1 6/5/3 =1.2→
s2 0 3 0 5/3 0 1 X 0 4/5/3=1.8
←Δj
Z = -2 – 2M
0 0 X 0
x1 -2 3/5 01 0 1/5 0 X X
x2 -1 6/5 0 1 -3/5 0 X X
s2 0 1 0 0 1 1 X X
Z = -12/5 0 0 1/5 0 X X
Solution
SLPP
Max Z = 3x1 - x2 + 0s1 + 0s2 + 0s3 - M a1
Subject to
2x1 + x2 – s1+ a1= 2
x1 + 3x2 + s2 = 3
x2 + s3 = 4
x1 , x2 , s1, s2, s3, a1 ≥ 0
16
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Cj → 3 -1 0 0 0 -M
Basic Min ratio
CB XB X1 X2 S1 S2 S3 A1
Variables XB /Xk
a1 -M 2 2 1 -1 0 0 1 2 / 2 = 1→
s2 0 3 1 3 0 1 0 0 3/1=3
s3 0 4 0 1 0 0 1 0 -
↑
Z = -2M -2M-3 -M+1 M 0 0 0 ←Δj
x1 3 1 1 1/2 -1/2 0 0 X -
s2 0 2 0 5/2 1/2 1 0 X 2/1/2 = 4→
s3 0 4 0 1 0 0 1 X -
↑
Z=3 0 5/2 -3/2 0 0 X ←Δj
x1 3 3 1 3 0 1/2 0 X
s1 0 4 0 5 1 2 0 X
s3 0 4 0 1 0 0 1 X
Z=9 0 10 0 3/2 0 X
17
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Example 3
Min Z = 2x1 + 3x2
Subject to
x1 + x2 ≥ 5
x1 + 2x2 ≥ 6
and x1 ≥ 0, x2 ≥ 0
Solution
SLPP
׳
Min Z = Max Z2- = x1 - 3x2 + 0s1 + 0s2 - M a1 - M a2
Subject to
x1 + x2 – s1+ a1= 5
x1 + 2x2 – s2+ a2= 6
x1 , x2 , s1, s2, a1, a2 ≥ 0
Cj → -2 -3 0 0 -M -M
Basic Min ratio
CB XB X1 X2 S1 S2 A1 A2
Variables XB /Xk
a1 -M 5 1 1 -1 0 1 0 5 /1 = 5
a2 -M 6 1 2 0 -1 0 1 6 / 2 = 3→
↑
׳
Z11- = M -2M + 2 -3M+3 M M 0 0 ←Δj
a1 -M 2 1/2 0 -1 1/2 1 X 2/1/2 = 4→
x2 -3 3 1/2 1 0 -1/2 0 X 3/1/2 =6
↑
׳
Z2- = M-9 (-M+1) / 2 0 M (-M+3)/2 0 X ←Δj
x1 -2 4 1 0 -2 1 X X
x2 -3 1 0 1 1 -1 X X
׳
Z11- = 0 0 1 1 X X
Example 4
Max Z =3x1 + 2x2 + x3
18
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Subject to
2x1 + x2 + x3 = 12
3x1 + 4x2 = 11
and x1 is unrestricted
x2 ≥ 0, x3 ≥ 0
Solution
SLPP
' ''
Max Z = 3(x1 - x1 ) + 2x2 + x3 - M a1 - M a2
Subject to
' ''
2(x1 - x1 ) + x2 + x3 + a1= 12
' ''
3(x1 - x1 ) + 4x2 + a2 = 11
' ''
x1 , x1 , x2 , x3, a1, a2 ≥ 0
' ''
Max Z = 3x1 - 3x1 + 2x2 + x3 - M a1 - M a2
Subject to
' ''
2x1 - 2x1 + x2 + x3 + a1= 12
' ''
3x1 - 3x1 + 4x2 + a2 = 11
' ''
x1 , x1 , x2 , x3, a1, a2 ≥ 0
Cj → 3 -3 2 1 -M -M
Basic Min ratio
CB XB X1 ' X1
'' X2 X3 A1 A2
Variables XB /Xk
a1 -M 12 2 -2 1 1 1 0 12 /2 = 6
a2 -M 11 3 -3 4 0 0 1 11/3 =3.6→
↑
Z = -23M -5M-3 5M+3 -5M-2 -M-1 0 0 ←Δj
a1 ' -M 14/3 0 0 -5/3 1 1 X 14/3/1 = 14/3→
x1 3 11/3 1 -1 4/3 0 0 X -
↑
0 -6 5/3M+1 -M-1 0 X ←Δj
x3 ' 1 14/3 0 0 -5/3 1 X X
x1 3 11/3 1 -1 4/3 0 X X
Z = 47/3 0 0 1/3 0 X X
19
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Example 5
Max Z = 8x2
Subject to
x1 - x2 ≥ 0
2x1 + 3x2 ≤ -6
and x1 , x2 unrestricted
Solution
SLPP
' ''
Max Z = 8 (x2 – x2 ) + 0s1 + 0s2 - M a1 - M a2
Subject to
' '' ' ''
(x1 - x1 ) - (x2 – x2 ) – s1+ a1= 0
' '' ' ''
-2(x1 - x1 ) - 3(x2 – x2 ) - s2 + a2 = 6
' '' ' ''
x1 , x1 , x2 , x2 , s1, a1, a2 ≥ 0
' ''
Max Z = 8x2 – 8x2 + 0s1 + 0s2 - M a1 - M a2
Subject to
' '' ' ''
x1 - x1 - x2 + x2 – s1+ a1= 0
' '' ' ''
-2x1 + 2x1 - 3x2 + 3x2 - s2 + a2 = 6
' '' ' ''
x1 , x1 , x2 , x2 , s1, a1, a2 ≥ 0
20
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
' ''
x2 = 0, x2 = 6/5
' ''
x2 = x2 – x2 = 0 – 6/5 = -6/5
Phase I – In this phase, the simplex method is applied to a specially constructed auxiliary
linear programming problem leading to a final simplex table containing a basic feasible
solution to the original problem.
Step 1 – Assign a cost -1 to each artificial variable and a cost 0 to all other variables
in the objective function.
Step 2 – Construct the Auxiliary LPP in which the new objective function Z* is to
be maximized subject to the given set of constraints.
Step 3 – Solve the auxiliary problem by simplex method until either of the
following three possibilities do arise
i. Max Z* < 0 and atleast one artificial vector appear in the optimum basis
at a positive level (Δj ≥ 0). In this case, given problem does not possess
any feasible solution.
ii. Max Z* = 0 and at least one artificial vector appears in the optimum
basis at a zero level. In this case proceed to phase-II.
iii. Max Z* = 0 and no one artificial vector appears in the optimum basis. In
this case also proceed to phase-II.
Phase II – Now assign the actual cost to the variables in the objective function and a zero cost
to every artificial variable that appears in the basis at the zero level. This new objective
function is now maximized by simplex method subject to the given constraints.
Simplex method is applied to the modified simplex table obtained at the end of phase-I, until
an optimum basic feasible solution has been attained. The artificial variables which are non-
basic at the end of phase-I are removed.
21
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
and x1 ≥ 0, x2 ≥ 0
Solution
Standard LPP
Max Z = 3x1 - x2
Subject to
2x1 + x2 – s1+ a1= 2
x1 + 3x2 + s2 = 2
x2 + s3 = 4
x1 , x2 , s1, s2, s3,a1 ≥ 0
Auxiliary LPP
Max Z* = 0x1 - 0x2 + 0s1 + 0s2 + 0s3 -1a1
Subject to
2x1 + x2 – s1+ a1= 2
x1 + 3x2 + s2 = 2
x2 + s3 = 4
x1 , x2 , s1, s2, s3,a1 ≥ 0
Phase I
Cj → 0 0 0 0 0 -1
Basic Min ratio
CB XB X1 X2 S1 S2 S3 A1
Variables XB /Xk
a1 -1 2 2 1 -1 0 0 1 1→
s2 0 2 1 3 0 1 0 0 2
s3 0 4 0 1 0 0 1 0 -
↑
Z* = -2 -2 -1 1 0 0 0 ←Δj
x1 0 1 1 1/2 -1/2 0 0 X
s2 0 1 0 5/2 1/2 1 0 X
s3 0 4 0 1 0 0 1 X
Z* = 0 0 0 0 0 0 X ←Δj
Since all Δj ≥ 0, Max Z* = 0 and no artificial vector appears in the basis, we proceed to phase II.
Phase II
Cj → 3 -1 0 0 0
22
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Z=6 0 10 0 3 0 ←Δj
Example 2
Max Z = 5x1 + 8x2
Subject to
3x1 + 2x2 ≥ 3
x1 + 4x2 ≥ 4
x1 + x2 ≤ 5
and x1 ≥ 0, x2 ≥ 0
Solution
Standard LPP
Max Z = 5x1 + 8x2
Subject to
3x1 + 2x2 – s1+ a1 = 3
x1 + 4x2 – s2+ a2 = 4
x1 + x2 + s3 = 5
x1 , x2 , s1, s2, s3, a1, a2 ≥ 0
Auxiliary LPP
Max Z* = 0x1 + 0x2 + 0s1 + 0s2 + 0s3 -1a1 -1a2
Subject to
3x1 + 2x2 – s1+ a1 = 3
x1 + 4x2 – s2+ a2 = 4
x1 + x2 + s3 = 5
x1 , x2 , s1, s2, s3, a1, a2 ≥ 0
Phase I
Cj → 0 0 0 0 0 -1 -1
23
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Z* = 0 0 0 0 0 0 X X ←Δj
Since all Δj ≥ 0, Max Z* = 0 and no artificial vector appears in the basis, we proceed to phase II.
Phase II
Cj → 5 8 0 0 0
Basic Min ratio
24
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
CB XB X1 X2 S1 S2 S3
Variables XB /Xk
x1 5 2/5 1 0 -2/5 1/5 0 2→
x2 8 9/10 0 1 1/10 -3/10 0 -
s3 0 37/10 0 0 3/10 1/10 1 37
↑
Z = 46/5 0 0 -6/5 -7/5 0 ←Δj
s2 0 2 5 0 -2 1 0 -
x2 8 3/2 3/2 1 -1/2 0 0 -
s3 0 7/2 -1/2 0 1/2 0 1 7→
↑
Z = 12 7 0 -4 0 0 ←Δj
s2 0 16 3 0 0 1 2
x2 8 5 1 1 0 0 1/2
s1 0 7 -1 0 1 0 2
Z = 40 3 0 0 0 4
25
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Example 3
Max Z = -4x1 - 3x2 - 9x3
Subject to
2x1 + 4x2 + 6x3 ≥ 15
6x1 + x2 + 6x3 ≥ 12
and x1 ≥ 0, x2 ≥ 0, x3 ≥ 0
Solution
Standard LPP
Max Z = -4x1 - 3x2 - 9x3
Subject to
2x1 + 4x2 + 6x3 - s1+ a1= 15
6x1 + x2 + 6x3 - s2 + a2 = 12
x1 , x2 , s1, s2, a1, a2 ≥ 0
Auxiliary LPP
Max Z* = 0x1 - 0x2 - 0x3 + 0s1 + 0s2 -1a1 -1a2
Subject to
2x1 + 4x2 + 6x3 - s1+ a1= 15
6x1 + x2 + 6x3 - s2 + a2 = 12
x1 , x2 , s1, s2, a1, a2 ≥ 0
Phase I
Cj → 0 0 0 0 0 -1 -1
Basic Min ratio
26
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
CB XB X1 X2 X3 S1 S2 A1 A2
Variables XB /Xk
a1 -1 15 2 4 6 -1 0 1 0 15/6
a2 -1 12 6 1 6 0 -1 0 1 2→
↑
Z* = -27 -8 -5 -12 1 1 0 0 ←Δj
a1 -1 3 -4 3 0 -1 1 1 X 1→
x3 0 2 1 1/6 1 0 -1/6 0 X 12
↑
Z* = -3 4 -3 0 1 -1 0 X ←Δj
x2 0 1 -4/3 1 0 -1/3 1/3 X X
x3 0 11/6 22/18 0 1 1/18 -4/18 X X
Z* = 0 0 0 0 0 0 X X
27
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Since all Δj ≥ 0, Max Z* = 0 and no artificial vector appears in the basis, we proceed to phase II.
Phase II
Cj → -4 -3 -9 0 0
Basic Min ratio
CB XB X1 X2 X3 S1 S2
Variables XB /Xk
x2 -3 1 -4/3 1 0 -1/3 1/3 -
x3 -9 11/6 22/18 0 1 1/18 -4/18 3/2→
↑
Z = -39/2 -3 0 0 1/2 1 ←Δj
x2 -3 3 0 1 12/11 -3/11 1/11
x1 -4 3/2 1 0 18/22 1/22 -4/22
Example 4
Min Z = 4x1 + x2
Subject to
3x1 + x2 = 3
4x1 + 3x2 ≥ 6
x1 + 2x2 ≤ 4
and x1 ≥ 0, x2 ≥ 0
Solution
Standard LPP
'
Min Z = Max Z = – 4x1 – x2
Subject to
3x1 + x2 + a1 = 3
4x1 + 3x2 – s1+ a2 = 6
x1 + 2x2 + s2 = 4
x1 , x2 , s1, s2, a1, a2 ≥ 0
Auxiliary LPP
Max Z* = 0x1 – 0x2 + 0s1 + 0s2 –1a1 –1a2
Subject to
3x1 + x2 + a1 = 3
4x1 + 3x2 – s1+ a2 = 6
x1 + 2x2 + s2 = 4
x1 , x2 , s1, s2, a1, a2 ≥ 0
28
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Phase I
Cj → 0 0 0 0 -1 -1
Basic Min ratio
CB XB X1 X2 S1 S2 A1 A2
Variables XB /Xk
a1 -1 3 3 1 0 0 1 0 1→
a2 -1 6 4 3 -1 0 0 1 6/4
s2 0 4 1 2 0 1 0 0 4
↑
Z* = -9 -7 -4 1 0 0 0
x1 0 1 1 1/3 0 0 X 0 3
a2 -1 2 0 5/3 -1 0 X 1 6/5→
s2 0 3 0 5/3 0 1 X 0 9/5
↑
Z* = -2 0 -5/3 1 0 X 0
x1 0 3/5 1 0 1/5 0 X X
x2 0 6/5 0 1 -3/5 0 X X
s2 0 1 0 0 1 1 X X
Z* = 0 0 0 0 0 X X
Since all Δj ≥ 0, Max Z* = 0 and no artificial vector appears in the basis, we proceed to phase II.
Phase II
Cj → -4 -1 0 0
Basic Min ratio
CB XB X1 X2 S1 S2
Variables XB /Xk
x1 -4 3/5 1 0 1/5 0 3
x2 -1 6/5 0 1 -3/5 0 -
s2 0 1 0 0 1 1 1→
↑
'
Z = -18/5 0 0 -1/5 0 ←Δj
x1 -4 2/5 1 0 0 -1/5
x2 -1 9/5 0 1 0 3/5
s1 0 1 0 0 1 1
'
Z = -17/5 0 0 0 1/5 ←Δj
29
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
'
Therefore the solution is Max Z = -17/5
Min Z = 17/5, x1 = 2/5, x2 = 9/5
Example 5
Min Z = x1 –2x2– 3x3
Subject to
–2x1 + x2 + 3x3= 2
2x1 + 3x2 + 4x3= 1
and x1 ≥ 0, x2 ≥ 0
Solution
Standard LPP
'
Min Z = Max Z = –x1 +2x2+ 3x3
Subject to
–2x1 + x2 + 3x3 + a1 = 2
2x1 + 3x2 + 4x3+ a2 = 1
x1 , x2 , a1, a2 ≥ 0
Auxiliary LPP
Max Z* = 0x1 + 0x2 + 0x3 –1a1 –1a2
Subject to
–2x1 + x2 + 3x3 + a1 = 2
2x1 + 3x2 + 4x3+ a2 = 1
x1 , x2 , a1, a2 ≥ 0
Phase I
Cj→ 0 0 0 -1 -1
Basic Min Ratio
CB XB X1 X2 X3 A1 A2
Variables XB / XK
a1 -1 2 -2 1 3 1 0 2/3
a2 -1 1 2 3 4 0 1 1/4→
↑
Z* = -3 0 -4 -7 0 0 ←Δj
a1 -1 5/4 -7/4 -5/4 0 1 X
x3 0 1/4 1/2 3/4 1 0 X
30
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Since for all Δj ≥ 0, optimum level is achieved. At the end of phase-I Max Z* < 0 and one of the
artificial variable a1 appears at the positive optimum level. Hence the given problem does not
posses any feasible solution.
3.1.1 Degeneracy
The concept of obtaining a degenerate basic feasible solution in a LPP is known as degeneracy.
The degeneracy in a LPP may arise At the initial stage when at least one basic variable
is zero in the initial basic feasible solution.
At any subsequent iteration when more than one basic variable is eligible to leave
the basic and hence one or more variables becoming zero in the next iteration and
the problem is said to degenerate. There is no assurance that the value of the objective
function will improve, since the new solutions may remain degenerate. As a result, it is
possible to repeat the same sequence of simplex iterations endlessly without
improving the solutions. This concept is known as cycling or circling.
Example 1
Solution
Standard LPP
Max Z = 3x1 + 9x2 + 0s1 + 0s2
Subject to
x1 + 4x2 + s1 = 8
x1 + 2x2 + s2 = 4
x1 , x2 , s1, s2 ≥ 0
Cj→ 3 9 0 0
Basic
CB XB X1 X2 S1 S2 XB / XK S 1 / X2
Variables
31
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
s1 0 8 1 4 1 0 1/4
s2 0 4 1 2 0 1 0/2→
↑
Z=0 -3 -9 0 0 ←Δj
s1 0 0 -1 0 1 -1
x2 9 2 1/2 1 0 1/2
Note – Since a tie in minimum ratio (degeneracy), we find minimum of s1 /xk for these rows
for which the tie exists.
Example 2
Max Z = 2x1 + x2
Subject to
4x1 + 3x2 ≤ 12
4x1 + x2 ≤ 8
4x1 - x2 ≤ 8
and x1 ≥ 0, x2 ≥ 0
Solution
Standard LPP
Max Z = 2x1 + x2 + 0s1 + 0s2 + 0s3
Subject to
4x1 + 3x2 + s1 = 12
4x1 + x2 + s2 = 8
4x1 - x2 + s3 = 8
x1 , x2 , s1, s2, s3 ≥ 0
32
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Cj → 2 1 0 0 0
Basic
CB XB X1 X2 S1 S2 S3 XB / XK S1 / X1 S2 / X1
Varibles
s1 0 12 4 3 1 0 0 12/4=3
s2 0 8 4 1 0 1 0 8/4=2 4/0=0 1/4
s3 0 8 4 -1 0 0 1 8/4=2 4/0=0 0/4=0→
↑
Z=0 -2 -1 0 0 0 ←Δj
1 0 4 0 4 1 0 -1 4/4=1
s2 0 0 0 2 0 1 -1 0→
x1 2 2 1 -1/4 0 0 1/4 -
↑
Z=4 0 -3/2 0 0 1/2 ←Δj
s1 0 4 0 0 1 -2 1 0→
x2 1 0 0 1 0 1/2 -1/2 -
x1 2 2 1 0 0 1/8 1/8 16
↑
Z=4 0 0 0 3/4 -1/4 ←Δj
s3 0 4 0 0 1 -2 1
x2 1 2 0 1 1/2 -1/2 0
x1 2 3/2 1 0 -1/8 3/8 0
Example
Max Z = 3x1 +2x2
Subject to
2x1 + x2 ≤ 2
3x1 + 4x2 ≥
12 and x1 ≥ 0, x2
≥0
33
Department Of Mechanical Engineering, GMIT, Bharathinagara
OPERATION RESEARCH [15ME81] 2019-2020
Solution
Standard LPP
Max Z = 3x1 +2 x2 + 0s1 + 0s2 – Ma1
Subject to
2x1 + x2 + s1 = 2
3x1 + 4x2 - s2 + a1 = 12
x1 , x2 , s1, s2, s3 ≥ 0
Cj → 3 2 0 0 -M
Basic Min Ratio
CB XB X1 X2 S1 S2 A1
Variables XB / XK
s1 0 2 2 1 1 0 0 2/1=2→
a1 -M 12 3 4 0 -1 1 12/4=3
↑
Z= -12M -3M-3 -4M-2 0 M 0 ←Δj
x2 2 2 2 1 1 0 0
a1 -M 4 -5 0 -4 -1 1
34
Department Of Mechanical Engineering, GMIT, Bharathinagara