Operations Research PDF
Operations Research PDF
Faculty of Commerce
Department of Statistics
Operations Research
LECTURE NOTES
Prepared by
0
List of Contents
1
Chapter One
What Is Operations Research?
(1-1): Introduction:
2
To formalize this observation, the decision variables;
namely, the width and height of the rectangle are identified by
defining the width and height algebraically as:
1. ( )
2.
( )
3
The general OR model can be organized in the following
general format:
Subject to:
Constraints
4
Step 3. Formulate a Mathematical Model of the Problem
Step 4. Verify the Model and Use the Model for Prediction
If the decision maker(s) has accepted the study, the analyst aids
in implementing the recommendations. The system must be
constantly monitored (and updated dynamically as the
environment changes) to ensure that the recommendations are
enabling decision maker(s) to meet her/his/their objectives.
An objective function
6
An optimal solution to the model is the identification of a
set of variable values which are feasible (satisfy all the
constraints) and which lead to the optimal value of the
objective function.
Availability (tons)
M1 6 4 24
M2 1 2 6
Profit per 5 4
ton($1000)
8
Let Z represent the total daily profit (in thousands of dollars),
the objective (or goal) of Reddy Mikks is expressed as
( ) ( )
In a similar manner,
( )
( )
( )
9
The second restriction limits the daily demand of interior paint
to 2 tons – that is,
( )
10
Mine Cost per day £ Production (tons/day)
Solution
Variables
Constraints
High 6x + 1y ≥ 12
Medium 3x + 1y ≥ 8
Low 4x + 6y ≥ 24
11
days per week constraint - we cannot work more than a
certain maximum number of days a week e.g. for a 5 day
week we have
x≤5
y≤5
Objective
Subject to
6x + y ≥ 12
3x + y ≥ 8
4x + 6y ≥ 24
X ≤5
Y ≤5
x, y ≥ 0
12
LINEAR PROGRAMMING
Work scheduling
Capital budgeting
Financial planning
Distribution
o Inventory model
o Financial models
o Work scheduling
1.4 FORMULATING LP
Answer
14
Decision variables completely describe the decisions to be
made (in this case, by Giapetto). Giapetto must decide how
many soldiers and trains should be manufactured each week.
With this in mind, we define:
15
assume nonnegative values (Giapetto can not manufacture
negative number of soldiers or trains!)
(LP) model
Dorian makes luxury cars and jeeps for high-income men and
women. It wishes to advertise with 1 minute spots in comedy
shows and football games. Each comedy spot costs $50K and is
seen by 7M high-income women and 2M high-income men.
Each football spot costs $100K and is seen by 2M high-income
women and 12M high-income men. How can Dorian reach
28M high-income women and 24M high- income men at the
16
least cost?
Answer
Min z = 50 x1 + 100 x2
st 7x1 + 2x2 ≥ 28
2x1 +12x2 ≥ 24
x1, x2 ≥ 0
17
Profit/Kg P1 P2 Total available Machine
Machine 1 3 2 600
Machine 2 3 5 800
Machine 3 5 6 1100
Solution:
Let x1 = amount of P1
x2 = amount of P2
30x1 + 40x2
x1 ≥ 0 ; x2 ≥ 0
Maximize
30x1 + 40x2
Subject to:
x1≥ 0, x2 ≥ 0
Solution:
Let us define x1 and x2 are the mills A and B. Here the objective
is to minimize the cost of the machine runs and to satisfy the
contract order. The linear programming problem is given by
Minimize
2000x1 + 1500x2
Subject to:
6x1 + 2x2≥ 8
2x1 + 4x2≥ 12
4x1 + 12x2≥ 24
x1≥ 0, x2≥ 0
20
Chapter two
General Linear Programming Problem
Where cj, bi, and aij (i = 1, 2 , .... m; j = 1.2 .... , n) are constants
determined from the technology or the problem, and xj are the
decision variables. Only one sign (≤, =, or ≥) holds for each
constraint. Although all variables are declared nonnegative, the
preceding section shows that every unrestricted variable can be
converted equivalently to nonnegative variables. The non-
negative restriction is essential for the development of the
solution method for linear programming.
Subject to :
∑∑
22
2- All constraints are of the (≤) type.
Minimize Z = c1 x1 + c2 x2 + ..... + cn xn
Is équivalent to
a 1 x 1 + a2 x 2 ≤ b and a1 x 1 + a2 x 2 ≥ b
│a1x1 + a2 x2 │ ≤ b
is equivalent to
a 1 x 1 + a2 x 2 ≥ - b and a1x1 + a2 x2 ≤ b
Similarly, for q ≥ 0.
│p1x1 + p2 x2 │ ≥ q
is equivalent to either
p1x1 + p2 x2 ≥ q or p 1 x 1 + p2 x 2 ≤ - q
24
* Example:
Subject to
x1 + x2 + 3x2 ≤ 40
x1 + 9 x2 - 7x3 ≥ 50
5x1 + 3x2 = 20
x1 ≥ 0 , x2 ≥ 0
x3 is unconstrained in sign
x1 + 9 x2 - 7x3 ≥ 50
is equivalent to
5x1 + 3x2 = 20
is equivalent to
5x1 + 3x2 20
25
And -5x1 - 3x2 -20
is equivalent to
̿̿̿ ̅̅̅
Subject to :
x1 + x2 + 3 (̿̿̿ ̅̅̅) ≤ 40
5x1 + 3x2 ≤ 20
A1x1 + a2x2 ≥ b , b ≥ 0
p1x1 + p2x2 ≤ q , q ≥ 0
is changed to
P1x1 + p2x2 + S2 = q
Where S2 ≥ 0
Subject to :
∑∑
Subject to :
∑∑
Subject to :
∑∑
Subject to :
∑∑
29
The standard form basically reduces the linear program to
a set of (m) equations in (m + n) unknowns. A solution to these
equations is of interest only if it is feasible, that is , if it satisfies
the nonnegativity constraints xj ≥ 0 and Si ≥ 0. for all i and j . the
optimum solution is then given by the feasible solution that
maximizes x0.
30
Canonical form Standard form
Note that :
1- If │ax1 ± bx2 │≤ C
31
then : - C ≤ y ≤ C
c ≤ ax1 ± bx2 ≤ -c
32
o If the constraint is of ‘≥’ type, we subtract a variable
called Surplus variable and convert it to an
equation. For example
2x1 +x2 ≥ 15
2x1 +x2 – s2 = 15
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 = ̿̿̿ - ̅̅̅ ׳where ̿̿̿, ̅̅̅ ≥ 0
Therefore the inequality takes the form of equation as
x1 + (̿̿̿ - ̅̅̅) + s1 = 3
Using the above steps, we can write the GLPP in the form of
SLPP.
Write the Standard LPP (SLPP) of the following
Example
Maximize Z = 3x1 + x2
Subject to
2 x 1 + x2 ≤ 2
3 x1 + 4 x2 ≥ 12
and x1 ≥ 0, x2 ≥ 0
SLPP
33
Maximize Z = 3x1 + x2
Subject to
2 x 1 + x 2 + s1 = 2
3 x1 + 4 x2 – s2 = 12
x1≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0
Example
Minimize Z = 4x1 + 2 x2
Subject to
3x1 + x2 ≥ 2
x1 + x2 ≥ 21
x1 + 2x2 ≥ 30
and x1 ≥ 0, x2 ≥ 0
SLPP
Minimize Z = 4x1 + 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
Minimize Z = x1 + 2 x2 + 3x3
Subject to
2x1 + 3x2 + 3x3 ≥ – 4
3x1 + 5x2 + 2x3 ≤ 7
and x1 ≥ 0, x2 ≥ 0, x3 is unrestricted in sign
34
SLPP put : ̿̿̿ ̅̅̅
Example :
Subject to :
x1 + 2x2 + x3 ≤ 35
3x1 – 2x2 = 5
x1 , x3 ≥ 0
x2 is unrestricted in sign.
Required :
35
Solution :
̿̿̿ ̅̅̅
Subject to :
x1 + 2(̿̿̿ ̅̅̅) + x3
( (̿̿̿ ̅̅̅) )
(̿̿̿ ̅̅̅)
(̿̿̿ ̅̅̅)
( (̿̿̿ ̅̅̅) )
36
Since there is no intersection between the two arrows, (i.e.,
there is no feasible solution between the two arrows) , then ,
we have to take only one of the two sub constraints
(̿̿̿ ̅̅̅)
Or :
- (̿̿̿ ̅̅̅)
Subject to:
( (̿̿̿ ̅̅̅) )
(̿̿̿ ̅̅̅)
(̿̿̿ ̅̅̅)
( (̿̿̿ ̅̅̅) )
37
We have to take only one of the two sub constraints
(̿̿̿ ̅̅̅)
Or :
(̿̿̿ ̅̅̅)
38
It is a basic feasible solution which has exactly ‘m’ positive
xi, where i=1, 2 … m. In other words all ‘m’ basic variables
are positive and remaining ‘n’ variables are zero.
Optimum basic feasible solution
A basic feasible solution is said to be optimum if it optimizes
(max / min) the objective function.
Introduction to Simplex Method
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.
(1) If all the constraints of the LPM are in the relation form
less than or equal to, i.e., (≤) , then we must use the
ordeinary simplex method.
a. The M – technique.
: :
42
Note that, the special arrangement of the above tableau
provides useful information. The columns associated with the
starting basic variable always appear immediately to the left of
the solution column, and their constraints coefficients will
constitute an identity matrix.
3- Solution improvement :
43
b) Determine the pivot row (Leaving variable):
45
Feasibility condition. The leaving variable is the basic
variable corresponding to the smallest ratio of the current
value of the basic variables to the positive constraint coefficient
of the entering variable.
46
The following examples are now introduced to sum up the
basic features of the ordinary simplex method.
Computational Procedure of Simplex Method
Consider an example
Maximize Z = 3x1 + 2x2
Subject to
x 1 + x2 ≤ 4
x1 – x2 ≤ 2
and x1 ≥ 0, x2 ≥ 0
Solution
Step 1 Write the given GLPP in the form of SLPP
Maximize Z = 3x1 + 2x2 + 0s1 + 0s2
Subject to
x1 + x2+ s1= 4
x1 – x2 + s2= 2
x1≥ 0, x2 ≥ 0, s1 ≥ 0, s2 ≥ 0
Step 2 Present the constraints in the matrix form
x 1 + x 2+ s 1 = 4
x 1 – x 2 + s2 = 2
( )[ ] * +
47
Step 3 – Construct the starting simplex table using the notations
(TABLEAU 1)
Coef.Z Cj 3 2 0 0 Min
Ratio.
(B.V) CB XB X1 X2 S1 S2
0 4 1 1 1 0 4
0 2 1 -1 0 1 2
Zj=CBXj Z=0 0 0 0 0
Δj= Cj - Zj 3 2 0 0
Max Z= CBXB =0
49
Any new table will form:
The row that replaces the outgoing row (V)
= pivot row / pivot element
The remaining rows:
New row = old row – corresponding element at pivotcolumn * V
(TABLEAU 2)
Coef.Z Cj 3 2 0 0 Min
Ratio.
(B.V) CB XB X1 X2 S1 S2
0 2 0 2 1 -1 1
3 2 1 -1 0 1 -2
Zj=CBXj Z=6 3 -3 0 3
Δj= Cj - Zj 0 5 0 -3
Max Z= CBXB =6
50
(TABLEAU 3)
Coef.Z Cj 3 2 0 0 Min
Ratio.
(B.V) CB XB X1 X2 S1 S2
2 1 0 1 ⁄ ⁄
3 3 1 0 ⁄ ⁄
Zj=CBXj Z=11 3 1 ⁄ ⁄
Δj= Cj - Zj 0 0 ⁄ ⁄
Max Z= CBXB =11
( ) [ ]
[ ]
(TABLEAU 1)
Coef.Z Cj 5 7 0 0 0 Min
Ratio.
(B.V) CB XB X1 X2 S1 S2 S3
0 4 1 1 1 0 0 4
0 24 3 -8 0 1 0 -3
0 35 10 7 0 0 1 5
Zj=CBXj Z=0 0 0 0 0 0
Δj= Cj - Zj 5 7 0 0 0
Max Z= CBXB =0
52
(TABLEAU 2)
Coef.Z Cj 5 7 0 0 0 Min
Ratio.
(B.V) CB XB X1 X2 S1 S2 S3
7 4 1 1 1 0 0
0 56 11 0 8 1 0
0 7 3 0 -7 0 1
Zj=CBXj Z=0 7 7 7 0 0
Δj= Cj - Zj -2 0 -7 0 0
Max Z= CBXB =0
53
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
(TABLEAU 1)
Coef.Z Cj 2 -3 1 0 0 0 Min
Ratio.
(B.V) CB XB X Y Z S1 S2 S3
0 6 3 6 1 1 0 0 2
0 4 4 2 1 0 1 0 1
0 3 1 -1 1 0 0 1 3
Zj=CBXj Z=0 0 0 0 0 0 0
Δj= Cj - Zj 2 -3 1 0 0 0
Max Z= CBXB =0
(TABLEAU 2)
54
Coef.Z Cj 2 -3 1 0 0 0 Min
Ratio.
(B.V) CB XB X Y Z S1 S2 S3
0 3 0 1 0
12
2 1 1 0 0
4
0 2 0 0 1
Zj=CBXj Z=2 2 1 0 0
Δj= Cj - Zj 0 -4 0 0
Max Z= CBXB =0
(TABLEAU 3)
Coef.Z Cj 2 -3 1 0 0 0 Min
Ratio.
(B.V) CB XB X Y Z S1 S2 S3
0 0 5 0 1
2 1 1 0 0
1 0 -2 1 0
55
Zj=CBXj Z= 2 0 1 0
Δj= Cj - Zj 0 -3 0
Max Z= CBXB =0
56
( ) [ ]
[ ]
(TABLEAU 1)
Coef.Z Cj 1 -3 2 0 0 0 Min
Ratio.
(B.V) CB XB S1 S2 S3
0 7 3 -1 3 1 0 0 _
0 12 -2 44 0 0 1 0 3
0 10 -4 3 8 0 0 1
Zj=CBXj Z=0 0 0 0 0 0 0
Δj= Cj - Zj 1 -3 2 0 0 0
Max Z= CBXB =0
(TABLEAU 2)
Coef.Z Cj 1 -3 2 0 0 0 Min
Ratio.
(B.V) CB XB S1 S2 S3
0 10 3 1 0
4
57
-3 3 0 0 0
-
0 1 8 0 1
Zj=CBXj Z -3 0 0 0
=-9
Δj= Cj - Zj 0 0 0
Max Z= CBXB =-9
(TABLEAU 3)
Coef.Z Cj 1 -3 2 0 0 0 Min
Ratio.
(B.V) CB XB S1 S2 S3
1 4 1 0 0
-3 0 1
5
0 0 0 11 1
11
Zj=CBXj Z= 1 -3
Δj= Cj - Zj 0 0 0
Max Z= CBXB =-11
Since all Δj , optimal basic feasible solution is obtained
Therefore the solution is MinZ = -11, x1 = 4, x2 = 5, x3 = 0
Exercise
58
Solve by simplex method
1. Maximize Z = 5x1 + 3x2
Subject to
3x1 + 5x2 ≤ 15
5x1 + 2x2 ≤ 10
and x1 ≥ 0, x2 ≥ 0
[Ans. Max Z = 235/19, x1= 20/19, x2= 45/19]
2. Maximize Z = 5x1 + 7x2
Subject to
x1 + x 2 ≤ 4
3x1 - 8x2 ≤ 24
10x1 + 7x2 ≤ 35
and x1 ≥ 0, x2 ≥ 0
[Ans. Max Z = 28, x1= 0, x2= 4]
3. Maximize Z = 2x1 + 4x2 + x3+ x4
Subject to
x1 + 3x2 + x4 ≤ 4
2x1 + x2 ≤ 3
x2+ 4x3 + x4 ≤ 3
and x1 ≥ 0, x2 ≥ 0, x3 ≥ 0, x4 ≥ 0
[Ans. Max Z = 13/2, x1= 1, x2= 1, x3= 1/2, x4= 0]
4. Maximize Z = 7x1 + 5x2
Subject to
-x1 - 2x2 ≥-6
4x1 + 3x2 ≤ 12
and x1 ≥ 0, x2 ≥ 0
[Ans. Max Z = 21, x1= 3, x2= 0]
59
Computational Procedure of Big – M Method (Charne’s
Penalty Method)
Step 1 – Express the problem in the standard form.
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. An artificial variable is added to form
a starting solution similar to the all slack basic solution
.However ,the artificial variables are not part of the original
problem ,and a modeling (trick) is needed to force them to zero
by the time ,the optimum iteration is reached . The desired goal
is achieved by penalizing these variables in the objective
function using the following rule :
Penalty rule for artificial variables
Given M , a sufficiently large positive value (mathematically
,M ) ,The objective coefficient of an artificial represents an
appropriate penalty if :
60
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.
Subject to:
2 x 1 + x2 2
3x1 + 4x2 ≥ 12
x1 ≥ 0, x2 ≥ 0
Solution :
Max Z = 3x1+2x2 +0 -0
61
Subject to:
2x1 + x2 + s1= 2
x1 , x2 , s1, s2 ≥ 0
Max Z = 3x1+2x2 +0 + 0 -M
2 x1 + x2 + s1 =2
( ) * +
[ ]
Step 3 : construct the initial solution tableau by considering
that the slack variable corresponding to the constraint in the
form (≤) is the basic variables, and the artificial variable
corresponding to the constraint are in either the form (= or ≥) is
62
the basic variable for these constraints. Then, we have the
following tableau:
Coef.Z Cj 3 2 0 0 -M Min
Ratio.
(B.V) CB XB S1 S2 a1
0 2 2 1 1 0 0 2
-M 3 4 0 4
12 -1 1
1st iprovement :
2nd tableau
Coef.Z Cj 3 2 0 0 -M Min
Ratio.
(B.V) CB XB S1 S2 a1
2 2 2 1 1 0 0
-M -5 0 -4
4 -1 1
63
Δj= Cj - Zj -1-5M 0 -2-4M -M 0
Max Z=4-4M
Example :
Z = 4x1 + x2 minimize
Subject to :
64
3x1 + x2 = 3
4x1 + 3x2 ≥ 6
x1 + 2x2 ≤ 4
x 1 , x2 ≥ 0
Solution :
Z = 4x1 + x2 (Min)
Subject to :
3x1 + x2 = 3
4x1 + 3x2 - S1 = 6
x1 + 2x2 + S2 = 4
x 1 , x 2 , S 1 , S2 ≥ 0
65
variable with (M) coefficient in case of minimization value
of f (x), where M > 0 : then , we have the following model :
Subject to :
3x1 + x2 + = 3
4x1 + 3x2 – S1 + = 6
x1 + 2x2 + S2 = 4
x1 ≥ 0 for i = 1 : 6
( ) [ ]
[ ]
Step 3 : construct the initial solution tableau by considering
that the slack variable corresponding to the constraint in the
form (≤) is the basic variables, and the artificial variable
corresponding to the constraint are in either the form (= or ≥) is
the basic variable for these constraints.
Coef.Z Cj 4 1 0 0 M M Min
Ratio.
(B.V) CB XB S1 S2
66
M 3 3 1 0 0 1 0 1
M 4 3 -1 0 ⁄
6 0 1
0 1 2 0 1 0 4
4 0
Zj=CBXj Z= 7M 4M -M 0 M M
Δj= Cj - Zj M 0 M
Min Z= 9M
2nd tableau
Coef.Z Cj 4 1 0 0 M M Min
Ratio.
(B.V) CB XB S1 S2
4 1 1 ⁄ 0 0 ⁄ 0
3
M 0 ⁄ -1 0 ⁄
2 ⁄ 1
0 0 ⁄ 0 1 0 ⁄
3 ⁄
Zj=CBXj Z= 4 ⁄ -M 0 ⁄ M
67
Δj= Cj - Zj 0 M 0 0
⁄ ⁄
Min Z= 4+2M
3nd tableau
Coef.Z Cj 4 1 0 0 M M Min
Ratio.
(B.V) CB XB S1 S2
4 ⁄ 1 0 ⁄ 0 ⁄ ⁄
3
1 0 1 ⁄ 0
⁄ ⁄ ⁄ -
0 0 0 1 1 -1 1
1 1
Zj=CBXj ⁄ 4 1 ⁄ 0 ⁄ ⁄
Δj= Cj - Zj 0 0 0
⁄ ⁄ ⁄
Min Z= ⁄
68
4th tableau
Coef.Z Cj 4 1 0 0 M M Min
Ratio.
(B.V) CB XB S1 S2
4 ⁄ 1 0 0 ⁄ ⁄ 0
1 0 1 0 ⁄
⁄ ⁄ 0
0 0 0 1 1 -1
1 1
Zj=CBXj ⁄ 4 1 0 ⁄ ⁄ 0
Δj= Cj - Zj 0 0 0 M
⁄ ⁄
Min Z= ⁄
⁄ ⁄ ⁄
Example
Max Z =3x1 + 2x2 + x3
Subject to:
2x1 + x2 + x3 = 12
\3x1 + 4x2 = 11
x1 is unrestricted, x2 ≥ 0, x3 ≥ 0
Solution :
̿̿̿
̅̅̅
( ) * +
[ ]
Step 3 : construct the initial solution tableau.
Coef.Z Cj 3 -3 2 1 -M -M Min
Ratio.
(B.V) CB XB ̿̿̿ ̅̅̅ X3
-M 12 2 -2 1 1 1 0 6
-M 3 -3 4 0 ⁄
11 0 1
71
2nd tableau
Coef.Z Cj 3 -3 2 1 -M - Min
M Ratio.
-M ⁄ 0 0 ⁄ 1 1 x ⁄
̿̿̿ 3 1 -1 ⁄ 0 -
⁄ 0 x
Zj=CBXj 3 -3 -M -M x
Δj= Cj - Zj 0 0 1+M 0 x
Max Z= -23M
3nd tableau
Coef.Z Cj 3 -3 2 1 -M - Min
M Ratio.
1 ⁄ 0 0 ⁄ 1 x x ⁄
72
̿̿̿ 3 1 -1 ⁄ 0 -
⁄ x x
Zj=CBXj ⁄ 3 -3 ⁄ 1 x x
Δj= Cj - Zj 0 0 0 x x
MaxZ= ⁄
̿̿̿̿ ̅̅̅̅
Then ( )
Example
Max Z = 8x2
Subject to
x1 - x2 ≥ 0
2x1 + 3x2 ≤ - 6
x1 , x2 unrestricted
Solution :
̿̿̿
̅̅̅
̿̿̿
( ) ̅̅̅ * +
[ ]
Step 3 : construct the initial solution tableau.
74
Starting (or initial ) solution tableau:
Coef.Z Cj 0 0 8 -8 0 0 -M -M Min
Ratio.
(B.V) CB XB ̿̿̿ ̅̅̅ ̿̿̿ ̅̅̅
-M 0 1 -1 -1 1 -1 0 1 0 0
-M -2 2 -3 3 2
6 0 -1 0 1
2nd tableau
Coef.Z Cj 0 0 8 -8 0 0 -M -M Min
Ratio.
(B.V) CB XB ̿̿̿ ̅̅̅ ̿̿̿ ̅̅̅
̅̅̅ -8 0 1 -1 -1 1 -1 0 x 0 -
-M -5 5 0 0 6/5
6 3 -1 x 1
75
3nd tableau
Coef.Z Cj 0 0 8 -8 0 0 -M -M Min
Ratio.
(B.V) CB XB ̿̿̿ ̅̅̅ ̿̿̿ ̅̅̅
̅̅̅ -8 ⁄ 0 0 -1 1 ⁄ ⁄ X X
̅̅̅ 0 ⁄ -1 1 0 0 ⁄ ⁄ X X
Zj=CBXj ⁄ 0 0 8 -8 ⁄ ⁄ X X
Δj= Cj - Zj 0 0 0 0 x x
⁄ ⁄
Max Z=
Then ( ) , ( )
76
Exercise
Solve by Big-M method
1. Min Z = 4x1 + 2x2
Subject to: 3x1 + x2 ≥ 27
x1 + x2 ≥ 21
and x1 ≥ 0, x2 ≥ 0
[Ans. Min Z = 48, x1 = 3, x2 =18]
2. Min Z = x1 + x2 + 3x3
Subject to: 3x1 + 2x2 + x3 ≤ 3
2x1 + x2+ 2x3 ≥ 3
and x1 ≥ 0, x2 ≥ 0, x3 ≥ 0
[Ans. Min Z = 3, x1 = 3/4, x2 =0, x3 = 3/4]
3. Max Z = 3x1 - x2
Subject to
2x1 + x2 ≥ 2
x1 + 3x2 ≤ 2
x2≤ 4
and x1 ≥ 0, x2 ≥ 0
[Ans. Max Z = 6, x1 = 2, x2 =0]
4. Max Z = 5x1 - 2x2 +3x3
Subject to: 2x1 + 2x2 - x3 ≥ 2
3x1 - 4x2 ≤ 3
x2 + 3x3 ≤ 5
and x1 ≥ 0, x2 ≥ 0, x3 ≥ 0
[Ans. Max Z = 85/3, x1 = 23/3, x2 =5, x3 =0]
77
Duality in LPP
Every LPP called the primal is associated with another LPP
called dual. Either of the problems is primal with the other one
as dual. The optimal solution of either problem reveals the
information about the optimal solution of the other.
Let the primal problem be
Max Zx = c1x1 + c2x2 + … +cnxn
Subject to restrictions
a11x1 + a12x2 + … + a1nxn ≤ b1
a21x1 + a22x2 + … + a2nxn ≤ b2
.
.
am1x1 + am2x2 + … + amnxn≤ bn
and
x1 ≥ 0, x2 ≥ 0,…, xn≥ 0
The corresponding dual is defined as
Min Zw = b1w1 + b2w2 + … + bmwm
Subject to restrictions
a11w1 + a21w2 + … + am1wm ≥ c1
a12w1 + a22w2 + … + am2wm ≥ c2
.
.
a1nw1 + a2nw2 + ……….+amnwm ≥ cn
and
w1, w2, …,wm≥ 0
78
Matrix Notation
Primal Dual
Max Zx = CX Min Zy = bTY
Subject to Subject to
AX ≤ b , X ≥ 0 AT Y ≥ CT , Y ≥ 0
80
o All constraints have ‘≤’ sign, where the objective
function is of maximization form.
o All constraints have ‘≥’ sign, where the objective
function is of minimization from.
Rules for Converting any Primal into its Dual
1. Minimize the objective function instead of maximizing it,
and Maximize the objective function instead of
minimizing it.
2. Transpose the rows and columns of the constraint co-
efficient.
3. Transpose the co-efficient (c1,c2,…cn) of the objective
function and the right side constants (b1,b2,…bn)
4. Change the inequalities from ‘≤’ to ‘≥’ sign.
Example Problems
Write the dual of the given problems
Example :
Min Zx= 2x2 + 5x3
Subject to:
x1+x2 ≥ 2
2x1+x2+6x3 ≤ 6
x1- x2 +3x3 = 4
x1, x2 , x3 ≥ 0
Solution
Primal
Min Zx= 2x2 + 5x3
81
Subject to:
x1+ x2 ≥2 y1
-2x1- x2- 6x3 ≥ - 6 y2
-x1-+ x2 - 3x3 ≥ - 4 y3
x1- x2 + 3x3 ≥ 4 y4
x1, x2 , x3 ≥ 0
Dual
Max Zw = 2 y1 - 6 y2 - 4 y3 + 4 y4
Subject to: y1 – 2y2 - y3 +y4 ≤ 0
y1 - y2 + y3 -y4 ≤ 2
-6y2 – 3y3 +3y4 ≤ 5
y1, y2, y3, y4 ≥ 0
Example :
Min Zx= 3x1 - 2x2 + 4x3
Subject to: 3x1+5x2 + 4x3 ≥ 7
6x1+x2+3x3 ≥ 4
7x1 - 2x2 -x3 ≥ 10
x1- 2x2 + 5x3 ≥ 3
4x1 + 7x2 - 2x3 ≥ 2
x1, x2 , x3 ≥ 0
Solution
Primal
Min Zx= 3x1 - 2x2 + 4x3
Subject to:
82
3x1+ 5x2 + 4x3 ≥ 7 y1
6x1+x2+ 3x3 ≥ 4 y2
7x1 - 2x2 - x3 ≥ 10 y3
x1- 2x2 + 5x3 ≥ 3 y4
4x1 + 7x2 - 2x3 ≥ 2 y5
x1, x2 , x3 ≥ 0
Dual
Max Zw = 7y1 +4y2 + 10y3 + 3y4 +2y5
Subject to: 3y1 + 6y2 + 7y3 +y4 + 4y5 ≤ 3
5y1 +y2 – 2y3 – 2y4 –+7y5 ≤ - 2
4y1 + 3y2 - y3 +5y4 – 2y5 ≤ - 4
y1, y2, y3, y4, y5 ≥ 0
Example :
Max Z= 2x1+ 3x2 + x3
Subject to:
4x1+ 3x2 + x3 = 6
x1+ 2x2 + 5x3 = 4
x1, x2 ≥ 0
Solution
Primal
Max Zx= 2x1+ 3x2 + x3
Subject to :
83
4x1+ 3x2 + x3 ≤ 6 y1
-4x1 - 3x2 - x3 ≤ -6 y2
x1+ 2x2 + 5x3 ≤ 4 y3
-x1- 2x2 - 5x3 ≤ -4 y4
x1, x2 ,x3≥ 0
Dual
Min Zw = 6y1 – 6y2 + 4y3 –4y4
Subject to
4y1 – 4y2 + y3 –y 4 ≥ 2
3y1 – 3y2 + 2y3 – 2y4 ≥ 3
y1 - y2 + 5y3 – 5y4 ≥ 1
y1, y2, y3, y4≥ 0
Example ;
Min Zx= x1+ x2 + x3
Subject to: x1 - 3x2 + 4x3 = 5
x1- 2x2 ≤ 3
2x2 - x3 ≥ 4
x1, x2 ≥ 0 ,x3 is unrestricted in sign
Solution
Primal
Let (̿̿̿ ̅̅̅)
84
X1- 3X2 + 4(̿̿̿ ̅̅̅) ≥ 5 Y1
-X1+ 3X2 - 4(̿̿̿ ̅̅̅) ≥ -5 Y2
-X1+ 2X2 ≥ -3 Y3
2X2 - (̿̿̿ ̅̅̅) ≥ 4 Y4
X1, X2 , ̿̿̿ ̅̅̅
Dual
Max Zw = 5Y1 – 5Y2 - 3Y3 + 4Y4
Subject to
Y1 - Y2 - Y3 ≤1
-3Y1 + 3Y2 + 2Y3 – 2Y4 ≤ 1
4Y1 – 4Y2 -Y4 ≤ 1
-4Y1 + 4Y2 +Y4 ≤ -1
Y1, Y2, Y3, Y4, ≥ 0
Primal –Dual Relationship
Weak duality property
If x is any feasible solution to the primal problem and w is any
feasible solution to the dual problem then CX ≤ bT Y.
i.e. ZX ≤ Zy
Strong duality property
If x* is an optimal solution for the primal problem and y* is the
optimal solution for the dual problem then CX* = bT Y*
i.e. ZX = Zy
85
Complementary optimal solutions property
At the final iteration, the simplex method simultaneously
identifies an optimal solution x* for primal problem and a
complementary optimal solution w* for the dual problem where
ZX = Zy.
Symmetry property
For any primal problem and its dual problem, all relationships
between them must be symmetric because dual of dual is primal.
Complementary basic solutions property
Each basic solution in the primal problem has a complementary
basic solution in the dual problem where ZX = Zy.
Complementary slackness property
The variables in the primal basic solution and the
complementary dual basic solution satisfy the complementary
slackness relationship as shown in the table.
Primal variable Associated dual variable
Decision variable (xj) Cj– Zj (surplus variable) j = 1, 2, ..n
Slack variable (Si) Yi (decision variable) i = 1, 2, .. n
Duality and Simplex Method
Example ; If you have the following LPM:
Subject to :
3X1 + x2 ≥ 3
4x1 + 3x2 ≥ 6
86
x1+ 2x2 ≤ 3
x 1 , x2 ≥ 0
Solution :
Subject to :
3X1 + x2 ≥ 3 (y1)
x 1 , x2 ≥ 0
87
Then, the dual problem becomes: find y1 , y2 and y3 by which :
Subject to :
3y1 + 4y2- y3 ≤ 2
y1 + 3y2 – 2y3 ≤ 1
y1 , y2, y3 ≥ 0
Subject to :
3y1 + 4y2- y3 + S1 = 2
y1 + 3y2 –y3 + S2 = 1
( ) * +
[ ]
88
1st tableau
Coef.Z Cj 3 6 -3 0 0 Min
Ratio.
(B.V) CB XB
0 2 3 4 -1 1 0 ⁄
0 1 -2 ⁄
1 3 0 1
Zj=CBXj 0 0 0 0 0 0
Δj= Cj - Zj 3 6 -3 0 0
Max Z= 0
Coef.Z Cj 3 6 -3 0 0 Min
Ratio.
(B.V) CB XB
0 ⁄ 0 ⁄ 1 ⁄ ⁄
⁄
6 ⁄ ⁄ 1 ⁄ ⁄
0 1
Zj=CBXj 2 1 0 1 0 -2
Δj= Cj - Zj 0 0 -2
1 1
Max Z= 2
89
Coef.Z Cj 3 6 -3 0 0 Min
Ratio.
(B.V) CB XB
3 ⁄ 1 0 1 ⁄ ⁄
6 ⁄ 0 1 -1 ⁄ ⁄
Zj=CBXj ⁄ 3 6 -3 ⁄ ⁄
Δj= Cj - Zj 0 0 0
⁄ ⁄
Max Z=
3-From the optimal solution tableau for the dual problem, we can get
the optimal solution tableau for the primal problem as follows :
Since the optimal values for the decision variable to the primal problem
can be derived from the optimality condition row (Δj) to the dual
problem, specifically under the slack variables for the dual problem, i.e.,
X*1 = the coefficient of the 1st slack variable in the dual problem
X*2 = the coefficient of the 2nd slack variable in the Δj row = 6/5
90
And , the optimal values for the slack variables in the primal problem are
the same coefficient for the decision variables in the dual problem, i.e.,
= the coefficient of the 1st decision variable in the Δj for the dual
problem.
Coef.Z Cj 2 1 0 0 0 Min
Ratio.
(B.V) CB XB
2 ⁄ 1 0
0
1 ⁄ 0 1
0
0 0 0 0 1
Zj=CBXj 2 1 ⁄ ⁄
0
Δj= Cj - Zj 0 0
⁄ ⁄
Min Z= 0
91
4-Determining the shadow prices for the different resources
(constraints) to the primal problem by using two different
methods:
The 1st shadow price = the coefficient of the 1st slack variable in
the optimality row condition Δj
the 2nd shadow price = the coefficient of the 2nd slack variable
Secondly : One can determine the shadow prices for the set of
resources (constraints) to the primal problem form the optimal
solution for the dual problem from the optimal value for the
decision variable (solution) for it. Therefore;
The 1st shadow price for the 1st resource = the optimum value
for the 1st decision variable. i.e.,
92
The 3rd shadow price = y*3 = 0
.Exercise
1. Max f(x) = 5x1+ 4x2
Subject to:
6 x1 + 4x2 24
x1 + 2x2 6
- x 1 + x2 1
x2 2
x1 ≥ 0, x2 ≥ 0
93
Chapter Three
Transportation Problem
Introduction
The transportation problem is to transport various
amounts of a single homogeneous commodity that are
initially stored at various origins, to different
destinations in such a way that the total transportation
cost is a minimum.
It can also be defined as to ship goods from various
origins to various destinations in such a manner that the
transportation cost is a minimum.
The availability as well as the requirements is finite. It
is assumed that the cost of shipping is linear.
Mathematical Formulation
Let there be m origins, ith origin possessing ai units of a
certain product
Let there be n destinations, with destination j requiring
bj units of a certain product
97
Let cij be the cost of shipping one unit from ith source to
jth destination
Let xij be the amount to be shipped from ith source to jth
destination
It is assumed that the total availabilities Σai satisfy the
total requirements Σbji.e.
Σai = Σbj(i = 1, 2, 3 … m and j = 1, 2, 3 …n)
The problem now, is to determine non-negative of xij
satisfying both the availability constraints
Tabular Representation
Let ‘m’ denote number of factories (F1, F2 … Fm)
Let ‘n’ denote number of warehouse (W1, W2 … Wn)
98
W→
F Capacities
W1 W2 … Wn
↓ (Availability)
F1 c11 c12 … c1n a1
F2 c21 c22 … c2n a2
. . . . . .
. . . . . .
Fm cm1 cm2 … cmn am
Required b1 b2 … bn Σai = Σbj
W→
F Capacities
W1 W2 … Wn
↓ (Availability)
F1 x11 x12 … x1n a1
F2 x21 x22 … x2n a2
. . . . . .
. . . . . .
Fm xm1 xm2 … xmn am
Required b1 b2 … bn Σai = Σbj
102
Solution
W1 W2 W3 W5 Availability
5 2
F1 7 2 0
(19) (30)
6 3
F2 9 3 0
(30) (40)
4 14
F3 18 14 0
(70) (20)
5 8 7 14
Requirement 0 6 4 0
0 0
Problem (2):
103
D1 D2 D3 D4 Supply
O1 1 5 3 3 34
O2 3 3 1 2 15
O3 0 2 2 3 12
O4 2 7 2 4 19
Demand 21 25 17 17 80
Solution
D1 D2 D3 D4 Supply
21 13
O1 (1) (5) 34 13 0
12 3
O2 (3) (1) 15 3 0
12
O3 (2) 12 0
2 17
O4 (2) (4) 19 17
Demand 21 25 17 17
0 12 14 0
0 2
0
Initial Basic Feasible Solution
104
x11 = 21, x12 = 13, x22 = 12, x23 = 3, x33 = 12, x43 = 2, x44 =
17
The transportation cost is 21 (1) + 13 (5) + 12 (3) + 3 (1) +
12 (2) + 2 (2) + 17 (4) = 221
Problem (3):
From To Supply
2 11 10 3 7 4
1 4 7 2 1 8
3 1 4 8 12 9
Demand 3 3 4 5 6
Solution
From To Supply
3 1
(2) (11) 4 1 0
2 4 2
8 6 2
(4) (7) (2)
0
3 6
(8) (12) 9 6 0
3 3 4 5 6
Demand 0 2 0 3 0
0 0
105
x11 = 3, x12 = 1, x22 = 2, x23 = 4, x24 = 2, x34 = 3, x35 = 6
The transportation cost is 3 (2) + 1 (11) + 2 (4) + 4 (7) + 2
(2) + 3 (8) + 6 (12) = 153
106
new minimum cost in the first row. Cross out the
first row and move to second row.
Step 3
Repeat steps 1 and 2 for the reduced transportation table
until all the requirements are satisfied
Example:
Determine the initial basic feasible solution using Row
Minima Method
W1 W2 W3 W4 Availability
F1 19 30 50 10 7
F2 70 30 40 60 9
F3 40 80 70 20 18
Requirement 5 8 7 14
Solution
W1 W2 W3 W4
7
X
F1 (19) (30) (50) (10)
F2 9
(70) (30) (40) (60)
F3 18
(40) (80) (70) (20)
5 8 7 7
107
W1 W2 W3 W4
7
X
F1 (19) (30) (50) (10)
8
F2 (70) (30) (40) (60) 1
F3 18
(40) (80) (70) (20)
5 X 7 7
W1 W2 W3 W4
7
X
F1 (19) (30) (50) (10)
8 1
X
F2 (70) (30) (40) (60)
109
W1 W2 W3 W4 Availability
F1 19 30 50 10 7
F2 70 30 40 60 9
F3 40 8 70 20 18
Requirement 5 8 7 14
Solution
W1 W2 W3 W4
F1 7
(19) (30) (50) (10)
F2 9
(70) (30) (40) (60)
F3 8 10
(40) (8) (70) (20)
5 X 7 14
W1 W2 W3 W4
7
F1 X
(19) (30) (50) (10)
F2 9
(70) (30) (40) (60)
8
F3 10
(40) (8) (70) (20)
5 X 7 7
W1 W2 W3 W4
110
7
F1 X
(19) (30) (50) (10)
F2 9
(70) (30) (40) (60)
8 7
F3 3
(40) (8) (70) (20)
5 X 7 X
W1 W2 W3 W4
7
F1 X
(19) (30) (50) (10)
F2 9
(70) (30) (40) (60)
3 8 7
F3 X
(40) (8) (70) (20)
2 X 7 X
W1 W2 W3 W4
7
F1 X
(19) (30) (50) (10)
2 7
F2 X
(70) (30) (40) (60)
3 8 7
F3 X
(40) (8) (70) (20)
X X X X
111
Initial Basic Feasible Solution
Step1: For each row of the table, identify the smallest and
the next to smallest cost. Determine the difference
between them for each row. These are called penalties. Put
them aside by enclosing them in the parenthesis against the
respective rows. Similarly compute penalties for each
column.
113
Example:
Problem (1):
W1 W2 W3 W4 Availability
F1 19 30 50 10 7
F2 70 30 40 60 9
F3 40 8 70 20 18
Requirement 5 8 7 14
Solution
W1 W2 W3 W4 Av. Penalty
F1 19 30 50 10 7 19-10=9
F2 70 30 40 60 9 40-30=10
F3 40 8 70 20 18 20-8=12
Req. 5 8 7 14
Penalty 40-19=21 30-8=22 50-40=10 20-10=10
W1 W2 W3 W4 Availability Penalty
F1 (19) (30) (50) (10) 7 9
F2 (70) (30) (40) (60) 9 10
114
F3 (40) 8(8) (70) (20) 18/10 12
Requirement 5 8/0 7 14
Penalty 21 22 10 10
W1 W2 W3 W4 Availability Penalty
F1 5(19) (30) (50) (10) 7/2 9
F2 (70) (30) (40) (60) 9 20
F3 (40) 8(8) (70) (20) 18/10 20
Requirement 5/0 X 7 14
Penalty 21 X 10 10
W1 W2 W3 W4 Availability Penalty
F1 5(19) (30) (50) (10) 7/2 40
F2 (70) (30) (40) (60) 9 20
F3 (40) 8(8) (70) 10(20) 18/10/0 50
Requirement X X 7 14/4
Penalty X X 10 10
W1 W2 W3 W4 Availability Penalty
F1 5(19) (30) (50) 2(10) 7/2/0 40
F2 (70) (30) (40) (60) 9 20
F3 (40) 8(8) (70) 10(20) X X
115
Requirement X X 7 14/4/2
Penalty X X 10 50
W1 W2 W3 W4 Availability Penalty
F1 5(19) (30) (50) 2(10) X X
F2 (70) (30) 7(40) 2(60) X X
F3 (40) 8(8) (70) 10(20) X X
Requirement X X X X
Penalty X X X X
Solution
116
Stores Availability Penalty
I II III IV
A (21) (16) (15) (13) 11 2
Warehouse B (17) (18) (14) (23) 13 3
C (32) (27) (18) (41) 19 9
Requirement 6 10 12 15
Penalty 4 2 1 10
I II III IV
Penalty 15 9 4 18
I II III IV
117
B 6(17) (18) (14) 4(23) 13/9/3 3
C (32) (27) (18) (41) 19 9
Req. 6/0 10 12 X
Penalty 15 9 4 X
118
The transportation cost is 11 (13) + 6 (17) + 3 (18) + 4 (23)
+ 7 (27) + 12 (18) = 796
Exercise
1. Determine an initial basic feasible solution to the
following transportation problem using north-west corner
rule.
I II III IV Supply
A 13 11 15 20 2000
From B 17 14 12 13 6000
C 18 18 15 12 7000
Demand 3000 3000 4000 5000
119
D1 D2 D3 D4 D5 Capacity
O1 2 11 10 3 7 4
From O2 1 4 7 2 1 8
O3 3 9 4 8 12 9
Demand 3 3 4 5 6 21
[Ans. Min cost = 68]
4. Determine the minimum cost to the following
transportation problem using matrix minima method and
vogel’s method
D1 D2 D3 D4 Capacity
O1 1 2 1 4 30
From O2 3 3 2 1 50
O3 4 2 5 9 20
Demand 20 40 30 10
[Ans. Min cost = 180]
Examining the Initial Basic Feasible Solution for
Non-Degeneracy
Examine the initial basic feasible solution for non-
degeneracy. If it is said to be non-degenerate then it has
the following two properties
Initial basic feasible solution must contain
exactly m + n – 1 number of individual
allocations.
120
These allocations must be in independent
positions
Independent Positions
Non-Independent Positions
121
Transportation Algorithm for Minimization
Problem (MODI Method)
Step1: Construct the transportation table entering the origin
capacities ai, the destination requirement bj and the cost cij
Step 2: Find an initial basic feasible solution by vogel’s
method or by any of the given method.
Step 3: For all the basic variables xij, solve the system of
equations ui + vj = cij, for all i, j for which cell (i, j) is in the
basis, starting initially with some ui = 0, calculate the
values of ui and vj on the transportation table
Step 4: Compute the cost differences dij = cij – ( ui + vj ) for
all the non-basic cells
Step 5: Apply optimality test by examining the sign of each
dij
If all dij ≥ 0, the current basic feasible solution is
optimal
If at least one dij< 0, select the variable x rs (most
negative) to enter the basis.
Solution under test is not optimal if any dij is
negative and further improvement is required by
repeating the above process.
122
Step 6: Let the variable xrs enter the basis. Allocate an
unknown quantity Ө to the cell (r, s). Then construct a loop
that starts and ends at the cell (r, s) and connects some of
the basic cells. The amount Ө is added to and subtracted
from the transition cells of the loop in such a manner that
the availabilities and requirements remain satisfied.
Step 7: Assign the largest possible value to the Ө in such a
way that the value of at least one basic variable becomes
zero and the other basic variables remain non-negative. The
basic cell whose allocation has been made zero will leave
the basis.
Step 8: Now, return to step 3 and repeat the process until
an optimal solution is obtained.
Worked Examples
Example 1
Find an optimal solution
W1 W2 W3 W4 Availability
F1 19 30 50 10 7
F2 70 30 40 60 9
F3 40 8 70 20 18
Requirement 5 8 7 14
Solution
123
1. Applying vogel’s approximation method for finding the
initial basic feasible solution
W1 W2 W3 W4 Availability Penalty
F1 5(19) (30) (50) 2(10) X X
F2 (70) (30) 7(40) 2(60) X X
F3 (40) 8(8) (70) 10(20) X X
Requirement X X X X
Penalty X X X X
ui
(19) (10) u1= -10
(40) (60) u2 = 40
(8) (20) u3 = 0
vj v1= 29 v2 = 8 v3 = 0 v4 = 20
Assign a ‘u’ value to zero. (Convenient rule is to select the
ui, which has the largest number of allocations in its row)
124
Let u3 = 0, then
u3 + v4= 20 which implies 0 + v 4 = 20, so v4 = 20
u2 + v4= 60 which implies u2 + 20 = 60, so u2 = 40
u1 + v4= 10 which implies u1 + 20 = 10, so u1 = -10
u2 + v3= 40 which implies 40 + v3 = 40, so v3 = 0
u3 + v2= 8 which implies 0 + v2 = 8, so v2 = 8
u1 + v1= 19 which implies -10 + v1= 19, so v1 = 29
dij = cij – ( ui + vj )
32 60
1 -18
11 70
5.Optimality test
dij< 0 i.e. d22 = -18
so x22 is entering the basis
125
6.Construction of loop and allocation of unknown
quantity Ө
ui
(19) (10) u1= -10
(30) (40) u2 = 22
(8) (20) u3 = 0
vj v1= 29 v2 = 8 v3 = 18 v4 = 20
126
cij ui + vj
(30) (50) -2 8
(70) (60) 51 42
(40) (70) 29 18
dij = cij – ( ui + vj )
32 42
19 18
11 52
Since dij> 0, an optimal solution is obtained with minimal
cost 743
Example 2
Solve by lowest cost entry method and obtain an optimal
solution for the following problem
Available
50 30 220 1
From 90 45 170 3
250 200 50 4
Required 4 2 2
Solution
By lowest cost entry method
127
Available
1(30) 1/0
From 2(90) 1(45) 3/2/0
2(250) 2(50) 4/2/0
Required 4/2/2 2/1/0 2/0
ui
(30) u1= -15
(90) (45) u2 = 0
(250) (50) u3 = 160
vj v1= 90 v2 = 45 v3 = -110
128
Calculation of cost differences for non-basic cellsdij = cij – ( ui + vj
)
cij u i + vj
50 220 75 -125
170 -110
200 205
dij = cij – ( ui + vj )
-25 345
280
-5
Optimality test
dij< 0 i.e. d11 = -25 is most negative
So x11 is entering the basis
Construction of loop and allocation of unknown
quantity Ө
1(90) 2(45)
2(250) 2(50)
ui
(50) u1= -40
(90) (45) u2 = 0
(250) (50) u3 = 160
vj v1= 90 v2 = 45 v3 = -110
cij u i + vj
30 220 5 -150
170 -110
200 205
dij = cij – ( ui + vj )
25 370
130
280
-5
Optimality test
dij< 0 i.e. d32 = -5
So x32 is entering the basis
Construction of loop and allocation of unknown
quantity Ө
1(50)
3(90) 0(45)
2(200) 2(50)
131
III Iteration
Calculation of ui and vjby equationui + vj = cij
ui
(50) u1= -40
(90) (45) u2 = 0
(200) (50) u3 = 155
vj v1= 90 v2 = 45 v3 = -105
u i + vj
cij
30 220 5 -145
170 -105
250 245
dij = cij – ( ui + vj )
25 365
275
5
Example 3
132
Is x13 = 50, x14 = 20, x21 = 55, x31 = 30, x32 = 35, x34 = 25 an
optimal solution to the transportation problem.
Available
6 1 9 3 70
From 11 5 2 8 55
10 12 4 7 90
Required 85 35 50 45
Solution
Available
50(9) 20(3) X
From 55(11) X
30(10) 35(12) 25(7) X
Required X X X X
133
ui
(9) (3) u1= -4
(11) u2 = 1
(10) (12) (7) u3 = 0
vj v1= 10 v2 = 12 v3 = 13 v4 = 7
Calculation of cost differences for non-basic cells
dij = cij – ( ui + vj )
cij ui + vj
6 1 6 8
5 2 8 13 14 8
4 13
dij = cij – ( ui + vj )
0 -7
-8 -12 0
-9
Optimality test
dij< 0 i.e. d23 = -12 is most negative
So x23 is entering the basis
Construction of loop and allocation of unknown
quantity Ө
134
min(50-Ө, 55-Ө, 25-Ө) = 25 which gives Ө = 25. Therefore
x34 is outgoing as it becomes zero.
25(9) 45(3)
30(11) 25(2)
55(10) 35(12)
Minimum transportation cost is 25 (9) + 45 (3) + 30 (11) +
25 (2) + 55 (10) + 35 (12) = Rs. 1710
II iteration
Calculation of ui and vjby equationui + vj = cij
ui
(9) (3) u1= 8
(11) (2) u2 = 1
(10) (12) u3 = 0
vj v1= 10 v2 = 12 v3 = 1 v4 = -5
135
cij ui + vj
6 1 20
18
-
5 8 13
4
-
4 7 1
5
dij = cij – ( ui + vj )
-12 -19
-8 12
3 12
Optimality test
dij< 0 i.e. d12 = -19 is most negative
So x12 is entering the basis
Construction of loop and allocation of unknown
quantity Ө
136
25(1) 45(3)
5(11) 50(2)
80(10) 10(12)
Minimum transportation cost is 25 (1) + 45 (3) + 5 (11) +
50 (2) + 80 (10) + 10 (12) = 1235
III Iteration
Calculation of ui and vjby equationui + vj = cij
ui
(1) (3) u1= -11
(11) (2) u2 = 1
(10) (12) u3 = 0
vj v1= 10 v2 = 12 v3 = 1 v4 = 14
5 8
13 15
4 7 1
14
137
dij = cij – ( ui + vj )
7 19
-8 -7
3 -7
Optimality test
dij< 0 i.e. d22 = -8 is most negative
So x22 is entering the basis
Construction of loop and allocation of unknown
quantity Ө
25(1) 45(3)
5(5) 50(2)
85(10) 5(12)
Minimum transportation cost is 25 (1) + 45 (3) + 5 (5) + 50
(2) + 85 (10) + 5 (12) = Rs. 1195
IV Iteration
Calculation of ui and vj : - ui + vj = cij
138
ui
(1) (3) u1= -11
(5) (2) u2 = -7
(10) (12) u3 = 0
vj v1= 10 v2 = 12 v3 = 9 v4 = 14
Calculation of cost differences for non-basic cells
dij = cij – ( ui + vj )
cij ui + vj
-
6 9 -2
1
8 3 7
11
4 7 9
14
dij = cij – ( ui + vj )
7 11
8 1
-5 -7
Optimality test
dij< 0 i.e. d34 = -7 is most negative
So x34 is entering the basis
Construction of loop and allocation of unknown
quantity Ө
139
min(5-Ө, 45-Ө) = 5 which gives Ө = 5. Therefore x 32 is
outgoing as it becomes zero.
30(1) 40(3)
5(5) 50(2)
85(10) 5(7)
Minimum transportation cost is 30 (1) + 40 (3) + 5 (5) + 50
(2) + 85 (10) + 5 (7) = Rs. 1160
V Iteration
Calculation of ui and vjby equationui + vj = cij
ui
(1) (3) u1= -4
(5) (2) u2 = 0
(10) (7) u3 = 0
vj v1= 10 v2 = 5 v3 = 2 v4 = 7
140
cij ui + vj
6 9 6 -2
11 8 10 7
12 4 5 2
dij = cij – ( ui + vj )
0 11
1 1
7 2
Since dij> 0, an optimal solution is obtained with minimal
cost Rs.1160. Furthermore d11 = 0 which indicates that
alternative optimal solution also exists.
Exercises
To Supply
2 3 11 7 6
1 0 6 1 1
From
5 8 15 10 10
Demand 7 5 3 2 17
141
[Ans. x12 = 5, x13 = 1, x24 = 1, x31 = 7, x33 = 2, x34 = 1 Min
cost = Rs 102]
2. Using North-West Corner rule for initial basic
feasible solution, obtain an optimum basic feasible
solution to the following problem
To Available
7 3 4 2
From 2 1 3 3
3 4 6 5
Demand 4 1 5 10
To Supply
10 7 3 6 3
1 6 7 3 5
From
7 4 5 3 7
Demand 3 2 6 4
[Ans. x13 = 3, x21 = 3, x24 = 2, x32 = 2, x33 = 3, x34 = 2, Min
cost = 47]
142
Chapter Four
Network Analysis
Shortest Route Problem
The criterion of this method is to find the shortest
distance between two nodes with minimal cost.
Example 1
Find the shortest path
Solution
143
Solved
nodes
Closest
Last connection
nth nearest node
Min. distance
directly Total
connected
n connected distance
unsolved
to involved
node
unsolved
nodes
1 a c 7 c 7 a-c
a b 13 b 13 a-b
2
c e 7+6 =13 e 13 c-e
b d 13+5 =18 d 18 b-d
3 c f 7+11 =18 f 18 c-f
e h 13+8 =21 - - -
e h 13+8 =21 h 21 e-h
4 d g 18+9 =27 - - -
f h 18+5 =23 - - -
e g 13+10 =23 g 23 e-g
5 h i 21+10 =31 - - -
d g 18+9 =27 - - -
g i 23+6 =29 i 29 g-i
6
h i 21+10 =31 - - -
144
The shortest path from a to i is a → c →e →g → i
Distance = 7 + 6 + 10 + 6 = 29 units
Example 2
145
Solution:
Solved
nodes
Closest
directly Total nth
connected Min. Last
n connected distance nearest
unsolved distance connection
to involved node
node
unsolved
nodes
1 1 3 1 3 1 1-3
1 2 5 - - -
2
3 2 1+2 =3 2 3 3-2
2 5 3+1 =4 5 4 2-5
3
3 4 1+6 =7 - - -
2 6 3+6 =9 - - -
4 3 4 1+6 =7 4 7 3-4
5 4 4+3 =7 4 7 5-4
2 6 3+6 =9 6 9 2-6
5 4 6 7+4 =11 - - -
5 6 4+5 =9 6 9 5-6
4 7 7+6 =13 - - -
6 5 7 4+9 =13 - - -
6 7 9+2 =11 7 11 6-7
146
1 →3 → 2 → 6 →7
Total distance is 11 units
1 → 3 → 2 →5 → 6 →7
Total distance = 11 units
Minimal Spanning Tree Problem
A tree is defined to be an undirected, acyclic and
connected graph. A spanning tree is a sub graph of G
(undirected, connected graph), is a tree and contains all
147
the vertices of G. A minimum spanning tree is a
spanning tree but has weights or lengths associated
with edges and the total weight is at the minimum.
Prim’s Algorithm
It starts at any vertex (say A) in a graph and finds
the least cost vertex (say B) connected to the start
vertex.
Now either from A or B, it will find the next least
costly vertex connection, without creating cycle
(say C)
Now either from A, B or C find the next least
costly vertex connection, without creating a cycle
and so on.
Eventually all the vertices will be connected
without any cycles and a minimum spanning tree
will be the result.
Example 1
Suppose it is desired to establish a cable
communication network that links major cities, which
is shown in the figure. Determine how the cities are
148
connected such that the total cable mileage is
minimized.
Solution
C = {LA} C' = {SE, DE, DA, EH, NY, DC}
C = {LA, SE} C' = {DE, DA, EH, NY, DC}
C = {LA, SE, DE} C' = {DA, EH, NY, DC}
C = {LA, SE, DE, DA} C' = {EH, NY, DC}
C = {LA, SE, DE, DA, EH} C' = {NY, DC}
C = {LA, SE, DE, DA, EH, NY} C' = {DC}
C = {LA, SE, DE, DA, EH, NY, DC} C' = { }
149
Thus the total cable mileage is 1100 + 1300 + 780 +
900 + 800 + 200 = 5080
Example 2
For the following graph obtain the minimum spanning
tree. The numbers on the branches represent the cost.
Solution
C = {A} C' = {B, C, D, E, F, G}
C = {A, D} C' = {B, C, E, F, G}
C = {A, D, B} C' = {C, E, F, G}
150
C = {A, D, B, C} C' = {E, F, G}
C = {A, D, B, C, G} C' = {E, F}
C = {A, D, B, C, G, F} C' = {E}
C = {A, D, B, C, G, F, E} C' = { }
The resultant network is
Cost = 2 + 1 + 4 + 3 + 3 + 5 = 18 units
Example 3
Solve the minimum spanning problem for the given
network. The numbers on the branches represent in
terms of miles.
151
Solution
C = {1} C' = {2, 3, 4, 5, 6}
C = {1, 2} C' = {3, 4, 5, 6}
C = {1, 2, 5} C' = {3, 4, 6}
C = {1, 2, 5, 4} C' = {3, 6}
C = {1, 2, 5, 4, 6} C' = {3}
C = {1, 2, 5, 4, 6, 3} C' = {}
The resultant network is
1 + 4 + 5+ 3 + 3 = 16 miles
152
Maximal Flow Problem
Algorithm
Step1: Find a path from source to sink that can
accommodate a positive flow of material. If no path
exists go to step 5
Step2: Determine the maximum flow that can be
shipped from this path and denote by ‘k’ units.
Step3:Decrease the direct capacity (the capacity in the
direction of flow of k units) of each branch of this path
‘k’ and increase the reverse capacity k1. Add ‘k’ units
to the amount delivered to sink.
Step4: Go to step1
Step5: The maximal flow is the amount of material
delivered to the sink. The optimal shipping schedule is
determined by comparing the original network with the
final network. Any reduction in capacity signifies
shipment.
Example 1
Consider the following network and determine the
amount of flow between the networks.
153
Solution
Iteration 1: 1 – 3 – 5
Iteration 2: 1 – 2 – 3 – 4 – 5
154
Iteration 3: 1 – 4 – 5
Iteration 4: 1 – 2 – 5
155
Iteration 5: 1 – 3 – 2 – 5
156
Example 2
Solve the maximal flow problem
Solution
Iteration 1: O – A – D – T
157
Iteration 2: O – B – E – T
Iteration 3: O – A – B – D – T
158
Iteration 4: O – C – E – D – T
Iteration 5: O – C – E – T
159
Iteration 6: O – B – D – T
160
Exercise
1. Find the shortest path
161
Introduction to CPM / PERT Techniques
CPM/PERT or Network Analysis as the technique is
sometimes called, developed along two parallel
streams, one industrial and the other military.
CPM (Critical Path Method)was the discovery of
M.R.Walker of E.I.Du Pont de Nemours & Co. and
J.E.Kelly of Remington Rand, circa 1957. The
computation was designed for the UNIVAC-I
computer. The first test was made in 1958, when CPM
was applied to the construction of a new chemical
plant. In March 1959, the method was applied to
maintenance shut-down at the Du Pont works in
Louisville, Kentucky. Unproductive time was reduced
from 125 to 93 hours.
PERT (Project Evaluation and Review
Technique)was devised in 1958 for the POLARIS
missile program by the Program Evaluation Branch of
the Special Projects office of the U.S.Navy, helped by
the Lockheed Missile Systems division and the
Consultant firm of Booz-Allen & Hamilton. The
calculations were so arranged so that they could be
carried out on the IBM Naval Ordinance Research
Computer (NORC) at Dahlgren, Virginia.
162
The methods are essentially network-oriented
techniques using the same principle. PERT and CPM
are basically time-oriented methods in the sense that
they both lead to determination of a time schedule for
the project. The significant difference between two
approaches is that the time estimates for the different
activities in CPM were assumed to be deterministic
while in PERT these are described probabilistically.
These techniques are referred as project scheduling
techniques.
In CPM activities are shown as a network of
precedence relationships using activity-on-node
network construction
– Single estimate of activity time
– Deterministic activity times
USED IN: Production management - for the jobs of
repetitive in nature where the activity time estimates
can be predicted with considerable certainty due to the
existence of past experience.
In PERT activities are shown as a network of
precedence relationships using activity-on-arrow
network construction
– Multiple time estimates
– Probabilistic activity times
USED IN: Project management - for non-repetitive
jobs (research and development work), where the time
and cost estimates tend to be quite uncertain. This
technique uses probabilistic time estimates.
Benefits of PERT/CPM
Useful at many stages of project management
163
Mathematically simple
Give critical path and slack time
Provide project documentation
Useful in monitoring costs
Limitations of PERT/CPM
Clearly defined, independent and stable activities
Specified precedence relationships
Over emphasis on critical paths
165
Resource allocation usually incurs a compromise
and the choice of this compromise depends on the
judgment of managers.
4. Controlling
The final phase in project management is
controlling. Critical path methods facilitate the
application of the principle of management by
expectation to identify areas that are critical to
the completion of the project.
By having progress reports from time to time and
updating the network continuously, a better
financial as well as technical control over the
project is exercised.
Arrow diagrams and time charts are used for
making periodic progress reports. If required, a
new course of action is determined for the
remaining portion of the project.
The Framework for PERT and CPM
Essentially, there are six steps which are common to
both the techniques. The procedure is listed below:
I. Define the Project and all of its significant
activities or tasks. The Project (made up of
several tasks) should have only a single start
activity and a single finish activity.
II. Develop the relationships among the activities.
Decide which activities must precede and which
must follow others.
III. Draw the "Network" connecting all the activities.
Each Activity should have unique event numbers.
Dummy arrows are used where required to avoid
giving the same numbering to two activities.
166
IV. Assign time and/or cost estimates to each activity
V. Compute the longest time path through the
network. This is called the critical path.
VI. Use the Network to help plan, schedule, and
monitor and control the project.
The Key Concept used by CPM/PERT is that a small
set of activities, which make up the longest path
through the activity network control the entire project.
If these "critical" activities could be identified and
assigned to responsible persons, management resources
could be optimally used by concentrating on the few
activities which determine the fate of the entire project.
Non-critical activities can be replanned, rescheduled
and resources for them can be reallocated flexibly,
without affecting the whole project.
Five useful questions to ask when preparing an activity
network are:
Is this a Start Activity?
Is this a Finish Activity?
What Activity Precedes this?
What Activity Follows this?
What Activity is Concurrent with this?
Network Diagram Representation
In a network representation of a project certain
definitions are used
1. Activity
Any individual operation which utilizes resources and
has an end and a beginning is called activity. An arrow
is commonly used to represent an activity with its head
indicating the direction of progress in the project. These
are classified into four categories
167
1. Predecessor activity – Activities that must be
completed immediately prior to the start of
another activity are called predecessor activities.
2. Successor activity – Activities that cannot be
started until one or more of other activities are
completed but immediately succeed them are
called successor activities.
3. Concurrent activity – Activities which can be
accomplished concurrently are known as
concurrent activities. It may be noted that an
activity can be a predecessor or a successor to an
event or it may be concurrent with one or more of
other activities.
4. Dummy activity – An activity which does not
consume any kind of resource but merely depicts
the technological dependence is called a dummy
activity.
168
2. Event
An event represents a point in time signifying the
completion of some activities and the beginning of new
ones. This is usually represented by a circle in a
network which is also called a node or connector.
The events are classified in to three categories
1. Merge event – When more than one activity
comes and joins an event such an event is known
as merge event.
2. Burst event – When more than one activity
leaves an event such an event is known as burst
event.
3. Merge and Burst event – An activity may be
merge and burst event at the same time as with
respect to some activities it can be a merge event
and with respect to some other activities it may
be a burst event.
3. Sequencing
The first prerequisite in the development of network is
to maintain the precedence relationships. In order to
169
make a network, the following points should be taken
into considerations
What job or jobs precede it?
What job or jobs could run concurrently?
What job or jobs follow it?
What controls the start and finish of a job?
Since all further calculations are based on the network,
it is necessary that a network be drawn with full care.
Rules for Drawing Network Diagram
Rule 1
Each activity is represented by one and only one arrow
in the network
Rule 2
No two activities can be identified by the same end
events
Rule 3
In order to ensure the correct precedence relationship in
the arrow diagram, following questions must be
checked whenever any activity is added to the network
What activity must be completed immediately
before this activity can start?
170
What activities must follow this activity?
What activities must occur simultaneously with
this activity?
In case of large network, it is essential that certain good
habits be practiced to draw an easy to follow network
Try to avoid arrows which cross each other
Use straight arrows
Do not attempt to represent duration of activity
by its arrow length
Use arrows from left to right. Avoid mixing two
directions, vertical and standing arrows may be
used if necessary.
Use dummies freely in rough draft but final
network should not have any redundant dummies.
The network has only one entry point called start
event and one point of emergence called the end
event.
Common Errors in Drawing Networks
The three types of errors are most commonly observed
in drawing network diagrams
1. Dangling
To disconnect an activity before the completion of all
activities in a network diagram is known as dangling.
As shown in the figure activities (5 – 10) and (6 – 7)
are not the last activities in the network. So the diagram
is wrong and indicates the error of dangling
171
2. Looping or Cycling
Looping error is also known as cycling error in a
network diagram. Drawing an endless loop in a
network is known as error of looping as shown in the
following figure.
3. Redundancy
Unnecessarily inserting the dummy activity in network
logic is known as the error of redundancy as shown in
the following diagram
172
Advantages and Disadvantages
PERT/CPM has the following advantages:
A PERT/CPM chart explicitly defines and makes
visible dependencies (precedence relationships)
between the elements,
PERT/CPM facilitates identification of the
critical path and makes this visible,
PERT/CPM facilitates identification of early
start, late start, and slack for each activity,
PERT/CPM provides for potentially reduced
project duration due to better understanding of
dependencies leading to improved overlapping of
activities and tasks where feasible.
174
i. Earliest starting time of activity (i, j) is the
earliest event time of the tail end event i.e. (Es)ij
= Ei
ii. Earliest finish time of activity (i, j) is the earliest
starting time + the activity time i.e. (Ef)ij =
(Es)ij + Dijor (Ef)ij = Ei + Dij
iii. Earliest event time for event j is the maximum of
the earliest finish times of all activities ending in
to that event i.e. Ej = max [(Ef)ij for all
immediate predecessor of (i, j)] or Ej =max [Ei +
Dij]
2. Backward Pass computation (for latest allowable
time)
Step 1: For ending event assume E = L.
Remember that all E’s have been computed by
forward pass computations.
Step 2: Latest finish time for activity (i, j) is
equal to the latest event time of event j i.e. (Lf)ij =
Lj
Step 3: Latest starting time of activity (i, j) = the
latest completion time of (i, j) – the activity time
or (Ls)ij =(Lf)ij - Dij or (Ls)ij = Lj - Dij
175
Step 4: Latest event time for event ‘i’ is the
minimum of the latest start time of all activities
originating from that event i.e. Li = min [(Ls)ij for
all immediate successor of (i, j)] = min [(Lf) ij -
Dij] = min [Lj - Dij]
3. Determination of floats and slack times
There are three kinds of floats
Total float – The amount of time by which the
completion of an activity could be delayed
beyond the earliest expected completion time
without affecting the overall project duration
time.
Mathematically
(Tf)ij = (Latest start – Earliest start) for activity (i–j)
(Tf)ij = (Ls)ij - (Es)ij or (Tf)ij = (Lj - Dij) - Ei
Free float – The time by which the completion of
an activity can be delayed beyond the earliest
finish time without affecting the earliest start of a
subsequent activity.
Mathematically
(Ff)ij = (Earliest time for event j – Earliest time for
event i) – Activity time for ( i, j)
176
(Ff)ij = (Ej - Ei) - Dij
Independent float – The amount of time by
which the start of an activity can be delayed
without effecting the earliest start time of any
immediately following activities, assuming that
the preceding activity has finished at its latest
finish time.
Mathematically
(If)ij = (Ej - Li) - Dij
The negative independent float is always taken as
zero.
Event slack - It is defined as the difference
between the latest event and earliest event times.
Mathematically
Head event slack = Lj – Ej,
Tail event slack = Li - Ei
4. Determination of critical path
Critical event – The events with zero slack times
are called critical events. In other words the event
i is said to be critical if Ei = Li
Critical activity – The activities with zero total
float are known as critical activities. In other
177
words an activity is said to be critical if a delay in
its start will cause a further delay in the
completion date of the entire project.
Critical path – The sequence of critical activities
in a network is called critical path. The critical
path is the longest path in the network from the
starting event to ending event and defines the
minimum time required to complete the project.
Exercise
1. What is PERT and CPM?
178
Worked Examples on CPM
Example 1: Determine the early start and late start in
respect of all node points and identify critical path for
the following network.
179
Network Analysis Table
Earliest Time Latest Time Float
Normal
Activity(i, Finish Start Time
Time Start Finish
j) (Ei + (Li - (Li - Dij )
(Dij) (Ei) (Li)
Dij ) Dij ) - Ei
(1, 2) 10 0 10 0 10 0
(1, 3) 8 0 8 1 9 1
(1, 4) 9 0 9 1 10 1
(2, 5) 8 10 18 10 18 0
(4, 6) 7 9 16 10 17 1
(3, 7) 16 8 24 9 25 1
(5, 7) 7 18 25 18 25 0
(6, 7) 7 16 23 18 25 2
(5, 8) 6 18 24 18 24 0
(6, 9) 5 16 21 17 22 1
(7, 10) 12 25 37 25 37 0
(8, 10) 13 24 37 24 37 0
(9, 10) 15 21 36 22 37 1
From the table, the critical nodes are (1, 2), (2, 5), (5,
7), (5, 8), (7, 10) and (8, 10)
From the table, there are two possible critical paths
i. 1 → 2 → 5 → 8 → 10
ii. 1 → 2 → 5 → 7 → 10
180
Example 2: Find the critical path and calculate the
slack time for the following network
181
From the above table, the critical nodes are the
activities (1, 3), (3, 5) and (5, 9)
182
Construct the network and compute
1. TE and TL for each event
2. Float for each activity
3. Critical path and its duration
Solution
The network is
Event
1 2 3 4 5 6 7 8 9 10
No.:
TE: 0 4 1 5 7 11 15 17 18 25
TL: 0 12 1 13 7 16 15 17 18 25
183
TL (Head
Activity Duration TE (Tail event) Float
event)
(1 – 2) 4 0 12 8
(1 – 3) 1 0 1 0
(2 – 4) 1 4 13 8
(3 – 4) 1 1 13 11
(3 – 5) 6 1 7 0
(4 – 9) 5 5 18 8
(5 – 6) 4 7 16 5
(5 – 7) 8 7 15 0
(6 – 8) 1 11 17 5
(7 – 8) 2 15 17 0
(8 – 9) 1 17 18 0
(8 – 10) 8 17 25 0
(9 – 10) 7 18 25 0
The resultant network shows the critical path
186
Task: A B C D E F G H I J K
Least time: 4 5 8 2 4 6 8 5 3 5 6
Greatest
8 10 12 7 10 15 16 9 7 11 13
time:
Most
likely 5 7 11 3 7 9 12 6 5 8 9
time:
187
Expected Start Finish Total
Task
time (te) Earliest Latest Earliest Latest float
A 5.33 0 0 5.33 5.33 0
B 7.17 0 8.83 7.17 16 8.83
C 10.67 5.33 5.33 16 16 0
D 3.5 0 10 3.5 13.5 10
E 7 16 16 23 23 0
F 9.5 3.5 13.5 13 23 10
G 12 3.5 18.5 15.5 30.5 15
H 6.33 23 23 29.33 29.33 0
I 5 23 25.5 28 30.5 2.5
J 8 28 30.5 36 38.5 2.5
K 9.17 29.33 29.33 31.5 38.5 0
The network is
Destinations B1 B2 B3 Supply
Sources units
A1 3 4 8 500
A2 2 6 10 600
A3 5 7 9 800
Demand 400 650 550
units
Required:
Max Z = 3x1-5x2
Subject to:
7x1- 2x2 ≥ -14
x1 ≤4
x2≤ 6
3x1+ 2x2 ≤ 18
and x1 ≥ 0, x2 ≥ 0
Required:
a) Put the LPM in its canonical form.
b) using the suitable simplex method, solve the
LPM.
c) Determine the dual problem and find the
optimum value for its objective function.
4) If you have the following LPM:
Max Z = 2x1 + 4x2
Subject to:
2x1+ x2 ≤ 400
2x1+ 2x2 ≤ 500
x1 ≤ 300
190
x2≤ 400
and x1 ≥ 0, x2 ≥ 0
Required:
a- By using the suitable simplex method, solve
the LPM.
b- Find the dual problem and determinates
optimum solution tableau .
c- By using two different methods, determine the
shadow prices.
5) If you have the following LPM:
Min Z = -x2
Subject to:
x 1 + x2 ≤ 2
| x1 - x2 | ≤ 1
and x1 ≥ 0, x2 ≥ 0
Required:
a) Put the LPM in its canonical form.
b) By using the suitable simplex method, solve
the LPM.
c) Determine the dual problem and find the
optimum value for its objective function.
191
References
Alrobh, M. Abd-elhamid.Operations Research: An
introduction. 2015, Faculty of commerce, South
Valley University
Hillier, Frederick S. and Lieberman, Gerald J.
Introduction to Operations Research.
SeventhEdition, 2000
Paul, A. Jensen and Jonathan, F. Bard.Operations
Research Models and Methods.John Wiley and
Sons 2003.
Reid, Dan R. and Sanders, Nada R. Operations
Management. 4th Edition, Wiley 2010
Taha,Hamdy A.Operations Research: An
Introduction. Eighth Edition, Pearson Prentice Hall
2007
192