Chapter 1
Introduction to
Quantitative Analysis
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Learning Objectives
After completing this chapter, students will be able to:
1. Describe the quantitative analysis approach
2. Understand the application of quantitative
analysis in a real situation
3. Describe the use of modeling in quantitative
analysis
4. Use computers and spreadsheet models to
perform quantitative analysis
5. Discuss possible problems in using
quantitative analysis
6. Perform a break-even analysis
© 2009 Prentice-Hall, Inc. 1–2
Chapter Outline
1.1 Introduction
1.2 What Is Quantitative Analysis?
1.3 The Quantitative Analysis Approach
1.4 How to Develop a Quantitative Analysis
Model
1.5 The Role of Computers and Spreadsheet
Models in the Quantitative Analysis
Approach
1.6 Possible Problems in the Quantitative
Analysis Approach
1.7 Implementation — Not Just the Final Step
© 2009 Prentice-Hall, Inc. 1–3
Introduction
Mathematical tools have been used for
thousands of years
Quantitative analysis can be applied to
a wide variety of problems
It’s not enough to just know the
mathematics of a technique
One must understand the specific
applicability of the technique, its
limitations, and its assumptions
© 2009 Prentice-Hall, Inc. 1–4
Examples of Quantitative Analyses
Taco Bell saved over $150 million using
forecasting and scheduling quantitative
analysis models
NBC television increased revenues by
over $200 million by using quantitative
analysis to develop better sales plans
Continental Airlines saved over $40
million using quantitative analysis
models to quickly recover from weather
delays and other disruptions
© 2009 Prentice-Hall, Inc. 1–5
What is Quantitative Analysis?
Quantitative analysis is a scientific approach
to managerial decision making whereby raw
data are processed and manipulated
resulting in meaningful information
Quantitative Meaningful
Raw Data Analysis Information
© 2009 Prentice-Hall, Inc. 1–6
What is Quantitative Analysis?
Quantitative factors might be different
investment alternatives, interest rates,
inventory levels, demand, or labor cost
Qualitative factors such as the weather,
state and federal legislation, and
technology breakthroughs should also be
considered
Information may be difficult to quantify but
can affect the decision-making process
© 2009 Prentice-Hall, Inc. 1–7
The Quantitative Analysis Approach
Defining the Problem
Developing a Model
Acquiring Input Data
Developing a Solution
Testing the Solution
Analyzing the Results
Figure 1.1 Implementing the Results
© 2009 Prentice-Hall, Inc. 1–8
Defining the Problem
Need to develop a clear and concise
statement that gives direction and meaning
to the following steps
This may be the most important and difficult
step
It is essential to go beyond symptoms and
identify true causes
May be necessary to concentrate on only a few
of the problems – selecting the right problems
is very important
Specific and measurable objectives may have to
be developed
© 2009 Prentice-Hall, Inc. 1–9
Developing a Model
Quantitative analysis models are realistic,
solvable, and understandable mathematical
representations of a situation
+ b 1X
$ Sales
b
Y= 0
$ Advertising
There are different types of models
Scale Schematic
models models
© 2009 Prentice-Hall, Inc. 1 – 10
Developing a Model
Models generally contain variables
(controllable and uncontrollable) and
parameters
Controllable variables are generally the
decision variables and are generally
unknown
Parameters are known quantities that
are a part of the problem
© 2009 Prentice-Hall, Inc. 1 – 11
Acquiring Input Data
Input data must be accurate – GIGO rule
Garbage
In
Process
Garbage
Out
Data may come from a variety of sources such as
company reports, company documents, interviews,
on-site direct measurement, or statistical sampling
© 2009 Prentice-Hall, Inc. 1 – 12
Developing a Solution
The best (optimal) solution to a problem
is found by manipulating the model
variables until a solution is found that is
practical and can be implemented
Common techniques are
Solving equations
Trial and error – trying various approaches
and picking the best result
Complete enumeration – trying all possible
values
Using an algorithm – a series of repeating
steps to reach a solution
© 2009 Prentice-Hall, Inc. 1 – 13
Testing the Solution
Both input data and the model should be
tested for accuracy before analysis and
implementation
New data can be collected to test the model
Results should be logical, consistent, and
represent the real situation
© 2009 Prentice-Hall, Inc. 1 – 14
Analyzing the Results
Determine the implications of the solution
Implementing results often requires change
in an organization
The impact of actions or changes needs to
be studied and understood before
implementation
Sensitivity analysis determines how much
the results of the analysis will change if
the model or input data changes
Sensitive models should be very thoroughly
tested
© 2009 Prentice-Hall, Inc. 1 – 15
Implementing the Results
Implementation incorporates the solution
into the company
Implementation can be very difficult
People can resist changes
Many quantitative analysis efforts have failed
because a good, workable solution was not
properly implemented
Changes occur over time, so even
successful implementations must be
monitored to determine if modifications are
necessary
© 2009 Prentice-Hall, Inc. 1 – 16
Modeling in the Real World
Quantitative analysis models are used
extensively by real organizations to solve
real problems
In the real world, quantitative analysis
models can be complex, expensive, and
difficult to sell
Following the steps in the process is an
important component of success
© 2009 Prentice-Hall, Inc. 1 – 17
How To Develop a Quantitative
Analysis Model
An important part of the quantitative
analysis approach
Let’s look at a simple mathematical
model of profit
Profit = Revenue – Expenses
© 2009 Prentice-Hall, Inc. 1 – 18
How To Develop a Quantitative
Analysis Model
Expenses can be represented as the sum of fixed and
variable costs and variable costs are the product of
unit costs times the number of units
Profit = Revenue – (Fixed cost + Variable cost)
Profit = (Selling price per unit)(number of units
sold) – [Fixed cost + (Variable costs per
unit)(Number of units sold)]
Profit = sX – [f + vX]
Profit = sX – f – vX
where
s = selling price per unit v = variable cost per unit
f = fixed cost X = number of units sold
© 2009 Prentice-Hall, Inc. 1 – 19
How To Develop a Quantitative
Analysis Model
Expenses can be represented as the sum of fixed and
variable costs and variable costs are the
The parameters product
of this modelof
unit costs times the number
are f, v,ofand
units
s as these are the
inputscost
Profit = Revenue – (Fixed inherent in the cost)
+ Variable model
Profit = (Selling priceThe
perdecision variable
unit)(number of of
units
sold) – [Fixedinterest
cost +is(Variable
X costs per
unit)(Number of units sold)]
Profit = sX – [f + vX]
Profit = sX – f – vX
where
s = selling price per unit v = variable cost per unit
f = fixed cost X = number of units sold
© 2009 Prentice-Hall, Inc. 1 – 20
Pritchett’s Precious Time Pieces
The company buys, sells, and repairs old clocks.
Rebuilt springs sell for $10 per unit. Fixed cost of
equipment to build springs is $1,000. Variable cost
for spring material is $5 per unit.
s = 10 f = 1,000 v=5
Number of spring sets sold = X
Profits = sX – f – vX
If sales = 0, profits = –$1,000
If sales = 1,000, profits = [(10)(1,000) – 1,000 – (5)(1,000)]
= $4,000
© 2009 Prentice-Hall, Inc. 1 – 21
Pritchett’s Precious Time Pieces
Companies are often interested in their break-even
point (BEP). The BEP is the number of units sold
that will result in $0 profit.
0 = sX – f – vX, or 0 = (s – v)X – f
Solving for X, we have
f = (s – v)X
f
X=
s–v
Fixed cost
BEP =
(Selling price per unit) – (Variable cost per unit)
© 2009 Prentice-Hall, Inc. 1 – 22
Pritchett’s Precious Time Pieces
Companies are often interested in their break-even
point (BEP). The BEP is the number of units sold
BEP for Pritchett’s Precious Time Pieces
that will result in $0 profit.
BEP
0= sX –= f$1,000/($10
– vX, or – 0$5) = –200
= (s v)Xunits
–f
Salesfor
Solving of less
X, wethan
have 200 units of rebuilt springs
will result in a loss
f = (s – v)X
Sales of over 200 units of rebuilt springs will
result in a profit X = f
s–v
Fixed cost
BEP =
(Selling price per unit) – (Variable cost per unit)
© 2009 Prentice-Hall, Inc. 1 – 23
Advantages of Mathematical Modeling
1. Models can accurately represent reality
2. Models can help a decision maker
formulate problems
3. Models can give us insight and information
4. Models can save time and money in
decision making and problem solving
5. A model may be the only way to solve large
or complex problems in a timely fashion
6. A model can be used to communicate
problems and solutions to others
© 2009 Prentice-Hall, Inc. 1 – 24
Models Categorized by Risk
Mathematical models that do not involve
risk are called deterministic models
We know all the values used in the model
with complete certainty
Mathematical models that involve risk,
chance, or uncertainty are called
probabilistic models
Values used in the model are estimates
based on probabilities
© 2009 Prentice-Hall, Inc. 1 – 25
Computers and Spreadsheet Models
QM for Windows
An easy to use
decision support
system for use in
POM and QM
courses
This is the main
menu of
quantitative
models
Program 1.1
© 2009 Prentice-Hall, Inc. 1 – 26
Computers and Spreadsheet Models
Excel QM’s Main Menu (2003)
Works automatically within Excel spreadsheets
Program 1.2A
© 2009 Prentice-Hall, Inc. 1 – 27
Computers and Spreadsheet Models
Excel QM’s
Main Menu
(2007)
Program 1.2B
© 2009 Prentice-Hall, Inc. 1 – 28
Computers and Spreadsheet Models
Excel QM
for the
Break-
Even
Problem
Program 1.3A
© 2009 Prentice-Hall, Inc. 1 – 29
Computers and Spreadsheet Models
Excel QM
Solution
to the
Break-
Even
Problem
Program 1.3B
© 2009 Prentice-Hall, Inc. 1 – 30
Computers and Spreadsheet Models
Using
Goal Seek
in the
Break-
Even
Problem
Program 1.4
© 2009 Prentice-Hall, Inc. 1 – 31
Possible Problems in the
Quantitative Analysis Approach
Defining the problem
Problems are not easily identified
Conflicting viewpoints
Impact on other departments
Beginning assumptions
Solution outdated
Developing a model
Fitting the textbook models
Understanding the model
© 2009 Prentice-Hall, Inc. 1 – 32
Possible Problems in the
Quantitative Analysis Approach
Acquiring input data
Using accounting data
Validity of data
Developing a solution
Hard-to-understand mathematics
Only one answer is limiting
Testing the solution
Analyzing the results
© 2009 Prentice-Hall, Inc. 1 – 33
Implementation –
Not Just the Final Step
Lack of commitment and resistance
to change
Management may fear the use of
formal analysis processes will
reduce their decision-making power
Action-oriented managers may want
“quick and dirty” techniques
Management support and user
involvement are important
© 2009 Prentice-Hall, Inc. 1 – 34
Implementation –
Not Just the Final Step
Lack of commitment by quantitative
analysts
An analysts should be involved with
the problem and care about the
solution
Analysts should work with users and
take their feelings into account
© 2009 Prentice-Hall, Inc. 1 – 35
Summary
Quantitative analysis is a scientific
approach to decision making
The approach includes
Defining the problem
Acquiring input data
Developing a solution
Testing the solution
Analyzing the results
Implementing the results
© 2009 Prentice-Hall, Inc. 1 – 36
Summary
Potential problems include
Conflicting viewpoints
The impact on other departments
Beginning assumptions
Outdated solutions
Fitting textbook models
Understanding the model
Acquiring good input data
Hard-to-understand mathematics
Obtaining only one answer
Testing the solution
Analyzing the results
© 2009 Prentice-Hall, Inc. 1 – 37
Summary
Implementation is not the final step
Problems can occur because of
Lack of commitment to the approach
Resistance to change
© 2009 Prentice-Hall, Inc. 1 – 38
Chapter 7
Linear Programming Models:
Graphical and Computer
Methods
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Introduction
◼ Many management decisions involve trying to
make the most effective use of limited resources
◼ Machinery, labor, money, time, warehouse space, raw
materials
◼ Linear programming (LP) is a widely used
mathematical modeling technique designed to
help managers in planning and decision making
relative to resource allocation
◼ Belongs to the broader field of mathematical
programming
◼ In this sense, programming refers to modeling and
solving a problem mathematically
© 2009 Prentice-Hall, Inc. 7–2
Requirements of a Linear
Programming Problem
◼ LP has been applied in many areas over
the past 50 years
◼ All LP problems have 4 properties in
common
1. All problems seek to maximize or minimize
some quantity (the objective function)
2. The presence of restrictions or constraints that
limit the degree to which we can pursue our
objective
3. There must be alternative courses of action to
choose from
4. The objective and constraints in problems
must be expressed in terms of linear equations
or inequalities
© 2009 Prentice-Hall, Inc. 7–3
Examples of Successful
LP Applications
◼ Development of a production schedule that will
◼ satisfy future demands for a firm’s production
◼ while minimizing total production and inventory costs
◼ Determination of grades of petroleum products to yield
the maximum profit
◼ Selection of different blends of raw materials to feed
mills to produce finished feed combinations at
minimum cost
◼ Determination of a distribution system that will
minimize total shipping cost from several warehouses
to various market locations
© 2009 Prentice-Hall, Inc. 7–4
LP Properties and Assumptions
PROPERTIES OF LINEAR PROGRAMS
1. One objective function
2. One or more constraints
3. Alternative courses of action
4. Objective function and constraints are linear
ASSUMPTIONS OF LP
1. Certainty
2. Proportionality
3. Additivity
4. Divisibility
5. Nonnegative variables
Table 7.1
© 2009 Prentice-Hall, Inc. 7–5
Basic Assumptions of LP
◼ We assume conditions of certainty exist and numbers
in the objective and constraints are known with
certainty and do not change during the period being
studied
◼ We assume proportionality exists in the objective and
constraints
▪ constancy between production increases and resource
utilization – if 1 unit needs 3 hours then 10 require 30
hours
◼ We assume additivity in that the total of all activities
equals the sum of the individual activities
◼ We assume divisibility in that solutions need not be
whole numbers
◼ All answers or variables are nonnegative as we are
dealing with real physical quantities
© 2009 Prentice-Hall, Inc. 7–6
Formulating LP Problems
◼ Formulating a linear program involves
developing a mathematical model to represent
the managerial problem
◼ The steps in formulating a linear program are
1. Completely understand the managerial
problem being faced
2. Identify the objective and constraints
3. Define the decision variables
4. Use the decision variables to write
mathematical expressions for the objective
function and the constraints
© 2009 Prentice-Hall, Inc. 7–7
Formulating LP Problems
◼ One of the most common LP applications is the
product mix problem
◼ Two or more products are produced using limited
resources such as personnel, machines, and raw
materials
◼ The profit that the firm seeks to maximize is
based on the profit contribution per unit of each
product
◼ The company would like to determine how many
units of each product it should produce so as to
maximize overall profit given its limited
resources
© 2009 Prentice-Hall, Inc. 7–8
Flair Furniture Company
◼ The Flair Furniture Company produces
inexpensive tables and chairs
◼ Processes are similar in that both require a certain
amount of hours of carpentry work and in the
painting and varnishing department
◼ Each table takes 4 hours of carpentry and 2 hours
of painting and varnishing
◼ Each chair requires 3 of carpentry and 1 hour of
painting and varnishing
◼ There are 240 hours of carpentry time available
and 100 hours of painting and varnishing
◼ Each table yields a profit of $70 and each chair a
profit of $50
© 2009 Prentice-Hall, Inc. 7–9
Flair Furniture Company
◼ The company wants to determine the best
combination of tables and chairs to produce to
reach the maximum profit
HOURS REQUIRED TO
PRODUCE 1 UNIT
(T) (C) AVAILABLE HOURS
DEPARTMENT TABLES CHAIRS THIS WEEK
Carpentry 4 3 240
Painting and varnishing 2 1 100
Profit per unit $70 $50
Table 7.2
© 2009 Prentice-Hall, Inc. 7 – 10
Flair Furniture Company
◼ The objective is to
Maximize profit
◼ The constraints are
1. The hours of carpentry time used cannot
exceed 240 hours per week
2. The hours of painting and varnishing time
used cannot exceed 100 hours per week
◼ The decision variables representing the actual
decisions we will make are
T = number of tables to be produced per week
C = number of chairs to be produced per week
© 2009 Prentice-Hall, Inc. 7 – 11
Flair Furniture Company
◼ We create the LP objective function in terms of T
and C
Maximize profit = $70T + $50C
◼ Develop mathematical relationships for the two
constraints
◼ For carpentry, total time used is
(4 hours per table)(Number of tables produced)
+ (3 hours per chair)(Number of chairs produced)
◼ We know that
Carpentry time used ≤ Carpentry time available
4T + 3C ≤ 240 (hours of carpentry time)
© 2009 Prentice-Hall, Inc. 7 – 12
Flair Furniture Company
◼ Similarly
Painting and varnishing time used
≤ Painting and varnishing time available
2 T + 1C ≤ 100 (hours of painting and varnishing time)
This means that each table produced
requires two hours of painting and
varnishing time
◼ Both of these constraints restrict production
capacity and affect total profit
© 2009 Prentice-Hall, Inc. 7 – 13
Flair Furniture Company
◼ The values for T and C must be nonnegative
T ≥ 0 (number of tables produced is greater
than or equal to 0)
C ≥ 0 (number of chairs produced is greater
than or equal to 0)
◼ The complete problem stated mathematically
Maximize profit = $70T + $50C
subject to
4T + 3C ≤ 240 (carpentry constraint)
2T + 1C ≤ 100 (painting and varnishing constraint)
T, C ≥ 0 (nonnegativity constraint)
© 2009 Prentice-Hall, Inc. 7 – 14
Graphical Solution to an LP Problem
◼ The easiest way to solve a small LP
problems is with the graphical solution
approach
◼ The graphical method only works when
there are just two decision variables
◼ When there are more than two variables, a
more complex approach is needed as it is
not possible to plot the solution on a two-
dimensional graph
◼ The graphical method provides valuable
insight into how other approaches work
© 2009 Prentice-Hall, Inc. 7 – 15
Graphical Representation of a
Constraint
C
100 –
– This Axis Represents the Constraint T ≥ 0
80 –
Number of Chairs
–
60 –
–
40 – This Axis Represents the
– Constraint C ≥ 0
20 –
–
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.1 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 16
Graphical Representation of a
Constraint
◼ The first step in solving the problem is to
identify a set or region of feasible
solutions
◼ To do this we plot each constraint
equation on a graph
◼ We start by graphing the equality portion
of the constraint equations
4T + 3C = 240
◼ We solve for the axis intercepts and draw
the line
© 2009 Prentice-Hall, Inc. 7 – 17
Graphical Representation of a
Constraint
◼ When Flair produces no tables, the
carpentry constraint is
4(0) + 3C = 240
3C = 240
C = 80
◼ Similarly for no chairs
4T + 3(0) = 240
4T = 240
T = 60
◼ This line is shown on the following graph
© 2009 Prentice-Hall, Inc. 7 – 18
Graphical Representation of a
Constraint
C
Graph of carpentry constraint equation
100 –
–
80 –
Number of Chairs
(T = 0, C = 80)
–
60 –
–
40 –
–
(T = 60, C = 0)
20 –
–
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.2 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 19
Graphical Representation of a
Constraint
C ◼ Any point on or below the constraint
plot will not violate the restriction
100 –
◼ Any point above the plot will violate
– the restriction
80 –
Number of Chairs
–
60 –
–
(30, 40) (70, 40)
40 –
–
20 –
– (30, 20)
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.3 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 20
Graphical Representation of a
Constraint
◼ The point (30, 40) lies on the plot and
exactly satisfies the constraint
4(30) + 3(40) = 240
◼ The point (30, 20) lies below the plot and
satisfies the constraint
4(30) + 3(20) = 180
◼ The point (30, 40) lies above the plot and
does not satisfy the constraint
4(70) + 3(40) = 400
© 2009 Prentice-Hall, Inc. 7 – 21
Graphical Representation of a
Constraint
C
100 – (T = 0, C = 100)
–
80 –
Number of Chairs
Graph of painting and varnishing
– constraint equation
60 –
–
40 –
–
(T = 50, C = 0)
20 –
–
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.4 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 22
Graphical Representation of a
Constraint
◼ To produce tables and chairs, both
departments must be used
◼ We need to find a solution that satisfies both
constraints simultaneously
◼ A new graph shows both constraint plots
◼ The feasible region (or area of feasible
solutions) is where all constraints are satisfied
◼ Any point inside this region is a feasible
solution
◼ Any point outside the region is an infeasible
solution
© 2009 Prentice-Hall, Inc. 7 – 23
Graphical Representation of a
Constraint
C
◼ Feasible solution region for Flair Furniture
100 –
–
80 –
Number of Chairs
Painting/Varnishing Constraint
–
60 –
–
40 –
–
Carpentry Constraint
20 – Feasible
– Region
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.5 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 24
Graphical Representation of a
Constraint
◼ For the point (30, 20)
Carpentry 4T + 3C ≤ 240 hours available
constraint (4)(30) + (3)(20) = 180 hours used ✓
Painting 2T + 1C ≤ 100 hours available
constraint (2)(30) + (1)(20) = 80 hours used
✓
◼ For the point (70, 40)
Carpentry 4T + 3C ≤ 240 hours available
constraint (4)(70) + (3)(40) = 400 hours used
Painting 2T + 1C ≤ 100 hours available
constraint (2)(70) + (1)(40) = 180 hours used
© 2009 Prentice-Hall, Inc. 7 – 25
Graphical Representation of a
Constraint
◼ For the point (50, 5)
Carpentry 4T + 3C ≤ 240 hours available
constraint (4)(50) + (3)(5) = 215 hours used ✓
Painting 2T + 1C ≤ 100 hours available
constraint (2)(50) + (1)(5) = 105 hours used
© 2009 Prentice-Hall, Inc. 7 – 26
Isoprofit Line Solution Method
◼ Once the feasible region has been graphed, we need
to find the optimal solution from the many possible
solutions
◼ The speediest way to do this is to use the isoprofit
line method
◼ Starting with a small but possible profit value, we
graph the objective function
◼ We move the objective function line in the direction of
increasing profit while maintaining the slope
◼ The last point it touches in the feasible region is the
optimal solution
© 2009 Prentice-Hall, Inc. 7 – 27
Isoprofit Line Solution Method
◼ For Flair Furniture, choose a profit of $2,100
◼ The objective function is then
$2,100 = 70T + 50C
◼ Solving for the axis intercepts, we can draw the
graph
◼ This is obviously not the best possible solution
◼ Further graphs can be created using larger profits
◼ The further we move from the origin while
maintaining the slope and staying within the
boundaries of the feasible region, the larger the
profit will be
◼ The highest profit ($4,100) will be generated when
the isoprofit line passes through the point (30, 40)
© 2009 Prentice-Hall, Inc. 7 – 28
Isoprofit Line Solution Method
C
◼ Isoprofit line at $2,100
100 –
–
80 –
Number of Chairs
–
60 –
– $2,100 = $70T + $50C
(0, 42)
40 –
–
(30, 0)
20 –
–
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.6 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 29
Isoprofit Line Solution Method
C
◼ Four isoprofit lines
100 –
–
$3,500 = $70T + $50C
80 –
Number of Chairs
– $2,800 = $70T + $50C
60 –
– $2,100 = $70T + $50C
40 –
– $4,100 = $70T + $50C
20 –
–
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.7 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 30
Isoprofit Line Solution Method
C
◼ Optimal solution to the
100 – Flair Furniture problem
–
80 –
Number of Chairs
Maximum Profit Line
–
60 – Optimal Solution Point
– (T = 30, C = 40)
40 –
– $4,100 = $70T + $50C
20 –
–
|– | | | | | | | | | | |
0 20 40 60 80 100 T
Figure 7.8 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 31
Corner Point Solution Method
◼ A second approach to solving LP problems
employs the corner point method
◼ It involves looking at the profit at every
corner point of the feasible region
◼ The mathematical theory behind LP is that
the optimal solution must lie at one of the
corner points, or extreme point, in the
feasible region
◼ For Flair Furniture, the feasible region is a
four-sided polygon with four corner points
labeled 1, 2, 3, and 4 on the graph
© 2009 Prentice-Hall, Inc. 7 – 32
Corner Point Solution Method
C
◼ Four corner points of
100 – the feasible region
2 –
80 –
Number of Chairs
–
60 –
–
3
40 –
–
20 –
–
1 |– | | | | | | | | | | |
0 20 40 60 80 100 T
4
Figure 7.9 Number of Tables
© 2009 Prentice-Hall, Inc. 7 – 33
Corner Point Solution Method
Point 1 : (T = 0, C = 0) Profit = $70(0) + $50(0) = $0
Point 2 : (T = 0, C = 80) Profit = $70(0) + $50(80) = $4,000
Point 4 : (T = 50, C = 0) Profit = $70(50) + $50(0) = $3,500
Point 3 : (T = 30, C = 40) Profit = $70(30) + $50(40) = $4,100
◼ Because Point 3 returns the highest profit, this
is the optimal solution
◼ To find the coordinates for Point 3 accurately we
have to solve for the intersection of the two
constraint lines
◼ The details of this are on the following slide
© 2009 Prentice-Hall, Inc. 7 – 34
Corner Point Solution Method
◼ Using the simultaneous equations method, we
multiply the painting equation by –2 and add it to
the carpentry equation
4T + 3C = 240 (carpentry line)
– 4T – 2C = –200 (painting line)
C = 40
◼ Substituting 40 for C in either of the original
equations allows us to determine the value of T
4T + (3)(40) = 240 (carpentry line)
4T + 120 = 240
T = 30
© 2009 Prentice-Hall, Inc. 7 – 35
Summary of Graphical Solution
Methods
ISOPROFIT METHOD
1. Graph all constraints and find the feasible region.
2. Select a specific profit (or cost) line and graph it to find the slope.
3. Move the objective function line in the direction of increasing profit (or
decreasing cost) while maintaining the slope. The last point it touches in the
feasible region is the optimal solution.
4. Find the values of the decision variables at this last point and compute the
profit (or cost).
CORNER POINT METHOD
1. Graph all constraints and find the feasible region.
2. Find the corner points of the feasible reason.
3. Compute the profit (or cost) at each of the feasible corner points.
4. select the corner point with the best value of the objective function found in
Step 3. This is the optimal solution.
Table 7.3
© 2009 Prentice-Hall, Inc. 7 – 36
Solving Flair Furniture’s LP Problem
Using QM for Windows and Excel
◼ Most organizations have access to
software to solve big LP problems
◼ While there are differences between
software implementations, the approach
each takes towards handling LP is
basically the same
◼ Once you are experienced in dealing with
computerized LP algorithms, you can
easily adjust to minor changes
© 2009 Prentice-Hall, Inc. 7 – 37
Using QM for Windows
◼ First select the Linear Programming module
◼ Specify the number of constraints (non-negativity
is assumed)
◼ Specify the number of decision variables
◼ Specify whether the objective is to be maximized
or minimized
◼ For the Flair Furniture problem there are two
constraints, two decision variables, and the
objective is to maximize profit
© 2009 Prentice-Hall, Inc. 7 – 38
Using QM for Windows
◼ Computer screen for input of data
Program 7.1A
© 2009 Prentice-Hall, Inc. 7 – 39
Using QM for Windows
◼ Computer screen for input of data
Program 7.1B
© 2009 Prentice-Hall, Inc. 7 – 40
Using QM for Windows
◼ Computer screen for output of solution
Program 7.1C
© 2009 Prentice-Hall, Inc. 7 – 41
Using QM for Windows
◼ Graphical output of solution
Program 7.1D
© 2009 Prentice-Hall, Inc. 7 – 42
Solving Minimization Problems
◼ Many LP problems involve minimizing an
objective such as cost instead of
maximizing a profit function
◼ Minimization problems can be solved
graphically by first setting up the feasible
solution region and then using either the
corner point method or an isocost line
approach (which is analogous to the
isoprofit approach in maximization
problems) to find the values of the decision
variables (e.g., X1 and X2) that yield the
minimum cost
© 2009 Prentice-Hall, Inc. 7 – 43
Holiday Meal Turkey Ranch
◼ The Holiday Meal Turkey Ranch is considering
buying two different brands of turkey feed and
blending them to provide a good, low-cost diet for
its turkeys
Let
X1 = number of pounds of brand 1 feed purchased
X2 = number of pounds of brand 2 feed purchased
Minimize cost (in cents) = 2X1 + 3X2
subject to:
5X1 + 10X2 ≥ 90 ounces (ingredient constraint A)
4X1 + 3X2 ≥ 48 ounces (ingredient constraint B)
0.5X1 ≥ 1.5 ounces (ingredient constraint C)
X1 ≥ 0 (nonnegativity constraint)
X2 ≥ 0 (nonnegativity constraint)
© 2009 Prentice-Hall, Inc. 7 – 44
Holiday Meal Turkey Ranch
◼ Holiday Meal Turkey Ranch data
COMPOSITION OF EACH POUND
OF FEED (OZ.) MINIMUM MONTHLY
REQUIREMENT PER
INGREDIENT BRAND 1 FEED BRAND 2 FEED TURKEY (OZ.)
A 5 10 90
B 4 3 48
C 0.5 0 1.5
Cost per pound 2 cents 3 cents
Table 7.4
© 2009 Prentice-Hall, Inc. 7 – 45
Holiday Meal Turkey Ranch
X2
◼ Using the corner –
point method
◼ First we construct 20 –
the feasible Ingredient C Constraint
Pounds of Brand 2
solution region
15 –
◼ The optimal Feasible Region
solution will lie at a
10 –
on of the corners
as it would in a Ingredient B Constraint
maximization 5–
Ingredient A Constraint
b
problem
0 |– | | | c | |
5 10 15 20 25 X1
Figure 7.10 Pounds of Brand 1
© 2009 Prentice-Hall, Inc. 7 – 46
Holiday Meal Turkey Ranch
◼ We solve for the values of the three corner points
◼ Point a is the intersection of ingredient constraints
C and B
4X1 + 3X2 = 48
X1 = 3
◼ Substituting 3 in the first equation, we find X2 = 12
◼ Solving for point b with basic algebra we find X1 =
8.4 and X2 = 4.8
◼ Solving for point c we find X1 = 18 and X2 = 0
© 2009 Prentice-Hall, Inc. 7 – 47
Holiday Meal Turkey Ranch
◼ Substituting these value back into the objective
function we find
Cost = 2X1 + 3X2
Cost at point a = 2(3) + 3(12) = 42
Cost at point b = 2(8.4) + 3(4.8) = 31.2
Cost at point c = 2(18) + 3(0) = 36
◼ The lowest cost solution is to purchase 8.4
pounds of brand 1 feed and 4.8 pounds of brand 2
feed for a total cost of 31.2 cents per turkey
© 2009 Prentice-Hall, Inc. 7 – 48
Holiday Meal Turkey Ranch
X2
◼ Using the isocost –
approach Feasible Region
◼ Choosing an 20 –
initial cost of 54
Pounds of Brand 2
cents, it is clear
15 –
improvement is
possible
10 –
5–
(X1 = 8.4, X2 = 4.8)
0 |– | | | | |
5 10 15 20 25 X1
Figure 7.11 Pounds of Brand 1
© 2009 Prentice-Hall, Inc. 7 – 49
Holiday Meal Turkey Ranch
◼ QM for Windows can also be used to solve the
Holiday Meal Turkey Ranch problem
Program 7.3
© 2009 Prentice-Hall, Inc. 7 – 50
Employee Scheduling Applications
◼ Assignment Problems
◼ Involve determining the most efficient way to
assign resources to tasks
◼ Objective may be to minimize travel times or
maximize assignment effectiveness
◼ Assignment problems are unique because they
have a coefficient of 0 or 1 associated with
each variable in the LP constraints and the
right-hand side of each constraint is always
equal to 1
© 2009 Prentice-Hall, Inc. 7 – 51
Employee Scheduling Applications
◼ Ivan and Ivan law firm maintains a large staff of
young attorneys
◼ Ivan wants to make lawyer-to-client assignments
in the most effective manner
◼ He identifies four lawyers who could possibly be
assigned new cases
◼ Each lawyer can handle one new client
◼ The lawyers have different skills and special
interests
◼ The following table summarizes the lawyers
estimated effectiveness on new cases
© 2009 Prentice-Hall, Inc. 7 – 52
Employee Scheduling Applications
◼ Effectiveness ratings
CLIENT’S CASE
CORPORATE
LAWYER DIVORCE MERGER EMBEZZLEMENT EXHIBITIONISM
Adams 6 2 8 5
Brooks 9 3 5 8
Carter 4 8 3 4
Darwin 6 7 6 4
1 if attorney i is assigned to case j
Let Xij =
0 otherwise
where i = 1, 2, 3, 4 stands for Adams, Brooks,
Carter, and Darwin respectively
j = 1, 2, 3, 4 stands for divorce, merger,
embezzlement, and exhibitionism © 2009 Prentice-Hall, Inc. 7 – 53
Employee Scheduling Applications
◼ The LP formulation is
6X=
Maximize effectiveness 11 + 2X12 + 8X13 + 5X14 + 9X21 + 3X22
+ 5X23 + 8X24 + 4X31 + 8X32 + 3X33 + 4X
+ 6X41 + 7X42 + 6X43 + 4X44
subject to X11 + X21 + X31 + X41 = 1 (divorce case)
X12 + X22 + X32 + X42 = 1 (merger)
X13 + X23 + X33 + X43 = 1 (embezzlement)
X14 + X24 + X34 + X44 = 1 (exhibitionism)
X11 + X12 + X13 + X14 = 1 (Adams)
X21 + X22 + X23 + X24 = 1 (Brook)
X31 + X32 + X33 + X34 = 1 (Carter)
X41 + X42 + X43 + X44 = 1 (Darwin)
© 2009 Prentice-Hall, Inc. 7 – 54
Employee Scheduling Applications
◼ Solving Ivan and Ivan’s assignment scheduling
LP problem using QM for Windows
Program 8.4
© 2009 Prentice-Hall, Inc. 7 – 55
Four Special Cases in LP
◼ Four special cases and difficulties arise at
times when using the graphical approach
to solving LP problems
◼ Infeasibility
◼ Unboundedness
◼ Redundancy
◼ Alternate Optimal Solutions
© 2009 Prentice-Hall, Inc. 7 – 56
Four Special Cases in LP
◼ No feasible solution
◼ Exists when there is no solution to the
problem that satisfies all the constraint
equations
◼ No feasible solution region exists
◼ This is a common occurrence in the real world
◼ Generally one or more constraints are relaxed
until a solution is found
© 2009 Prentice-Hall, Inc. 7 – 57
Four Special Cases in LP
◼ A problem with no feasible
X1+2X2 <=6
solution 2X1+X2 <=8
X2
X1 >=7
8–
–
6–
–
Region Satisfying
4– Third Constraint
–
2–
–
0– | | | | | | | | | |
2 4 6 8 X1
Figure 7.12 Region Satisfying First Two Constraints
© 2009 Prentice-Hall, Inc. 7 – 58
Four Special Cases in LP
◼ Unboundedness
◼ Sometimes a linear program will not have a
finite solution
◼ In a maximization problem, one or more
solution variables, and the profit, can be made
infinitely large without violating any
constraints
◼ In a graphical solution, the feasible region will
be open ended
◼ This usually means the problem has been
formulated improperly
© 2009 Prentice-Hall, Inc. 7 – 59
Four Special Cases in LP
◼ A solution region unbounded to the right
X2
X1 ≥ 5
15 –
X2 ≤ 10
10 –
Feasible Region
5–
X1 + 2X2 ≥ 15
0 |– | | | |
5 10 15 X1
Figure 7.13
© 2009 Prentice-Hall, Inc. 7 – 60
Four Special Cases in LP
◼ Redundancy
◼ A redundant constraint is one that does not
affect the feasible solution region
◼ One or more constraints may be more binding
◼ This is a very common occurrence in the real
world
◼ It causes no particular problems, but
eliminating redundant constraints simplifies
the model
© 2009 Prentice-Hall, Inc. 7 – 61
Four Special Cases in LP
◼ A problem with X2
a redundant 30 –
constraint
25 –
2X1 + X2 ≤ 30
20 –
Redundant
Constraint
15 –
X1 ≤ 25
10 –
X1 + X2 ≤ 20
Feasible
5–
Region
0– | | | | | |
Figure 7.14 5 10 15 20 25 30 X1
© 2009 Prentice-Hall, Inc. 7 – 62
Four Special Cases in LP
◼ Alternate Optimal Solutions
◼ Occasionally two or more optimal solutions
may exist
◼ Graphically this occurs when the objective
function’s isoprofit or isocost line runs
perfectly parallel to one of the constraints
◼ This actually allows management great
flexibility in deciding which combination to
select as the profit is the same at each
alternate solution
© 2009 Prentice-Hall, Inc. 7 – 63
Four Special Cases in LP
◼ Example of X2
alternate 8–
optimal
solutions 7–
6 –A
Optimal Solution Consists of All
5– Combinations of X1 and X2 Along
Maximize 3X1 + 2X2 the AB Segment
Subj. To: 6X1 + 4X2 < 24 4–
X1 <3 3– Isoprofit Line for $8
X1, X 2 > 0
2–
B Isoprofit Line for $12
1 – Feasible Overlays Line Segment AB
Region
0– | | | | | | | |
Figure 7.15 1 2 3 4 5 6 7 8 X1
© 2009 Prentice-Hall, Inc. 7 – 64
Sensitivity Analysis
◼ Optimal solutions to LP problems thus far have
been found under what are called deterministic
assumptions
◼ This means that we assume complete certainty in
the data and relationships of a problem
◼ But in the real world, conditions are dynamic and
changing
◼ We can analyze how sensitive a deterministic
solution is to changes in the assumptions of the
model
◼ This is called sensitivity analysis, postoptimality
analysis, parametric programming, or optimality
analysis
© 2009 Prentice-Hall, Inc. 7 – 65
Sensitivity Analysis
◼ Sensitivity analysis often involves a series of
what-if? questions concerning constraints,
variable coefficients, and the objective function
◼ What if the profit for product 1 increases by 10%?
◼ What if less advertising money is available?
◼ One way to do this is the trial-and-error method
where values are changed and the entire model is
resolved
◼ The preferred way is to use an analytic
postoptimality analysis
◼ After a problem has been solved, we determine a
range of changes in problem parameters that will not
affect the optimal solution or change the variables in
the solution without re-solving the entire problem
© 2009 Prentice-Hall, Inc. 7 – 66
Sensitivity Analysis
◼ Sensitivity analysis can be used to deal not
only with errors in estimating input
parameters to the LP model but also with
management’s experiments with possible
future changes in the firm that may affect
profits.
© 2009 Prentice-Hall, Inc. 7 – 67
High Note Sound Company
◼ The High Note Sound Company manufactures
quality CD players and stereo receivers
◼ Products require a certain amount of skilled
artisanship which is in limited supply
◼ The firm has formulated the following product mix
LP model
Maximize profit = $50X1 + $120X2
Subject to 2X1 + 4X2 ≤ 80 (hours of electrician’s
time available)
3X1 + 1X2 ≤ 60 (hours of audio
technician’s time
available)
X1, X2 ≥ 0
© 2009 Prentice-Hall, Inc. 7 – 68
High Note Sound Company
◼ The High Note Sound Company graphical solution
X2
(receivers)
60 –
– Optimal Solution at Point a
X1 = 0 CD Players
40 – X2 = 20 Receivers
Profits = $2,400
a = (0, 20) –
b = (16, 12)
20 –
Isoprofit Line: $2,400 = 50X1 + 120X2
10 –
0– | | | | | |
10 20 30 40 50 60 X1
Figure 7.16 c = (20, 0) (CD players)
© 2009 Prentice-Hall, Inc. 7 – 69
Changes in the
Objective Function Coefficient
◼ In real-life problems, contribution rates in the
objective functions fluctuate periodically
◼ Graphically, this means that although the feasible
solution region remains exactly the same, the
slope of the isoprofit or isocost line will change
◼ We can often make modest increases or
decreases in the objective function coefficient of
any variable without changing the current optimal
corner point
◼ We need to know how much an objective function
coefficient can change before the optimal solution
would be at a different corner point
© 2009 Prentice-Hall, Inc. 7 – 70
Changes in the
Objective Function Coefficient
◼ Changes in the receiver contribution coefficients
X2
40 –
Profit Line for 50X1 + 80X2
(Passes through Point b)
30 –
Profit Line for 50X1 + 120X2
(Passes through Point a)
20 –
b
a Profit Line for 50X1 + 150X2
10 – (Passes through Point a)
0– | | c | | | |
10 20 30 40 50 60 X1
Figure 7.17
© 2009 Prentice-Hall, Inc. 7 – 71
QM for Windows and Changes in
Objective Function Coefficients
◼ Input and sensitivity analysis for High Note Sound
data
Program 7.5A
Program 7.5B
© 2009 Prentice-Hall, Inc. 7 – 72
Changes in the
Technological Coefficients
◼ Changes in the technological coefficients often
reflect changes in the state of technology
◼ If the amount of resources needed to produce a
product changes, coefficients in the constraint
equations will change
◼ This does not change the objective function, but
it can produce a significant change in the shape
of the feasible region
◼ This may cause a change in the optimal solution
© 2009 Prentice-Hall, Inc. 7 – 73
Changes in the
Technological Coefficients
◼ Change in the technological coefficients for the
High Note Sound Company
(a) Original Problem (b) Change in Circled (c) Change in Circled
Coefficient Coefficient
X2 X2 X2
60 – 60 – 60 –
Stereo Receivers
3X1 + 1X2 ≤ 60 2 X1 + 1X2 ≤ 60 3X1 + 1X2 ≤ 60
40 – 40 – 40 –
Optimal Still Optimal
Solution Optimal Solution
20 –a 20 –a 20 –
d
b 2X1 + 4X2 ≤ 80 2X1 + 4X2 ≤ 80 16 f g 2X1 + 5 X2 ≤ 80
– |c | | – | e | | | – |c | |
0 20 40 X1 0 20 30 40 X1 0 20 40 X1
CD Players
Figure 7.18
© 2009 Prentice-Hall, Inc. 7 – 74
Changes in Resources or
Right-Hand-Side Values
◼ The right-hand-side values of the constraints
often represent resources available to the firm
◼ If additional resources were available, a higher
total profit could be realized
◼ Sensitivity analysis about resources will help
answer questions such as:
◼ How much should the company be willing to pay
for additional hours?
◼ Is it profitable to have some electricians work
overtime?
◼ Should we be willing to pay for more audio
technician time?
© 2009 Prentice-Hall, Inc. 7 – 75
Changes in Resources or
Right-Hand-Side Values
◼ If the right-hand side of a constraint is
changed, the feasible region will change
(unless the constraint is redundant) and
often the optimal solution will change
◼ The amount of change (increase or
decrease) in the objective function value
that results from a unit change in one of
the resources available is called the dual
price or dual value
© 2009 Prentice-Hall, Inc. 7 – 76
Changes in Resources or
Right-Hand-Side Values
◼ However, the amount of possible increase in the
right-hand side of a resource is limited
◼ If the number of hours increases beyond the
upper bound (or decreases below the lower
bound), then the objective function would no
longer increase (decrease) by the dual price.
◼ There may be excess (slack) hours of a resource
or the objective function may change by an
amount different from the dual price.
◼ Thus, the dual price is relevant only within limits.
◼ If the dual value of a constraint is zero
◼ The slack is positive, indicating unused resource
◼ Additional amount of resource will simply increase
the amount of slack.
◼ The upper limit of infinity indicates that adding
more hours would simply increase the amount of
slack.
© 2009 Prentice-Hall, Inc. 7 – 77
Changes in the Electrician's Time
for High Note Sound
X2 (a)
60 – If the electricians’ hours are changed from 80 to
100, the new optimal solution is (0,25) with profit
of $3,000. The extra 20 hours resulted in an
increase in profit of $600 or $30/hour
40 –
Constraint Representing 60 Hours of
Audio Technician’s Time Resource
a
25 –
20 – b Changed Constraint Representing 100 Hours
of Electrician’s Time Resource
– | c | | |
0 20 40 50 60 X1
Figure 7.19
© 2009 Prentice-Hall, Inc. 7 – 78
Changes in the Electrician's Time
for High Note Sound
If the hours are decreased to 60,
X2 (b)
the new solution is (0,15) and the
60 –
profit is $1,800. Reducing hours
by 20 results in a $600 decrease in
profit or $30/hour
40 –
Constraint Representing 60 Hours of
Audio Technician’s Time Resource
20 – Changed Constraint Representing 60 Hours
a of Electrician’s Time Resource
15 –
b
– c | | | |
0 20 30 40 60 X1
Figure 7.19
© 2009 Prentice-Hall, Inc. 7 – 79
Changes in the Electrician's Time
for High Note Sound
X2 (c)
60 –
Changed Constraint Representing 240 Hours
of Electrician’s Time Resource
40 –
Constraint
Representing
20 – 60 Hours of Audio
Technician’s
Time Resource
– | | | | | |
0 20 40 60 80 100 120
X1
Figure 7.19
© 2009 Prentice-Hall, Inc. 7 – 80
Changes in the Electrician's Time
for High Note Sound
◼ If total electrician time was increased to 240,
the optimal solution would be (0,60) with a
profit of $7,200. This is $2,400 (the original
solution) + $30 (dual price)*160 hours(240-80)
◼ If the hours increases beyond 240, then the
optimal solution would still be (0,60) and profit
would not increase. The extra time is slack during
which the electricians are not working
© 2009 Prentice-Hall, Inc. 7 – 81
QM for Windows and Changes in
Right-Hand-Side Values
◼ Input and sensitivity analysis for High Note Sound
data
Program 7.5A
Program 7.5B
© 2009 Prentice-Hall, Inc. 7 – 82
Flair Furniture – Sensitivity
Analysis
Dual/Value RHS change Solution change
Carpentry/1.5 240→241 30/40: 410→
29.5/41: 411.5
Painting/.5 100→101 30/40: 410→
31.5/38: 410.5
Profit per item Solution Range
Table 7 6.67 → 10
Chair 5 3.5→ 5.25
© 2009 Prentice-Hall, Inc. 7 – 83
Personal Mini Warehouses
1. For the optimal solution, how much is spent on
advertising?
2. For the optimal solution, how much square footage
will be used?
3. Would the solution change if the advertising budget
were only $300 instead of $400? Why?
4. What would the optimal solution be if the profit on
the large spaces were reduced from $50 to $45?
5. How much would earnings increase if the square
footage requirement were increased from 8,000 to
9,000?
1. 2*60 + 4*40 = 280 (note slack of 120 on advertising constraint)
2. 100*60 + 50*40 = 8,000 (no slack)
3. Advertising has a lower bound of 280 and upper bound of Infinity. 300 is within the
bounds so the solution does not change.
4. The profit on large spaces has a range of 40 to Infinity, 45 is within that range so the
solution would not change. However, earnings would be reduced from $3,800 to 45*60 +
20*40 = $3,500.
5. The dual price for square footage is .4 and the upper bound is 9,500. By increasing the
square footage by 1,000 we increase the earning by .4*1,000 = $400
© 2009 Prentice-Hall, Inc. 7 – 84
Chapter 9
Linear Programming:
The Simplex Method
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Learning Objectives
After completing this chapter, students will be able to:
1. Convert LP constraints to equalities with slack,
surplus, and artificial variables
2. Set up and solve LP problems with simplex
tableaus
3. Interpret the meaning of every number in a
simplex tableau
4. Recognize special cases such as infeasibility,
unboundedness, and degeneracy
5. Use the simplex tables to conduct sensitivity
analysis
6. Construct the dual problem from the primal
problem
© 2009 Prentice-Hall, Inc. 9–2
Chapter Outline
9.1 Introduction
9.2 How to Set Up the Initial Simplex
Solution
9.3 Simplex Solution Procedures
9.4 The Second Simplex Tableau
9.5 Developing the Third Tableau
9.6 Review of Procedures for Solving
LP Maximization Problems
9.7 Surplus and Artificial Variables
© 2009 Prentice-Hall, Inc. 9–3
Chapter Outline
9.8 Solving Minimization Problems
9.9 Review of Procedures for Solving
LP Minimization Problems
9.10 Special Cases
9.11 Sensitivity Analysis with the
Simplex Tableau
9.12 The Dual
9.13 Karmarkar’s Algorithm
© 2009 Prentice-Hall, Inc. 9–4
Introduction
With only two decision variables it is possible to
use graphical methods to solve LP problems
But most real life LP problems are too complex for
simple graphical procedures
We need a more powerful procedure called the
simplex method
The simplex method examines the corner points
in a systematic fashion using basic algebraic
concepts
It does this in an iterative manner until an optimal
solution is found
Each iteration moves us closer to the optimal
solution
© 2009 Prentice-Hall, Inc. 9–5
Introduction
Why should we study the simplex method?
It is important to understand the ideas used to
produce solutions
It provides the optimal solution to the decision
variables and the maximum profit (or minimum
cost)
It also provides important economic information
To be able to use computers successfully and to
interpret LP computer printouts, we need to know
what the simplex method is doing and why
© 2009 Prentice-Hall, Inc. 9–6
How To Set Up The Initial
Simplex Solution
Let’s look at the Flair Furniture Company from
Chapter 7
This time we’ll use the simplex method to solve
the problem
You may recall
T = number of tables produced
C = number of chairs produced
and
Maximize profit = $70T + $50C (objective function)
subject to 2T + 1C ≤ 100 (painting hours constraint)
4T + 3C ≤ 240 (carpentry hours constraint)
T, C ≥ 0 (nonnegativity constraint)
© 2009 Prentice-Hall, Inc. 9–7
Converting the Constraints
to Equations
The inequality constraints must be converted into
equations
Less-than-or-equal-to constraints (≤) are
converted to equations by adding a slack variable
to each
Slack variables represent unused resources
For the Flair Furniture problem, the slacks are
S1 = slack variable representing unused hours
in the painting department
S2 = slack variable representing unused hours
in the carpentry department
The constraints may now be written as
2T + 1C + S1 = 100
4T + 3C + S2 = 240
© 2009 Prentice-Hall, Inc. 9–8
Converting the Constraints
to Equations
If the optimal solution uses less than the
available amount of a resource, the unused
resource is slack
For example, if Flair produces T = 40 tables and
C = 10 chairs, the painting constraint will be
2T + 1C + S1 = 100
2(40) +1(10) + S1 = 100
S1 = 10
There will be 10 hours of slack, or unused
painting capacity
© 2009 Prentice-Hall, Inc. 9–9
Converting the Constraints
to Equations
Each slack variable must appear in every
constraint equation
Slack variables not actually needed for an
equation have a coefficient of 0
So
2T + 1C + 1S1 + 0S2 = 100
4T + 3C +0S1 + 1S2 = 240
T, C, S1, S2 ≥ 0
The objective function becomes
Maximize profit = $70T + $50C + $0S1 + $0S2
© 2009 Prentice-Hall, Inc. 9 – 10
Finding an Initial Solution
Algebraically
There are now two equations and four
variables
When there are more unknowns than
equations, you have to set some of the
variables equal to 0 and solve for the
others
In this example, two variables must be set
to 0 so we can solve for the other two
A solution found in this manner is called a
basic feasible solution
© 2009 Prentice-Hall, Inc. 9 – 11
Finding an Initial Solution
Algebraically
The simplex method starts with an initial feasible
solution where all real variables are set to 0
While this is not an exciting solution, it is a corner
point solution
Starting from this point, the simplex method will
move to the corner point that yields the most
improved profit
It repeats the process until it can further improve
the solution
On the following graph, the simplex method starts
at point A and then moves to B and finally to C,
the optimal solution
© 2009 Prentice-Hall, Inc. 9 – 12
Finding an Initial Solution
Algebraically
Corner points C
for the Flair
100 –
Furniture
–
Company B = (0, 80)
Number of Chairs
80 –
problem 2T + 1C ≤ 100
–
60 –
–
40 – C = (30, 40)
–
20 – 4T + 3C ≤ 240
– D = (50, 0)
(0, 0) A |– | | | |
0 20 40 60 80 T
Figure 9.1 Number of Tables
© 2009 Prentice-Hall, Inc. 9 – 13
The First Simplex Tableau
Constraint equations
It simplifies handling the LP equations if we
put them in tabular form
These are the constraint equations for the Flair
Furniture problem
QUANTITY
SOLUTION MIX T C S1 S2 (RIGHT-HAND SIDE)
S1 2 1 1 0 100
S2 4 3 0 1 240
© 2009 Prentice-Hall, Inc. 9 – 14
The First Simplex Tableau
The first tableau is is called a simplex tableau
ix
es
s
it
le
bl
un
ns iab
n on
ns ria
n r
n t
m ti
m e
m n
m a
m r
lu it p
lu l va
lu uc
lu k v
lu sta
co rod
co lac
co rof
co on
co ea
R
C
P
P
S
Cj SOLUTION $70 $50 $0 $0 Profit per
MIX T C S1 S2 QUANTITY unit row
$0 S1 2 1 1 0 100 Constraint
equation rows
$0 S2 4 3 0 1 240
Gross
Zj $0 $0 $0 $0 $0 profit row
Cj - Zj $70 $50 $0 $0 $0
Net profit row
Table 9.1
© 2009 Prentice-Hall, Inc. 9 – 15
The First Simplex Tableau
The numbers in the first row represent the
coefficients in the first constraint and the
numbers in the second the second constraint
At the initial solution, T = 0 and C = 0, so S1 = 100
and S2 = 240
The two slack variables are the initial solution
mix
The values are found in the QUANTITY column
The initial solution is a basic feasible solution
T 0
C 0
=
S1 100
S2 240
© 2009 Prentice-Hall, Inc. 9 – 16
The First Simplex Tableau
Variables in the solution mix, called the basis in
LP terminology, are referred to as basic variables
Variables not in the solution mix or basis (value
of 0) are called nonbasic variables
The optimal solution was T = 30, C = 40, S1 = 0,
and S2 = 0
The final basic variables would be
T 30
C 40
=
S1 0
S2 0
© 2009 Prentice-Hall, Inc. 9 – 17
The First Simplex Tableau
Substitution rates
The numbers in the body of the tableau are the
coefficients of the constraint equations
These can also be thought of as substitution
rates
Using the variable T as an example, if Flair
were to produce 1 table (T = 1), 2 units of S1
and 4 units of S2 would have to be removed
from the solution
Similarly, the substitution rates for C are 1 unit
of S1 and 3 units of S2
Also, for a variable to appear in the solution
mix, it must have a 1 someplace in its column
and 0s in every other place in that column
© 2009 Prentice-Hall, Inc. 9 – 18
The First Simplex Tableau
Adding the objective function
We add a row to the tableau to reflect the
objective function values for each variable
These contribution rates are called Cj and
appear just above each respective variable
In the leftmost column, Cj indicates the unit
profit for each variable currently in the solution
mix
Cj $70 $50 $0 $0
SOLUTION
S1 S2 QUANTITY
MIX T C
$0 S1 2 1 1 0 100
$0 S2 4 3 0 1 240
© 2009 Prentice-Hall, Inc. 9 – 19
The First Simplex Tableau
The Zj and Cj – Zj rows
We can complete the initial tableau by adding
two final rows
These rows provide important economic
information including total profit and whether
the current solution is optimal
We compute the Zj value by multiplying the
contribution value of each number in a column
by each number in that row and the jth column,
and summing
© 2009 Prentice-Hall, Inc. 9 – 20
The First Simplex Tableau
The Zj value for the quantity column provides the
total contribution of the given solution
Zj (gross profit) = (Profit per unit of S1) × (Number of units of S1)
+ (profit per unit of S2) × (Number of units of S2)
= $0 × 100 units + $0 × 240 units
= $0 profit
The Zj values in the other columns represent the
gross profit given up by adding one unit of this
variable into the current solution
Zj = (Profit per unit of S1) × (Substitution rate in row 1)
+ (profit per unit of S2) × (Substitution rate in row 2)
© 2009 Prentice-Hall, Inc. 9 – 21
The First Simplex Tableau
Thus,
Zj (for column T) = ($0)(2) + ($0)(4) = $0
Zj (for column C) = ($0)(1) + ($0)(3) = $0
Zj (for column S1) = ($0)(1) + ($0)(0) = $0
Zj (for column S2) = ($0)(0) + ($0)(1) = $0
We can see that no profit is lost by adding one
unit of either T (tables), C (chairs), S1, or S2
© 2009 Prentice-Hall, Inc. 9 – 22
The First Simplex Tableau
The Cj – Zj number in each column represents the
net profit that will result from introducing 1 unit of
each product or variable into the solution
It is computed by subtracting the Zj total for each
column from the Cj value at the very top of that
variable’s column
COLUMN
T C S1 S2
Cj for column $70 $50
$0 $0
Zj for column 0 0
0 0
Cj – Zj for column $70 $50
$0 $0
© 2009 Prentice-Hall, Inc. 9 – 23
The First Simplex Tableau
Obviously with a profit of $0, the initial solution is
not optimal
By examining the numbers in the Cj – Zj row in
Table 9.1, we can see that the total profits can be
increased by $70 for each unit of T and $50 for
each unit of C
A negative number in the number in the Cj – Zj row
would tell us that the profits would decrease if the
corresponding variable were added to the
solution mix
An optimal solution is reached when there are no
positive numbers in the Cj – Zj row
© 2009 Prentice-Hall, Inc. 9 – 24
Simplex Solution Procedures
After an initial tableau has been
completed, we proceed through a series of
five steps to compute all the numbers
needed in the next tableau
The calculations are not difficult, but they
are complex enough that even the
smallest arithmetic error can produce a
wrong answer
© 2009 Prentice-Hall, Inc. 9 – 25
Five Steps of the Simplex Method for
Maximization Problems
1. Determine the variable to enter the solution mix
next. One way of doing this is by identifying the
column, and hence the variable, with the largest
positive number in the Cj - Zj row of the preceding
tableau. The column identified in this step is
called the pivot column.
column
2. Determine which variable to replace. This is
accomplished by dividing the quantity column by
the corresponding number in the column selected
in step 1. The row with the smallest nonnegative
number calculated in this fashion will be replaced
in the next tableau. This row is often referred to as
the pivot row.
row The number at the intersection of
the pivot row and pivot column is the pivot
number.
number © 2009 Prentice-Hall, Inc. 9 – 26
Five Steps of the Simplex Method for
Maximization Problems
3. Compute new values for the pivot row. To do this,
we simply divide every number in the row by the
pivot column.
4. Compute the new values for each remaining row.
All remaining rows are calculated as follows:
(New row numbers) = (Numbers in old row)
Number above Corresponding number in
– or below x the new row, that is, the
pivot number row replaced in step 3
© 2009 Prentice-Hall, Inc. 9 – 27
Five Steps of the Simplex Method for
Maximization Problems
5. Compute the Zj and Cj - Zj rows, as demonstrated
in the initial tableau. If all the numbers in the Cj - Zj
row are 0 or negative, an optimal solution has
been reached. If this is not the case, return to step
1.
© 2009 Prentice-Hall, Inc. 9 – 28
The Second Simplex Tableau
We can now apply these steps to the Flair
Furniture problem
Step 1.1 Select the variable with the largest positive
Cj - Zj value to enter the solution next. In this case,
variable T with a contribution value of $70.
Cj $70 $50 $0 $0
SOLUTION QUANTITY
MIX T C S1 S2 (RHS)
$0 S1 2 1 1 0 100
$0 S2 4 3 0 1 240
Zj $0 $0 $0 $0 $0
Cj - Zj $70 $50 $0 $0 total profit
Pivot column
Table 9.2
© 2009 Prentice-Hall, Inc. 9 – 29
The Second Simplex Tableau
Step 2.
2 Select the variable to be replaced. Either S1 or
S2 will have to leave to make room for T in the basis.
The following ratios need to be calculated.
For the S1 row
100(hours of painting time available)
= 50 tables
2(hours required per table)
For the S2 row
240(hours of carpentry time available)
= 60 tables
4(hours required per table)
© 2009 Prentice-Hall, Inc. 9 – 30
The Second Simplex Tableau
We choose the smaller ratio (50) and this determines
the S1 variable is to be replaced. This corresponds to
point D on the graph in Figure 9.2.
Cj $70 $50 $0 $0
SOLUTION QUANTITY
MIX T C S1 S2 (RHS)
$0 S1 2 1 1 0 100
$0 S2 4 3 0 1 240
Pivot number Pivot row
Zj $0 $0 $0 $0 $0
Cj - Zj $70 $50 $0 $0
Pivot column
Table 9.3
© 2009 Prentice-Hall, Inc. 9 – 31
The Second Simplex Tableau
Step 3.
3 We can now begin to develop the second,
improved simplex tableau. We have to compute a
replacement for the pivot row. This is done by
dividing every number in the pivot row by the pivot
number. The new version of the pivot row is below.
2 1 1* 0 100
=1 = 0.5 = 0.5 =0 = 50
2 2 2 2 2
Cj SOLUTION MIX T C S1 S2 QUANTITY
$70 T 1 0.5 0.5 0 50
© 2009 Prentice-Hall, Inc. 9 – 32
The Second Simplex Tableau
Step 4.
4 Completing the rest of the tableau, the S2
row, is slightly more complicated. The right of the
following expression is used to find the left side.
Number in Number in
– Number Below Corresponding Number
New S2 Row = Old S Row Pivot Number × in the New T Row
2
0 = 4 – (4) × (1)
1 = 3 – (4) × (0.5)
–2 = 0 – (4) × (0.5)
1 = 1 – (4) × (0)
40 = 240 – (4) × (50)
Cj SOLUTION MIX T C S1 S2 QUANTITY
$70 T 1 0.5 0.5 0 50
$0 S2 0 1 –2 1 40
© 2009 Prentice-Hall, Inc. 9 – 33
The Second Simplex Tableau
1
The T column contains and the S2 column
0
0
contains , necessary conditions for variables to
1
be in the solution. The manipulations of steps 3 and
4 were designed to produce 0s and 1s in the
appropriate positions.
© 2009 Prentice-Hall, Inc. 9 – 34
The Second Simplex Tableau
Step 5.
5 The final step of the second iteration is to
introduce the effect of the objective function. This
involves computing the Cj - Zj rows. The Zj for the
quantity row gives us the gross profit and the other
Zj represent the gross profit given up by adding one
unit of each variable into the solution.
Zj (for T column) = ($70)(1) + ($0)(0) = $70
Zj (for C column) = ($70)(0.5) + ($0)(1) = $35
Zj (for S1 column) = ($70)(0.5) + ($0)(–2) = $35
Zj (for S2 column) = ($70)(0) + ($0)(1) = $0
Zj (for total profit) = ($70)(50) + ($0)(40) = $3,500
© 2009 Prentice-Hall, Inc. 9 – 35
The Second Simplex Tableau
COLUMN
T C S1 S2
Cj for column $70 $50
$0 $0
Zj for column $70 $35
$35 $0
Cj – Zj for column
Completed $0 $15
second simplex tableau
–$35 $0
Cj $70 $50 $0 $0
SOLUTION QUANTITY
MIX T C S1 S2 (RHS)
$0 T 1 0.5 0.5 0 50
$0 S2 0 1 –2 1 40
Zj $70 $35 $35 $0 $3,500
Cj - Zj $0 $15 –$35 $0
Table 9.4
© 2009 Prentice-Hall, Inc. 9 – 36
Interpreting the Second Tableau
Current solution
The solution point of 50 tables and 0 chairs
(T = 50, C = 0) generates a profit of $3,500. T is
a basic variable and C is a nonbasic variable.
This corresponds to point D in Figure 9.2.
Resource information
Slack variable S2 is the unused time in the
carpentry department and is in the basis. Its
value implies there is 40 hours of unused
carpentry time remaining. Slack variable S1 is
nonbasic and has a value of 0 meaning there is
no slack time in the painting department.
© 2009 Prentice-Hall, Inc. 9 – 37
Interpreting the Second Tableau
Substitution rates
Substitution rates are the coefficients in the
heart of the tableau. In column C, if 1 unit of C
is added to the current solution, 0.5 units of T
and 1 unit of S2 must be given up. This is
because the solution T = 50 uses up all 100
hours of painting time available.
Because these are marginal rates of
substitution, so only 1 more unit of S2 is
needed to produce 1 chair
In column S1, the substitution rates mean that if
1 hour of slack painting time is added to
producing a chair, 0.5 less of a table will be
produced
© 2009 Prentice-Hall, Inc. 9 – 38
Interpreting the Second Tableau
Net profit row
The Cj - Zj row is important for two reasons
First, it indicates whether the current solution
is optimal
When there are no positive values in the
bottom row, an optimal solution to a
maximization LP has been reached
The second reason is that we use this row to
determine which variable will enter the
solution next
© 2009 Prentice-Hall, Inc. 9 – 39
Developing the Third Tableau
Since the previous tableau is not optimal, we
repeat the five simplex steps
Step 1.
1 Variable C will enter the solution as its Cj - Zj
value of 15 is the largest positive value. The C
column is the new pivot column.
Step 2.
2 Identify the pivot row by dividing the number
in the quantity column by its corresponding
substitution rate in the C column.
50
For the T row : = 100 chairs
0.5
40
For the S2 row : = 40 chairs
1
© 2009 Prentice-Hall, Inc. 9 – 40
Developing the Third Tableau
These ratios correspond to the values of C at points
F and C in Figure 9.2. The S2 row has the smallest
ratio so S2 will leave the basis and will be replaced
by C.
Cj $70 $50 $0 $0
SOLUTION
MIX T C S1 S2 QUANTITY
$70 T 1 0.5 0.5 0 50
$0 S2 0 1 –2 1 40
Pivot number Pivot row
Zj $70 $35 $35 $0 $3,500
Cj - Zj $0 $15 –$35 $0
Pivot column
Table 9.5
© 2009 Prentice-Hall, Inc. 9 – 41
Developing the Third Tableau
Step 3.
3 The pivot row is replaced by dividing every
number in it by the pivot point number
0 1 −2 1 40
=0 =1 = −2 =1 = 40
1 1 1 1 1
The new C row is
Cj SOLUTION MIX T C S1 S2 QUANTITY
$5 C 0 1 –2 1 40
© 2009 Prentice-Hall, Inc. 9 – 42
Developing the Third Tableau
Step 4.
4 The new values for the T row may now be
computed
Number in = Number in Number above Corresponding number
new T row old T row – pivot number × in new C row
1 = 1 – (0.5) × (0)
0 = 0.5 – (0.5) × (1)
1.5 = 0.5 – (0.5) × (–2)
–0.5 = 0 – (0.5) × (1)
30 = 50 – (0.5) × (40)
Cj SOLUTION MIX T C S1 S2 QUANTITY
$70 T 1 0 1.5 –0.5 30
$50 C 0 1 –2 1 40
© 2009 Prentice-Hall, Inc. 9 – 43
Developing the Third Tableau
Step 5.
5 The Zj and Cj - Zj rows can now be calculated
Zj (for T column) = ($70)(1) + ($50)(0) = $70
Zj (for C column) = ($70)(0) + ($50)(1) = $50
Zj (for S1 column) = ($70)(1.5) + ($50)(–2)= $5
Zj (for S2 column) = ($70)(–0.5) + ($50)(1)= $15
Zj (for total profit) = ($70)(30) + ($50)(40) = $4,100
And the net profit per unit row is now
COLUMN
T C S1 S2
Cj for column $70 $50
$0 $0
Zj for column $70 $50
$5 $15
Cj – Zj for column $0 $0 –
© 2009 Prentice-Hall, Inc. 9 – 44
Developing the Third Tableau
Note that every number in the Cj - Zj row is 0 or
negative indicating an optimal solution has been
reached
The optimal solution is
T = 30 tables
C = 40 chairs
S1 = 0 slack hours in the painting department
S2 = 0 slack hours in the carpentry department
profit = $4,100 for the optimal solution
© 2009 Prentice-Hall, Inc. 9 – 45
Developing the Third Tableau
The final simplex tableau for the Flair Furniture
problem corresponds to point C in Figure 9.2
Cj $70 $50 $0 $0
SOLUTION
MIX T C S1 S2 QUANTITY
$70 T 1 0 1.5 –0.5 30
$50 C 0 1 –2 1 40
Zj $70 $50 $5 $15 $4,100
Cj - Zj $0 $0 –$5 –$15
Table 9.6
Arithmetic mistakes are easy to make
It is always a good idea to check your answer by going
back to the original constraints and objective function
© 2009 Prentice-Hall, Inc. 9 – 46
Review of Procedures for Solving
LP Maximization Problems
I. Formulate the LP problem’s objective function
and constraints
II. Add slack variables to each less-than-or-equal-
to constraint and to the objective function
III. Develop and initial simplex tableau with slack
variables in the basis and decision variables set
equal to 0. compute the Zj and Cj - Zj values for
this tableau.
IV. Follow the five steps until an optimal solution
has been reached
© 2009 Prentice-Hall, Inc. 9 – 47
Review of Procedures for Solving
LP Maximization Problems
1. Choose the variable with the greatest positive
Cj - Zj to enter the solution in the pivot column.
2. Determine the solution mix variable to be
replaced and the pivot row by selecting the row
with the smallest (nonnegative) ratio of the
quantity-to-pivot column substitution rate.
3. Calculate the new values for the pivot row
4. Calculate the new values for the other row(s)
5. Calculate the Zj and Cj - Zj values for this
tableau. If there are any Cj - Zj numbers greater
than 0, return to step 1. If not, and optimal
solution has been reached.
© 2009 Prentice-Hall, Inc. 9 – 48
Surplus and Artificial Variables
Greater-than-or-equal-to (≥) constraints are just as
common in real problems as less-than-or-equal-to
(≤) constraints and equalities
To use the simplex method with these constraints,
they must be converted to a special form similar
to that made for the less-than-or-equal-to (≤)
constraints
If they are not, the simplex technique is unable to
set up an initial solution in the first tableau
Consider the following two constraints
Constraint 1: 5X1 + 10X2 + 8X3 ≥ 210
Constraint 2: 25X1 + 30X2 = 900
© 2009 Prentice-Hall, Inc. 9 – 49
Surplus and Artificial Variables
Surplus variables
Greater-than-or-equal-to (≥) constraints
require a different approach than the less-
than-or-equal-to (≤) constraints we have seen
They involve the subtraction of a surplus
variable rather than the addition of a slack
variable
The surplus variable tells us how much the
solution exceeds the constraint amount
This is sometimes called negative slack
© 2009 Prentice-Hall, Inc. 9 – 50
Surplus and Artificial Variables
To convert the first constraint we subtract a
surplus variable, S1, to create an equality
Constraint 1 rewritten : 5 X 1 + 10 X 2 + 8 X 3 − S1 = 210
If we solved this for X1 = 20, X2 = 8, X3 = 5, S1 would
be
5 X 1 + 10 X 2 + 8 X 3 − S1 = 210
5(20) + 10(8) + 8(5) − S1 = 210
100 + 80 + 40 − S1 = 210
− S1 = 210 − 220
S1 = 10 surplus units
© 2009 Prentice-Hall, Inc. 9 – 51
Surplus and Artificial Variables
Artificial variables
There is one more step in this process
If a surplus variable is added by itself, it would
have a negative value in the initial tableau
where all real variables are set to zero
5(0) + 10(0) + 8(0) − S1 = 210
0 − S1 = 210
S1 = −210
But all variables in LP problems must be
nonnegative at all times
© 2009 Prentice-Hall, Inc. 9 – 52
Surplus and Artificial Variables
To resolve this we add in another variable called
an artificial variable
Constraint 1 completed : 5 X 1 + 10 X 2 + 8 X 3 − S1 + A1 = 210
Now X1, X2, X3, and S1 can all be 0 in the initial
solution and A1 will equal 210
The same situation applies in equality constraint
equations as well
Constraint 2 rewritten : 25 X 1 + 30 X 2 + A2 = 900
© 2009 Prentice-Hall, Inc. 9 – 53
Surplus and Artificial Variables
Artificial variables are inserted into equality
constraints so we can easily develop an initial
feasible solution
When a problem has many constraint equations
with many variables, it is not possible to “eyeball”
an initial solution
Using artificial variables allows us to use the
automatic initial solution of setting all the other
variables to 0
Unlike slack or surplus variables, artificial
variables have no meaning in the problem
formulation
They are strictly a computational tool, they will be
gone in the final solution
© 2009 Prentice-Hall, Inc. 9 – 54
Surplus and Artificial Variables
Surplus and artificial variables in the
objective function
Both types of variables must be included in
the objective function
Surplus variables, like slack variables, carry a
$0 cost coefficient
Since artificial variables must be forced out of
the solution, we assign an arbitrarily high cost
By convention we use the coefficient M (or –M
in maximization problems) which simply
represents a very large number
© 2009 Prentice-Hall, Inc. 9 – 55
Surplus and Artificial Variables
A problem with this objective function
Minimize cost = $5 X 1 + $9 X 2 + $7 X 3
And the constraint equations we saw before
would appear as follows:
Minimize cost = $5X1 + $9X2 + $7X3 + $0S1 + $MA1 + $MA2
subject to 5X1 + 10X2 + 8X3 – 1S1 + 1A1 + 0A2 = 210
25X1 + 30X2 + 0X3 + 0S1 + 0A1 + 1A2 = 900
© 2009 Prentice-Hall, Inc. 9 – 56
Solving Minimization Problems
Once the necessary equations are
developed for a minimization problem, we
can use the simplex method to solve for
an optimal solution
© 2009 Prentice-Hall, Inc. 9 – 57
The Muddy River Chemical
Corporation Example
The Muddy River Chemical Corporation must
produce exactly 1,000 pounds of a special
mixture of phosphate and potassium for a
customer
Phosphate costs $5 per pound and potassium $6
per pound
No more than 300 pounds of phosphate can be
used and at least 150 pounds of potassium must
be used
The company wants to find the least-cost blend
of the two ingredients
© 2009 Prentice-Hall, Inc. 9 – 58
The Muddy River Chemical
Corporation Example
The model formulation would be
Minimize cost = $5X1 + $6X2
subject to X1 + X2 = 1,000 lb
X1 ≤ 300 lb
X2 ≥ 150 lb
X1 , X2 ≥ 0
where
X1 = number of pounds of phosphate
X2 = number of pounds of potassium
© 2009 Prentice-Hall, Inc. 9 – 59
The Muddy River Chemical
Corporation Example
Graphical analysis
Because there are only two decision variables,
we can plot the constraints and the feasible
region as shown in Figure 9.3
Because X1 + X2 = 1,000 is an equality, the
optimal solution must lie on this line
It must also lie between points A and B
because of the X1 ≤ 300 constraint
It turns out the X2 ≥ 150 is redundant and
nonbinding
The optimal corner point is point B (300, 700)
for a total cost of $5,700
© 2009 Prentice-Hall, Inc. 9 – 60
The Muddy River Chemical
Corporation Example
X2
–
X1 ≤ 300
1,000 – A
800 –
B
600 –
X1 + X2 = 1,000
400 –
X2 ≥ 150
200 – F G H
100 –
E D| |C
0 |– | | |
200 400 600 800 1,000 X1
Figure 9.3
© 2009 Prentice-Hall, Inc. 9 – 61
The Muddy River Chemical
Corporation Example
Rarely will problems be this simple
The simplex method can be used to solve
much more complex problems
In this example, the simplex method will
start at coroner point E, move to point F,
then G and finally to point B which is the
optimal solution
© 2009 Prentice-Hall, Inc. 9 – 62
The Muddy River Chemical
Corporation Example
Converting the constraints and objective
function
The necessary artificial variables, slack
variables, and surplus variables need to be
added to the equations
The revised model is
Minimize cost = $5X1 + $6X2 + $0S1 + $0S2 + $MA1 + $MA2
subject to 1X1 + 1X2 + 0S1 + 0S2 + 1A1 + 0A2 = 1,000
1X1 + 0X2 + 1S1 + 0S2 + 0A1 + 0A2 = 300
0X1 + 1X2 + 0S1 – 1S2 + 0A1 + 1A2 = 150
X1, X2, S1, S2, A1, A2 ≥0
© 2009 Prentice-Hall, Inc. 9 – 63
Rules of the Simplex Method for
Minimization Problems
Minimization problems are quite similar to the
maximization problems tackled earlier
The significant difference is the Cj - Zj row
We will now choose the variable with the negative
Cj - Zj that gives the largest improvement
We select the variable that decreases costs the
most
In minimization problems, an optimal solution is
reached when all the numbers in the Cj - Zj are 0
or positive
All other steps in the simplex method remain the
same
© 2009 Prentice-Hall, Inc. 9 – 64
Steps for Simplex Minimization
Problems
1. Choose the variable with the greatest negative
Cj - Zj to enter the solution in the pivot column.
2. Determine the solution mix variable to be
replaced and the pivot row by selecting the row
with the smallest (nonnegative) ratio of the
quantity-to-pivot column substitution rate.
3. Calculate the new values for the pivot row
4. Calculate the new values for the other row(s)
5. Calculate the Zj and Cj - Zj values for this
tableau. If there are any Cj - Zj numbers less
than 0, return to step 1. if not, and optimal
solution has been reached.
© 2009 Prentice-Hall, Inc. 9 – 65
First Simplex Tableau for the Muddy
River Chemical Corporation Example
The initial tableau is set up in the same manner
as the in the maximization problem
The first three rows are
Note the costs for the artificial variables are $M
We simply treat this as a very large number which
forces the artificial variables out of the solution
quickly
Cj SOLUTION MIX X1 X2 S1 S2 A1 A2 QUANTITY
$M A1 1 1 0 0 1 0 1,000
$0 S1 1 0 1 0 0 0 300
$M A2 0 1 0 –1 0 1 150
© 2009 Prentice-Hall, Inc. 9 – 66
First Simplex Tableau for the Muddy
River Chemical Corporation Example
The numbers in the Zj are computed by
multiplying the Cj column on the far left of the
table times the corresponding numbers in each
other column
Zj (for X1 column) = $M(1) + $0(1) + $M(0) = $M
Zj (for X2 column) = $M(1) + $0(0) + $M(1) = $2M
Zj (for S1 column) = $M(0) + $0(1) + $M(0) = $0
Zj (for S2 column) = $M(0) + $0(0) + $M(–1) = –$M
Zj (for A1 column) = $M(1) + $0(0) + $M(0) = $M
Zj (for A2 column) = $M(0) + $0(0) + $M(1) = $M
Zj (for total cost) = $M(1,000) + $0(300) + $M(150) = $1,150M
© 2009 Prentice-Hall, Inc. 9 – 67
First Simplex Tableau for the Muddy
River Chemical Corporation Example
The Cj – Zj entires are determined as follows
COLUMN
X1 X2 S1 S2 A1 A2
Cj for column $ $
$5 $6 $0 $0 M M
Zj for column $M $2M –$M $ $
$0 M M
Cj – Zj for
–$M + $5 –$2M + $6 $0 $M $0 $0
column
© 2009 Prentice-Hall, Inc. 9 – 68
First Simplex Tableau for the Muddy
River Chemical Corporation Example
The initial solution was obtained by letting each
of the variables X1, X2, and S2 assume a value of 0
The current basic variables are A1 = 1,000, S1 =
150, and A2 = 150
The complete solution could be expressed in
vector form as
X1 0
X2 0
S1 300
=
0
S2
1,000
A1
150
A2
© 2009 Prentice-Hall, Inc. 9 – 69
First Simplex Tableau for the Muddy
River Chemical Corporation Example
The initial tableau
Cj $5 $6 $0 $0 $M $M
SOLUTION X1 X2 S1 S2 A1 A2 QUANTITY
MIX
$M A1 1 1 0 0 1 0 1,000
$0 S1 1 0 1 0 0 0 300
$M A2 0 1 0 –1 0 1 150
Pivot number Pivot row
Zj $M $M $0 –$M $M $M $1,150M
Cj – Zj –$M + $5 –2M + $6 $0 $M $0 $0
Pivot column
Table 9.7
© 2009 Prentice-Hall, Inc. 9 – 70
Developing the Second Tableau
In the Cj – Zj row there are two entries with
negative values, X1 and X2
This means an optimal solution does not yet exist
The negative entry for X2 indicates it has the will
result in the largest improvement, which means it
will enter the solution next
To find the variable that will leave the solution,
we divide the elements in the quantity column by
the respective pivot column substitution rates
© 2009 Prentice-Hall, Inc. 9 – 71
Developing the Second Tableau
1,000
For the A1 row = = 1,000
1
300 (this is an undefined ratio,
For the S1 row = so we ignore it)
0
150 (smallest quotient,
For the A2 row = = 150 indicating pivot row)
1
Hence the pivot row is the A2 row and the pivot
number is at the intersection of the X2 column and
the A2 row
© 2009 Prentice-Hall, Inc. 9 – 72
Developing the Second Tableau
The entering row for the next tableau is found by
dividing each element in the pivot row by the pivot
number
(New row numbers) = (Numbers in old row)
Number above or Corresponding number
– below pivot × in newly replaced row
number
A1 Row S1 Row
1 = 1 – (1)(0) 1 = 1 – (0)(0)
0 = 1 – (1)(1) 0 = 0 – (0)(1)
0 = 0 – (1)(0) 1 = 1 – (0)(0)
1 = 0 – (1)(–1) 0 = 0 – (0)(–1)
1 = 1 – (1)(0) 0 = 0 – (0)(0)
–1 = 0 – (1)(1) 0 = 0 – (0)(1)
850 = 1,000 – (1)(150) 300 = 300 – (0)(150)
© 2009 Prentice-Hall, Inc. 9 – 73
Developing the Second Tableau
The Zj and Cj – Zj rows are computed next
Zj (for X1) = $M(1) + $0(1) + $6(0) = $M
Zj (for X2) = $M(0) + $0(0) + $6(1) = $6
Zj (for S1) = $M(0) + $0(1) + $6(0) = $0
Zj (for S2) = $M(1) + $0(0) + $6(–1) = $M – 6
Zj (for A1) = $M(1) + $0(0) + $6(0) = $M
Zj (for A2) = $M(–1) + $0(0) + $6(1) = –$M + 6
Zj (for total cost) = $M(850) + $0(300) + $6(150) = $850M + 900
COLUMN
X1 X2 S1 S2 A1 A2
Cj for column $5 $6 $0 $0 $M $M
Zj for column –$M +
$M $6 $0 $M – 6 $M
6
Cj – Zj for –$M + $2M –
–$M + $5 $0 $0
column $0 6 6 9 – 74
© 2009 Prentice-Hall, Inc.
Developing the Second Tableau
Second simplex tableau
Cj $5 $6 $0 $0 $M $M
SOLUTION X1 X2 S1 S2 A1 A2
MIX QUANTITY
$M A1 1 0 0 1 1 –1 850
$0 S1 1 0 1 0 0 0 300
Pivot number Pivot row
$6 X2 0 1 0 –1 0 1 150
Zj $850M +
$M $6 $0 $M – 6 $M –$M + 6
$900
Cj – Zj –$M + $5 $0 $0 –$M + $6 $0 $2M – 6
Pivot column
Table 9.8
© 2009 Prentice-Hall, Inc. 9 – 75
Developing a Third Tableau
The new pivot column is the X1 column and we
check the quantity column-to-pivot column ratio
850
For the A1 row = = 850
1
300
For the S1 row = = 300 (smallest ratio)
1
150
For the X 2 row = = undefined
0
Hence variable S1 will be replaced by X1
© 2009 Prentice-Hall, Inc. 9 – 76
Developing a Third Tableau
To replace the pivot row we divide each number in
the S1 row by 1 leaving it unchanged
The other calculations are shown below
A1 Row S1 Row
0 = 1 – (1)(1) 0 = 0 – (0)(1)
0 = 0 – (1)(0) 1 = 1 – (0)(0)
–1 = 0 – (1)(1) 0 = 0 – (0)(1)
1 = 1 – (1)(0) –1 = –1 – (0)(0)
1 = 1 – (1)(0) 0 = 0 – (0)(0)
–1 = –1 – (1)(0) 1 = 1 – (0)(0)
550 = 850 – (1)(300) 150 = 150 – (0)(300)
© 2009 Prentice-Hall, Inc. 9 – 77
Developing a Third Tableau
The Zj and Cj – Zj rows are computed next
Zj (for X1) = $M(0) + $5(1) + $6(0) = $5
Zj (for X2) = $M(0) + $5(0) + $6(1) = $6
Zj (for S1) = $M(–1) + $5(1) + $6(0) = –$M + 5
Zj (for S2) = $M(1) + $5(0) + $6(–1) = $M – 6
Zj (for A1) = $M(1) + $5(0) + $6(0) = $M
Zj (for A2) = $M(–1) + $5(0) + $6(1) = –$M + 6
Zj (for total cost) = $M(550) + $5(300) + $6(150) = $550M + 2,400
COLUMN
X1 X2 S1 S2 A1 A2
Cj for column $5 $6 $0 $0 $M $M
Zj for column –$M +
$5 $6 –$M + 5 $M – 6 $M
6
Cj – Zj for $2M –
$0 $0 $M + 5 –$M + 6 $0
column 6
© 2009 Prentice-Hall, Inc. 9 – 78
Developing a Third Tableau
The third simplex tableau for the Muddy River
Chemical problem
Cj $5 $6 $0 $0 $M $M
SOLUTION X1 X2 S1 S2 A1 A2
MIX QUANTITY
$M A1 0 0 –1 1 1 –1 550
Pivot number Pivot row
$5 X1 1 0 1 0 0 0 300
$6 X2 0 1 0 –1 0 1 150
Zj $5 $6 –$M + 5 $M – 6 $M –$M + 6 $550M + 2,400
Cj – Zj $0 $0 $M – 5 –$M + 6 $0 $2M – 6
Pivot column
Table 9.9
© 2009 Prentice-Hall, Inc. 9 – 79
Fourth Tableau for Muddy River
The new pivot column is the S2 column
550
For the A1 row = = 550 (row to be replaced)
1
300
For the X 1 row = (undefined)
0
150 (not considered
For the X 2 row = because it is
−1
negative)
© 2009 Prentice-Hall, Inc. 9 – 80
Fourth Tableau for Muddy River
Each number in the pivot row is again divided by 1
The other calculations are shown below
X1 Row X2 Row
1 = 1 – (0)(0) 0 = 0 – (–1)(0)
0 = 0 – (0)(0) 1 = 1 – (–1)(0)
1 = 1 – (0)(–1) –1 = 0 – (–1)(–1)
0 = 0 – (0)(1) 0 = –1 – (–1)(1)
0 = 0 – (0)(1) 1 = 0 – (–1)(1)
0 = 0 – (0)(–1) 0 = 1 – (–1)(–1)
300 = 300 – (0)(550) 700 = 150 – (–1)(550)
© 2009 Prentice-Hall, Inc. 9 – 81
Fourth Tableau for Muddy River
Finally the Zj and Cj – Zj rows are computed
Zj (for X1) = $0(0) + $5(1) + $6(0) = $5
Zj (for X2) = $(0) + $5(0) + $6(1) = $6
Zj (for S1) = $0(–1) + $5(1) + $6(–1) = –$1
Zj (for S2) = $0(1) + $5(0) + $6(0) = $0
Zj (for A1) = $0(1) + $5(0) + $6(1) = $6
Zj (for A2) = $0(–1) + $5(0) + $6(0) = $0
Zj (for total cost) = $0(550) + $5(300) + $6(700) = $5,700
COLUMN
X1 X2 S1 S2 A1 A2
Cj for column $5 $6 $0 $0 $M $M
Zj for column $5 $6 –$1 $0 $6 $0
Cj – Zj for $M –
$0 $0 $1 $0 $M
column 6
© 2009 Prentice-Hall, Inc. 9 – 82
Fourth Tableau for Muddy River
Fourth and optimal tableau for the Muddy River
Chemical Corporation problem
Cj $5 $6 $0 $0 $M $M
SOLUTION X1 X2 S1 S2 A1 A2
MIX QUANTITY
$0 S2 0 0 –1 1 1 –1 550
$5 X1 1 0 1 0 0 0 300
$6 X2 0 1 –1 0 1 0 700
Zj $5 $6 –$1 $0 $6 $0 $5,700
Cj – Zj $0 $0 $1 $0 $M – 6 $M
Table 9.10
© 2009 Prentice-Hall, Inc. 9 – 83
Review of Procedures for Solving
LP Minimization Problems
I. Formulate the LP problem’s objective function
and constraints
II. Include slack variables to each less-than-or-
equal-to constraint and both surplus and
artificial variables to greater-than-or-equal-to
constraints and add all variables to the objective
function
III. Develop and initial simplex tableau with artificial
and slack variables in the basis and the other
variables set equal to 0. compute the Zj and
Cj - Zj values for this tableau.
IV. Follow the five steps until an optimal solution
has been reached
© 2009 Prentice-Hall, Inc. 9 – 84
Review of Procedures for Solving
LP Minimization Problems
1. Choose the variable with the negative Cj - Zj
indicating the greatest improvement to enter the
solution in the pivot column
2. Determine the row to be replaced and the pivot
row by selecting the row with the smallest
(nonnegative) quantity-to-pivot column
substitution rate ratio
3. Calculate the new values for the pivot row
4. Calculate the new values for the other row(s)
5. Calculate the Zj and Cj - Zj values for the tableau.
If there are any Cj - Zj numbers less than 0, return
to step 1. If not, and optimal solution has been
reached.
© 2009 Prentice-Hall, Inc. 9 – 85
Special Cases
We have seen how special cases arise
when solving LP problems graphically
They also apply to the simplex method
You remember the four cases are
Infeasibility
Unbounded Solutions
Degeneracy
Multiple Optimal Solutions
© 2009 Prentice-Hall, Inc. 9 – 86
Infeasibility
Infeasibility comes about when there is no
solution that satisfies all of the problem’s
constraints
In the simplex method, an infeasible solution is
indicated by looking at the final tableau
All Cj - Zj row entries will be of the proper sign to
imply optimality, but an artificial variable will still
be in the solution mix
A situation with no feasible solution may exist if
the problem was formulated improperly
© 2009 Prentice-Hall, Inc. 9 – 87
Infeasibility
Illustration of infeasibility
Cj $5 $8 $0 $0 $M $M
SOLUTION X1 X2 S1 S2 A1 A2
MIX QUANTITY
$5 X1 1 0 –2 3 –1 0 200
$8 X2 0 1 1 2 –2 0 100
$M A2 0 0 0 –1 –1 1 20
Zj $5 $8 –$2 $31 – M –$21 – M $M $1,800 + 20M
Cj – Zj $0 $0 $2 $M – 31 $2M + 21 $0
Table 9.11
© 2009 Prentice-Hall, Inc. 9 – 88
Unbounded Solutions
Unboundedness describes linear programs that
do not have finite solutions
It occurs in maximization problems when a
solution variable can be made infinitely large
without violating a constraint
In the simplex method this will be discovered
prior to reaching the final tableau
It will be manifested when trying to decide which
variable to remove from the solution mix
If all the ratios turn out to be negative or
undefined, it indicates that the problem is
unbounded
© 2009 Prentice-Hall, Inc. 9 – 89
Unbounded Solutions
Problem with an unbounded solution
Cj $6 $9 $0 $0
SOLUTION MIX X1 X2 S1 S2 QUANTITY
$9 X2 –1 1 2 0 30
$0 S2 –2 0 –1 1 10
Zj –$9 $9 $18 $0 $270
Cj - Zj $15 $0 –$18 $0
Pivot column
Table 9.12
© 2009 Prentice-Hall, Inc. 9 – 90
Unbounded Solutions
The ratios from the pivot column
30
Ratio for the X 2 row :
−1
Negative ratios
unacceptable
10
Ratio for the S2 row :
−2
Since both pivot column numbers are negative,
an unbounded solution is indicated
© 2009 Prentice-Hall, Inc. 9 – 91
Degeneracy
Degeneracy develops when three constraints
pass through a single point
For example, suppose a problem has only these
three constraints X1 ≤ 10, X2 ≤ 10, and X1 + X2 < 20
All three constraint lines will pass through the
point (10, 10)
Degeneracy is first recognized when the ratio
calculations are made
If there is a tie for the smallest ratio, this is a
signal that degeneracy exists
As a result of this, when the next tableau is
developed, one of the variables in the solution
mix will have a value of zero
© 2009 Prentice-Hall, Inc. 9 – 92
Degeneracy
Degeneracy could lead to a situation known as
cycling in which the simplex algorithm alternates
back and forth between the same nonoptimal
solutions
One simple way of dealing with the issue is to
select either row in question arbitrarily
If unlucky and cycling does occur, simply go
back and select the other row
© 2009 Prentice-Hall, Inc. 9 – 93
Degeneracy
Problem illustrating degeneracy
Cj $5 $8 $2 $0 $0 $0
SOLUTION X1 X2 X3 S1 S2 S3 QUANTITY
MIX
$8 X2 0.25 1 1 –2 0 0 10
$0 S2 4 0 0.33 –1 1 0 20
$0 S3 2 0 2 0.4 0 1 10
Zj $2 $8 $8 $16 $0 $0 $80
Cj - Zj $3 $0 –$6 –$16 $0 $0
Pivot column
Table 9.13
© 2009 Prentice-Hall, Inc. 9 – 94
Degeneracy
The ratios are computed as follows
10
For the X 2 row : = 40
0.25
20
For the S2 row : =5 Tie for the smallest
4 ratio indicates
degeneracy
10
For the S3 row : =5
2
© 2009 Prentice-Hall, Inc. 9 – 95
Multiple Optimal Solutions
In the simplex method, multiple, or alternate,
optimal solutions can be spotted by looking at
the final tableau
If the Cj – Zj value is equal to 0 for a variable that
is not in the solution mix, more than one optimal
solution exists
© 2009 Prentice-Hall, Inc. 9 – 96
Multiple Optimal Solutions
A problem with alternate optimal solutions
Cj $3 $2 $0 $0
SOLUTION MIX X1 X2 S1 S2 QUANTITY
$2 X2 1.5 1 1 0 6
$0 S2 1 0 0.5 1 3
Zj $3 $2 $2 $0 $12
Cj - Zj $0 $0 –$2 $0
Table 9.14
© 2009 Prentice-Hall, Inc. 9 – 97
Sensitivity Analysis with the
Simplex Tableau
Sensitivity analysis shows how the optimal
solution and the value of its objective function
change given changes in various inputs to the
problem
Computer programs handling LP problems of all
sizes provide sensitivity analysis as an important
output feature
Those programs use the information provided in
the final simplex tableau to compute ranges for
the objective function coefficients and ranges for
the RHS values
They also provide “shadow prices,” a concept we
will introduce in this section
© 2009 Prentice-Hall, Inc. 9 – 98
High Note Sound Company Revisited
You will recall the model formulation is
Maximize profit = $50X1 +
$120X2
subject to 2X1 + 4X2 ≤
80
(hours of
electrician time)
And the optimal solution is
3X1 + 1X2 ≤
60
X2 = 20 receivers Basic
(hours of
variablestime)
S2 = 40 hours slack in technician time technician
X1 = 0 CD players Nonbasic
variables
S1 = 0 hours slack in electrician time
© 2009 Prentice-Hall, Inc. 9 – 99
High Note Sound Company Revisited
High Note Sound Company graphical solution
X2
(receivers)
60 –
– Optimal Solution at Point a
40 – X1 = 0 CD Players
X2 = 20 Receivers
–
Profits = $2,400
a = (0, 20)
20 –
b = (16, 12)
10 –
Isoprofit Line: $2,400 = 50X1 + 120X2
0–
| | | | | |
10 20 30 40 50 60 X1
Figure 9.4 c = (20, 0) (CD players)
© 2009 Prentice-Hall, Inc. 9 – 100
Changes in the
Objective Function Coefficient
Optimal solution by the simplex method
Cj $50 $120 $0 $0
SOLUTION X1 X2 S1 S2 QUANTITY
MIX
$120 X2 0.5 1 0.25 0 20
$0 S2 2.5 0 –0.25 1 40
Zj $60 $120 $30 $0 $2,400
Cj - Zj –$10 $0 –$30 $0
Table 9.15
© 2009 Prentice-Hall, Inc. 9 – 101
Changes in the
Objective Function Coefficient
Nonbasic objective function coefficient
The goal is to find out how sensitive the
problem’s optimal solution is to changes in the
contribution rates of variables not currently in
the basis
How much would the objective function
coefficients have to change before X1 or S1
would enter the solution mix and replace one
of the basic variables?
The answer lies in the Cj – Zj row of the final
simplex tableau
© 2009 Prentice-Hall, Inc. 9 – 102
Changes in the
Objective Function Coefficient
This is a maximization problem so the basis will not
change unless the Cj – Zj value of one of the
nonbasic variables becomes greater than 0
The values in the basis will not change as long as Cj
≤ Zj
The solution will not change as long as X1 does not
exceed $60 and the contribution rate of S2 does not
exceed $30
These values can also be made smaller without limit
in this situation
So the range of insignificance for the nonbasic
variables is
− ∞ ≤ C j ( for X 1 ) ≤ $60 − ∞ ≤ C j ( for S1 ) ≤ $30
© 2009 Prentice-Hall, Inc. 9 – 103
Changes in the
Objective Function Coefficient
Basic objective function coefficient
Sensitivity analysis on objective function
coefficients of variables in the basis or
solution mix is slightly more complex
A change in the profit or cost of a basic
variable can affect the Cj – Zj values for all
nonbasic variables
That’s because the Cj value is in both the row
and column
This then impacts the Cj – Zj row
© 2009 Prentice-Hall, Inc. 9 – 104
Changes in the
Objective Function Coefficient
Consider a change in the profit contribution of
stereo receivers
The current coefficient is $120
The changed coefficient will be represented as ∆
The revised final tableau will then be
Cj $50 $120 + ∆ $0 $0
SOLUTION X1 X2 S1 S2
MIX QUANTITY
$120 + ∆ X2 0.5 1 0.25 0 20
$0 S2 2.5 0 –0.25 1 40
Zj $60 + 0.5∆ $120 + ∆ $30 + 0.25∆ $0 $2,400 + 20∆
Cj - Zj –$10 –
$0 –$30 – 0.25∆ $0
0.5∆
Table 9.16
© 2009 Prentice-Hall, Inc. 9 – 105
Changes in the
Objective Function Coefficient
The new Cj – Zj values in the table were
determined in the same way as previous
examples
How may the value of ∆ vary so that all Cj – Zj
entries remain negative?
To find out, solve for ∆ in each column
–10 – 0.5∆ ≤ 0
–10 ≤ 0.5∆
–20 ≤ ∆ or ∆ ≥ –20
This inequality means the optimal solution will not
change unless X2’s profit coefficient decreases by at
least $20, ∆ = –20
© 2009 Prentice-Hall, Inc. 9 – 106
Changes in the
Objective Function Coefficient
Variable X1 will not enter the basis unless the
profit per receiver drops to $100 or less
For the S1 column
–30 – 0.25∆ ≤ 0
–30 ≤ 0.25∆
–120 ≤ ∆ or ∆ ≥ –120
Since the first inequality is more binding, we can
say that the range of optimality for X2’s profit
coefficient is
$100 ≤ C j ( for X 2 ) ≤ ∞
© 2009 Prentice-Hall, Inc. 9 – 107
Changes in the
Objective Function Coefficient
In larger problems, we would use this procedure
to test for the range of optimality of every real
decision variable in the final solution mix
Using this procedure helps us avoid the time-
consuming process of reformulating and
resolving the entire LP problem each time a small
change occurs
Within the bounds, changes in profit coefficients
will not force a change in the optimal solution
The value of the objective function will change,
but this is a comparatively simple calculation
© 2009 Prentice-Hall, Inc. 9 – 108
Changes in Resources or RHS Values
Making changes in the RHS values of
constraints result in changes in the feasible
region and often the optimal solution
Shadow prices
How much should a firm be willing to pay for
one additional unit of a resource?
This is called the shadow price
Shadow pricing provides an important piece of
economic information
This information is available in the final
tableau
© 2009 Prentice-Hall, Inc. 9 – 109
Changes in Resources or RHS Values
Final tableau for High Note Sound
Cj $50 $120 $0 $0
SOLUTION X1 X2 S1 S2 QUANTITY
MIX
$120 X2 0.5 1 0.25 0 20
$0 S2 2.5 0 –0.25 1 40
Zj $60 $120 $30 $0 $2,400
Cj - Zj –$10 $0 –$30 $0
Objective function increases by $30
if 1 additional hour of electricians’
time is made available
Table 9.17
© 2009 Prentice-Hall, Inc. 9 – 110
Changes in Resources or RHS Values
An important property of the Cj – Zj row is that the
negatives of the numbers in its slack variable (Si)
columns provide us with shadow prices
A shadow price is the change in value of the
objective function from an increase of one unit of
a scarce resource
High Note Sound is considering hiring an extra
electrician at $22 per hour
In the final tableau we see S1 (electricians’ time) is
fully utilized and has a Cj – Zj value of –$30
They should hire the electrician as the firm will
net $8 (= $30 – $22)
© 2009 Prentice-Hall, Inc. 9 – 111
Changes in Resources or RHS Values
Should High Note Sound hire a part-time audio
technician at $14 per hour?
In the final tableau we see S2 (audio technician
time) has slack capacity (40 hours) a Cj – Zj value
of $0
Thus there would be no benefit to hiring an
additional audio technician
© 2009 Prentice-Hall, Inc. 9 – 112
Changes in Resources or RHS Values
Right-hand side ranging
We can’t add an unlimited amount of a
resource without eventually violating one of
the other constraints
Right-hand-side ranging tells us how much we
can change the RHS of a scarce resource
without changing the shadow price
Ranging is simple in that it resembles the
simplex process
© 2009 Prentice-Hall, Inc. 9 – 113
Changes in Resources or RHS Values
This table repeats some of the information from
the final tableau for High Note Sound and
includes the ratios
QUANTITY S1 RATIO
20 0.25 20/0.25 = 80
40 –0.25 40/–0.25 = –160
The smallest positive ratio (80 in this example)
tells us how many hours the electricians’ time
can be reduced without altering the current
solution mix
© 2009 Prentice-Hall, Inc. 9 – 114
Changes in Resources or RHS Values
The smallest negative ratio (–160) tells us the
number of hours that can be added to the
resource before the solution mix changes
In this case, that’s 160 hours
So the range over which the shadow price for
electricians’ time is valid is 0 to 240 hours
The audio technician resource is slightly different
There is slack in this resource (S2 = 40) so we can
reduce the amount available by 40 before a
shortage occurs
However, we can increase it indefinitely with no
change in the solution
© 2009 Prentice-Hall, Inc. 9 – 115
Changes in Resources or RHS Values
The substitution rates in the slack variable
column can also be used to determine the actual
values of the solution mix variables if the right-
hand-side of a constraint is changed using the
following relationship
New Original Substitution Change in
quantity = quantity + rate the RHS
© 2009 Prentice-Hall, Inc. 9 – 116
Changes in Resources or RHS Values
For example, if 12 more electrician hours were
made available, the new values in the quantity
column of the simplex tableau are found as
follows
ORIGINAL QUANTITY S1 NEW QUANTITY
20 0.25 20 + 0.25(12) = 23
40 –0.25 40 + (–0.25)(12) = 37
If 12 hours were added, X2 = 23 and S2 = 37
Total profit would be 50(0) + 120(23) = $2,760, an
increase of $360
This of course, is also equal to the shadow price
of $30 times the 12 additional hours
© 2009 Prentice-Hall, Inc. 9 – 117
Sensitivity Analysis by Computer
Solver in Excel has the capability of producing
sensitivity analysis that includes the shadow
prices of resources
The following slides present the solution to the
High Note Sound problem and the sensitivity
report showing shadow prices and ranges
© 2009 Prentice-Hall, Inc. 9 – 118
Sensitivity Analysis by Computer
Program 9.1a
© 2009 Prentice-Hall, Inc. 9 – 119
Sensitivity Analysis by Computer
Program 9.1b
© 2009 Prentice-Hall, Inc. 9 – 120
The Dual
Every LP problem has another LP problem
associated with it called the dual
The first way of stating a problem (what we have
done so far) is called the primal
The second way of stating it is called the dual
The solutions to the primal and dual are
equivalent, but they are derived through
alternative procedures
The dual contains economic information useful to
managers and may be easier to formulate
© 2009 Prentice-Hall, Inc. 9 – 121
The Dual
Generally, if the LP primal is a maximize profit
problem with less-than-or-equal-to resource
constraints, the dual will involve minimizing total
opportunity cost subject to greater-than-or-equal-
to product profit constraints
Formulating a dual problem is not complex and
once formulated, it is solved using the same
procedure as a regular LP problem
© 2009 Prentice-Hall, Inc. 9 – 122
The Dual
Illustrating the primal-dual relationship with the
High Note Sound Company data
The primal problem is to determine the best
production mix between CD players (X1) and
receivers (X2) to maximize profit
Maximize profit = $50X1 +
$120X2
subject to 2X1 + 4X2 ≤
80
(hours of
available
electrician time)
3X1 + 1X2 ≤
60
(hours© 2009
of Prentice-Hall,
audio Inc. 9 – 123
The Dual
The dual of this problem has the objective of
minimizing the opportunity cost of not using the
resources in an optimal manner
The variables in the dual are
U1 = potential hourly contribution of
electrician time, or the dual value of 1
hour of electrician time
U2 = the imputed worth of audio technician
time, or the dual of technician resource
Each constraint in the primal problem will have a
corresponding variable in the dual and each
decision variable in the primal will have a
corresponding constraint in the dual
© 2009 Prentice-Hall, Inc. 9 – 124
The Dual
The RHS quantities of the primal constraints
become the dual’s objective function coefficients
The total opportunity cost will be represented by
the function
Minimize opportunity cost = 80U1 + 60U2
The corresponding dual constraints are formed
from the transpose of the primal constraint
coefficients
2 U1 + 3 U2 ≥ 50 Primal profit coefficients
4 U1 + 1 U2 ≥ 120 Coefficients from the second
primal constraint
Coefficients from the first
primal constraint
© 2009 Prentice-Hall, Inc. 9 – 125
The Dual
The first constraint says that the total imputed
value or potential worth of the scarce resources
needed to produce a CD player must be at least
equal to the profit derived from the product
The second constraint makes an analogous
statement for the stereo receiver product
© 2009 Prentice-Hall, Inc. 9 – 126
Steps to Form the Dual
If the primal is a maximization problem in the
standard form, the dual is a minimization, and
vice versa
The RHS values of the primal constraints become
the dual’s objective coefficients
The primal objective function coefficients
become the RHS values of the dual constraints
The transpose of the primal constraint
coefficients become the dual constraint
coefficients
Constraint inequality signs are reversed
© 2009 Prentice-Hall, Inc. 9 – 127
Solving the Dual of the High Note
Sound Company Problem
The formulation can be restated as
Minimize
opportunity = 80U1 + 60U2 + 0S1 + 0S2 + MA1 + MA2
cost
subject to:
2U1 + 3U2 – 0S1 + 1A1 = 50
4U1 + 1U2 – 0S2 + 1A2 = 120
© 2009 Prentice-Hall, Inc. 9 – 128
Solving the Dual of the High Note
Sound Company Problem
The first and second tableaus
Cj 80 60 0 0 M M
SOLUTION U1 U2 S1 S2 A1 A2 QUANTITY
MIX
First A1
$M 2 3 –1 0 1 0 50
tableau
$M A2 4 1 0 –1 0 1 120
Zj $6M $4M –$M –$M $M $M $170M
Cj – Zj 80 – 6M 60 – 4M M M 0 0
Second U1
$80 1 1.5 –0.5 0 0.5 0 25
tableau
$M A2 0 –5 2 –1 –2 1 20
$120 – $2,000 + 20M
Xj $80 –$40 + 2M –$M $40 – 2M $M
5M
Cj – Xj 0 5M – 60 –2M + 40 M 3M – 40 0
Table 9.18
© 2009 Prentice-Hall, Inc. 9 – 129
Solving the Dual of the High Note
Sound Company Problem
Comparison of the primal and dual optimal
tableaus Primal’s Optimal Solution
Cj $50 $120 $0 $0
Solution Mix X1 X2 S1 S2 Quantity
$120 X2 0.5 1 0.25 0 20
$0 S2 2.5 0 –0.25 1 40
Zj 60 120 30 0 $2,400
Cj – Zj –10 0 –30 0
Dual’s Optimal Solution
Cj 80 60 0 0 M M
Solution Mix U1 U2 S1 S2 A1 A2 Quantity
80 U1 1 0.25 0 –0.25 0 0.5 30
0 S1 0 –2.5 1 –0.5 –1 0.25 10
Zj 80 20 0 –20 0 20 $2,400
Cj – Zj 0 40 0 20 M M – 20
Figure 9.5
© 2009 Prentice-Hall, Inc. 9 – 130
Solving the Dual of the High Note
Sound Company Problem
In the final simplex tableau of a primal problem,
the absolute values of the numbers in the Cj – Zj
row under the slack variables represent the
solutions to the dual problem
They are shadow prices in the primal solution
and marginal profits in the dual
The absolute value of the numbers of the Cj – Zj
values of the slack variables represent the
optimal values of the primal X1 and X2 variables
The maximum opportunity cost derived in the
dual must always equal the maximum profit
derived in the primal
© 2009 Prentice-Hall, Inc. 9 – 131
Karmakar’s Algorithm
In 1984, Narendra Karmakar developed a new
method of solving linear programming problems
called the Karmakar algorithm
The simplex method follows a path of points on
the outside edge of feasible space
Karmakar’s algorithm works by following a path a
points inside the feasible space
It is much more efficient than the simplex method
requiring less computer time to solve problems
It can also handle extremely large problems
allowing organizations to solve previously
unsolvable problems
© 2009 Prentice-Hall, Inc. 9 – 132
Chapter 8
Linear Programming
Applications
To accompany
Quantitative Analysis for Management, Eleventh Edition,
by Render, Stair, and Hanna
Power Point slides created by Brian Peterson
Marketing Applications
Linear programming models have been
used in the advertising field as a decision
aid in selecting an effective media mix.
Media selection problems can be
approached with LP from two perspectives:
Maximize audience exposure.
Minimize advertising costs.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-2
Win Big Gambling Club
The Win Big Gambling Club promotes gambling
junkets to the Bahamas.
It has $8,000 per week to spend on advertising.
Its goal is to reach the largest possible high-
potential audience.
Media types and audience figures are shown in
the following table.
It needs to place at least five radio spots per week.
No more than $1,800 can be spent on radio
advertising each week.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-3
Win Big Gambling Club
Advertising options
AUDIENCE COST PER MAXIMUM ADS
MEDIUM REACHED PER AD AD ($) PER WEEK
TV spot (1 minute) 5,000 800 12
Daily newspaper (full-
8,500 925 5
page ad)
Radio spot (30
2,400 290 25
seconds, prime time)
Radio spot (1 minute,
2,800 380 20
afternoon)
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-4
Win Big Gambling Club
The problem formulation is
X1 = number of 1-minute TV spots each week
X2 = number of daily paper ads each week
X3 = number of 30-second radio spots each week
X4 = number of 1-minute radio spots each week
Objective:
Maximize audience coverage = 5,000X1 + 8,500X2 + 2,400X3 + 2,800X4
Subject to X1 ≤ 12 (max TV spots/wk)
X2 ≤ 5 (max newspaper ads/wk)
X3 ≤ 25 (max 30-sec radio spots ads/wk)
X4 ≤ 20 (max newspaper ads/wk)
800X1 + 925X2 + 290X3 + 380X4 ≤ $8,000 (weekly advertising budget)
X3 + X4 ≥ 5 (min radio spots contracted)
290X3 + 380X4 ≤ $1,800 (max dollars spent on radio)
X 1, X 2, X 3, X 4 ≥ 0
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-5
Win Big Gambling Club
Solution in Excel 2010
Program 8.1
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-6
Manufacturing Applications
Production Mix
LP can be used to plan the optimal mix of
products to manufacture.
Company must meet a myriad of constraints,
ranging from financial concerns to sales
demand to material contracts to union labor
demands.
Its primary goal is to generate the largest profit
possible.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-7
Fifth Avenue Industries
Fifth Avenue Industries produces four varieties of
ties:
One is expensive all-silk
One is all-polyester
Two are polyester and cotton blends
The table on the below shows the cost and
availability of the three materials used in the
production process:
MATERIAL AVAILABLE PER
MATERIAL COST PER YARD ($) MONTH (YARDS)
Silk 24 1,200
Polyester 6 3,000
Cotton 9 1,600
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-8
Fifth Avenue Industries
The firm has contracts with several major
department store chains to supply ties.
Contracts require a minimum number of ties but
may be increased if demand increases.
Fifth Avenue’s goal is to maximize monthly profit
given the following decision variables.
X1 = number of all-silk ties produced per month
X2 = number all-polyester ties
X3 = number of blend 1 polyester-cotton ties
X4 = number of blend 2 silk-cotton ties
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-9
Fifth Avenue Industries Data
MATERIAL
SELLING MONTHLY REQUIRED
VARIETY OF PRICE PER CONTRACT MONTHLY PER TIE MATERIAL
TIE TIE ($) MINIMUM DEMAND (YARDS) REQUIREMENTS
All silk 19.24 5,000 7,000 0.125 100% silk
All polyester 8.70 10,000 14,000 0.08 100% polyester
Poly – cotton 50% polyester –
9.52 13,000 16,000 0.10
blend 1 50% cotton
Silk-cotton 60% silk - 40%
10.64 5,000 8,500 0.11
blend 2 cotton
Table 8.1
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-10
Fifth Avenue Industries
Fifth Avenue also has to calculate profit per tie
for the objective function.
SELLING MATERIAL MATERIAL
PRICE REQUIRED PER COST PER COST PER PROFIT
VARIETY OF TIE PER TIE ($) TIE (YARDS) YARD ($) TIE ($) PER TIE ($)
All silk $19.24 0.125 $24 $3.00 $16.24
All polyester $8.70 0.08 $6 $0.48 $8.22
Poly-cotton
$9.52 0.05 $6 $0.30
blend 1
0.05 $9 $0.45 $8.77
Silk – cotton
$10.64 0.06 $24 $1.44
blend 2
0.06 $9 $0.54 $8.66
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-11
Fifth Avenue Industries
The complete Fifth Avenue Industries Model
Objective function
Maximize profit = $16.24X1 + $8.22X2 + $8.77X3 + $8.66X4
Subject to 0.125X1+ 0.066X4 ≤ 1200 (yds of silk)
0.08X2 + 0.05X3 ≤ 3,000 (yds of polyester)
0.05X3 + 0.44X4 ≤ 1,600 (yds of cotton)
X1 ≥ 5,000 (contract min for silk)
X1 ≤ 7,000 (contract min)
X2 ≥ 10,000 (contract min for all polyester)
X2 ≤ 14,000 (contract max)
X3 ≥ 13,000 (contract mini for blend 1)
X3 ≤ 16,000 (contract max)
X4 ≥ 5,000 (contract mini for blend 2)
X4 ≤ 8,500 (contract max)
X 1, X 2, X 3, X 4 ≥ 0
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-12
Fifth Avenue Solution in Excel
2010
Program 8.3
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-13
Financial Applications
Portfolio Selection
Bank, investment funds, and insurance
companies often have to select specific
investments from a variety of alternatives.
The manager’s overall objective is generally to
maximize the potential return on the
investment given a set of legal, policy, or risk
restraints.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-14
International City Trust
International City Trust (ICT) invests in short-term
trade credits, corporate bonds, gold stocks, and
construction loans.
The board of directors has placed limits on how
much can be invested in each area:
INTEREST MAXIMUM INVESTMENT
INVESTMENT EARNED (%) ($ MILLIONS)
Trade credit 7 1.0
Corporate bonds 11 2.5
Gold stocks 19 1.5
Construction loans 15 1.8
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-15
International City Trust
ICT has $5 million to invest and wants to
accomplish two things:
Maximize the return on investment over the next six
months.
Satisfy the diversification requirements set by the
board.
The board has also decided that at least 55% of
the funds must be invested in gold stocks and
construction loans and no less than 15% be
invested in trade credit.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-16
International City Trust
The variables in the model are:
X1 = dollars invested in trade credit
X2 = dollars invested in corporate bonds
X3 = dollars invested in gold stocks
X4 = dollars invested in construction loans
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-17
International City Trust
Objective:
Maximize
dollars of
interest = 0.07X1 + 0.11X2 + 0.19X3 + 0.15X4
earned
subject to: X1 ≤ 1,000,000
X2 ≤ 2,500,000
X3 ≤ 1,500,000
X4 ≤ 1,800,000
X3 + X 4 ≥ 0.55(X1 + X2 + X3 + X4)
X1 ≥ 0.15(X1 + X2 + X3 + X4)
X1 + X2 + X3 + X 4 ≤ 5,000,000
X 1, X 2, X 3, X 4 ≥ 0
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-18
International City Trust
The optimal solution to the ICT is to make the
following investments:
X1 = $750,000
X2 = $950,000
X3 = $1,500,000
X4 = $1,800,000
The total interest earned with this plan is
$712,000.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-19
ICT Portfolio Solution in Excel
2010
Program 8.6
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 8-20
Chapter 10
Transportation and
Assignment Models
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Learning Objectives
After completing this chapter, students will be able to:
1. Structure special LP problems using the
transportation and assignment models
2. Use the northwest corner, VAM, MODI, and
stepping-stone methods
3. Solve facility location and other application
problems with transportation models
4. Solve assignment problems with the
Hungarian (matrix reduction) method
© 2009 Prentice-Hall, Inc. 10 – 2
Chapter Outline
10.1 Introduction
10.2 Setting Up a Transportation Problem
10.3 Developing an Initial Solution: Northwest
Corner Rule
10.4 Stepping-Stone Method: Finding a
Least-Cost Solution
10.5 MODI Method
10.6 Vogel’s Approximation Method: Another
Way to Find an Initial Solution
10.7 Unbalanced Transportation Problems
© 2009 Prentice-Hall, Inc. 10 – 3
Chapter Outline
10.8 Degeneracy in Transportation Problems
10.9 More Than One Optimal Solution
10.10 Maximization Transportation Problems
10.11 Unacceptable or Prohibited Routes
10.12 Facility Location Analysis
10.13 Assignment Model Approach
10.14 Unbalanced Assignment Problems
10.15 Maximization Assignment Problems
© 2009 Prentice-Hall, Inc. 10 – 4
Introduction
In this chapter we will explore two special
linear programming models
The transportation model
The assignment model
Because of their structure, they can be
solved more efficiently than the simplex
method
These problems are members of a category
of LP techniques called network flow
problems
© 2009 Prentice-Hall, Inc. 10 – 5
Introduction
Transportation model
The transportation problem deals with the
distribution of goods from several points of
supply (sources)
sources to a number of points of
demand (destinations)
destinations
Usually we are given the capacity of goods at
each source and the requirements at each
destination
Typically the objective is to minimize total
transportation and production costs
© 2009 Prentice-Hall, Inc. 10 – 6
Introduction
Example of a transportation problem in a network
format
Factories Warehouses
(Sources) (Destinations)
100 Units Des Moines Albuquerque 300 Units
300 Units Evansville Boston 200 Units
300 Units Fort Lauderdale Cleveland 200 Units
Capacities Shipping Routes Requirements
Figure 10.1
© 2009 Prentice-Hall, Inc. 10 – 7
Introduction
Assignment model
The assignment problem refers to the class of
LP problems that involve determining the most
efficient assignment of resources to tasks
The objective is most often to minimize total
costs or total time to perform the tasks at hand
One important characteristic of assignment
problems is that only one job or worker can be
assigned to one machine or project
© 2009 Prentice-Hall, Inc. 10 – 8
Introduction
Special-purpose algorithms
Although standard LP methods can be used to
solve transportation and assignment problems,
special-purpose algorithms have been
developed that are more efficient
They still involve finding and initial solution and
developing improved solutions until an optimal
solution is reached
They are fairly simple in terms of computation
© 2009 Prentice-Hall, Inc. 10 – 9
Introduction
Streamlined versions of the simplex method are
important for two reasons
1. Their computation times are generally 100 times faster
2. They require less computer memory (and hence can
permit larger problems to be solved)
Two common techniques for developing initial
solutions are the northwest corner method and
Vogel’s approximation
The initial solution is evaluated using either the
stepping-stone method or the modified
distribution (MODI) method
We also introduce a solution procedure called the
Hungarian method,
method Flood’s technique,
technique or the
reduced matrix method
© 2009 Prentice-Hall, Inc. 10 – 10
Setting Up a Transportation Problem
The Executive Furniture Corporation
manufactures office desks at three locations: Des
Moines, Evansville, and Fort Lauderdale
The firm distributes the desks through regional
warehouses located in Boston, Albuquerque, and
Cleveland
Estimates of the monthly production capacity of
each factory and the desks needed at each
warehouse are shown in Figure 10.1
© 2009 Prentice-Hall, Inc. 10 – 11
Setting Up a Transportation Problem
Production costs are the same at the three
factories so the only relevant costs are shipping
from each source to each destination
Costs are constant no matter the quantity
shipped
The transportation problem can be described as
how to select the shipping routes to be used and
the number of desks to be shipped on each route
so as to minimize total transportation cost
Restrictions regarding factory capacities and
warehouse requirements must be observed
© 2009 Prentice-Hall, Inc. 10 – 12
Setting Up a Transportation Problem
The first step is setting up the transportation
table
Its purpose is to summarize all the relevant data
and keep track of algorithm computations
Transportation costs per desk for Executive Furniture
TO
FROM ALBUQUERQUE BOSTON CLEVELAND
DES MOINES $5 $4 $3
EVANSVILLE $8 $4 $3
FORT LAUDERDALE $9 $7 $5
Table 10.1
© 2009 Prentice-Hall, Inc. 10 – 13
Setting Up a Transportation Problem
Geographical locations of Executive Furniture’s
factories and warehouses
Boston
Cleveland
Factory
Des Moines
Evanston Warehouse
Albuquerque
Fort Lauderdale
Figure 10.2
© 2009 Prentice-Hall, Inc. 10 – 14
Setting Up a Transportation Problem
Des Moines
Transportation table for Executive Furniture capacity
constraint
TO WAREHOUSE WAREHOUSE WAREHOUSE
AT AT AT FACTORY
FROM ALBUQUERQUE BOSTON CLEVELAND CAPACITY
DES MOINES $5 $4 $3
100
FACTORY
EVANSVILLE $8 $4 $3
300
FACTORY
FORT LAUDERDALE $9 $7 $5
300
FACTORY
WAREHOUSE
300 200 200 700
REQUIREMENTS
Cell representing a
Table 10.2 Total supply source-to-destination
Cost of shipping 1 unit from Cleveland (Evansville to Cleveland)
and demand
Fort Lauderdale factory to warehouse shipping assignment
Boston warehouse demand that could be made
© 2009 Prentice-Hall, Inc. 10 – 15
Setting Up a Transportation Problem
In this table, total factory supply exactly
equals total warehouse demand
When equal demand and supply occur, a
balanced problem is said to exist
This is uncommon in the real world and
we have techniques to deal with
unbalanced problems
© 2009 Prentice-Hall, Inc. 10 – 16
Developing an Initial Solution:
Northwest Corner Rule
Once we have arranged the data in a table, we
must establish an initial feasible solution
One systematic approach is known as the
northwest corner rule
Start in the upper left-hand cell and allocate units
to shipping routes as follows
1. Exhaust the supply (factory capacity) of each row
before moving down to the next row
2. Exhaust the demand (warehouse) requirements of
each column before moving to the right to the next
column
3. Check that all supply and demand requirements are
met.
In this problem it takes five steps to make the
initial shipping assignments © 2009 Prentice-Hall, Inc. 10 – 17
Developing an Initial Solution:
Northwest Corner Rule
1. Beginning in the upper left hand corner, we
assign 100 units from Des Moines to
Albuquerque. This exhaust the supply from Des
Moines but leaves Albuquerque 200 desks short.
We move to the second row in the same column.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
300
(E)
FORT LAUDERDALE $9 $7 $5
300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
© 2009 Prentice-Hall, Inc. 10 – 18
Developing an Initial Solution:
Northwest Corner Rule
2. Assign 200 units from Evansville to Albuquerque.
This meets Albuquerque’s demand. Evansville
has 100 units remaining so we move to the right
to the next column of the second row.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 300
(E)
FORT LAUDERDALE $9 $7 $5
300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
© 2009 Prentice-Hall, Inc. 10 – 19
Developing an Initial Solution:
Northwest Corner Rule
3. Assign 100 units from Evansville to Boston. The
Evansville supply has now been exhausted but
Boston is still 100 units short. We move down
vertically to the next row in the Boston column.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 100 300
(E)
FORT LAUDERDALE $9 $7 $5
300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
© 2009 Prentice-Hall, Inc. 10 – 20
Developing an Initial Solution:
Northwest Corner Rule
4. Assign 100 units from Fort Lauderdale to Boston.
This fulfills Boston’s demand and Fort
Lauderdale still has 200 units available.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 100 300
(E)
FORT LAUDERDALE $9 $7 $5
100 300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
© 2009 Prentice-Hall, Inc. 10 – 21
Developing an Initial Solution:
Northwest Corner Rule
5. Assign 200 units from Fort Lauderdale to
Cleveland. This exhausts Fort Lauderdale’s
supply and Cleveland’s demand. The initial
shipment schedule is now complete.
Table 10.3
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 100 300
(E)
FORT LAUDERDALE $9 $7 $5
100 200 300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
© 2009 Prentice-Hall, Inc. 10 – 22
Developing an Initial Solution:
Northwest Corner Rule
We can easily compute the cost of this shipping
assignment
ROUTE
FROM TO UNITS PER UNIT TOTAL
SHIPPED x COST ($) = COST ($)
D A 100 5 500
E A 200 8 1,600
E B 100 4 400
F B 100 7 700
F C 200 5 1,000
4,200
This solution is feasible but we need to check to
see if it is optimal
© 2009 Prentice-Hall, Inc. 10 – 23
Stepping-Stone Method:
Finding a Least Cost Solution
The stepping-stone method is an iterative
technique for moving from an initial
feasible solution to an optimal feasible
solution
There are two distinct parts to the process
Testing the current solution to determine if
improvement is possible
Making changes to the current solution to
obtain an improved solution
This process continues until the optimal
solution is reached
© 2009 Prentice-Hall, Inc. 10 – 24
Stepping-Stone Method:
Finding a Least Cost Solution
There is one very important rule
The number of occupied routes (or squares) must
always be equal to one less than the sum of the
number of rows plus the number of columns
In the Executive Furniture problem this means
the initial solution must have 3 + 3 – 1 = 5
squares used
Occupied shipping Number Number of
routes (squares) = of rows + columns –1
When the number of occupied rows is less than
this, the solution is called degenerate
© 2009 Prentice-Hall, Inc. 10 – 25
Testing the Solution for
Possible Improvement
The stepping-stone method works by
testing each unused square in the
transportation table to see what would
happen to total shipping costs if one unit
of the product were tentatively shipped on
an unused route
There are five steps in the process
© 2009 Prentice-Hall, Inc. 10 – 26
Five Steps to Test Unused Squares
with the Stepping-Stone Method
1. Select an unused square to evaluate
2. Beginning at this square, trace a closed path
back to the original square via squares that are
currently being used with only horizontal or
vertical moves allowed
3. Beginning with a plus (+) sign at the unused
square, place alternate minus (–) signs and plus
signs on each corner square of the closed path
just traced
© 2009 Prentice-Hall, Inc. 10 – 27
Five Steps to Test Unused Squares
with the Stepping-Stone Method
4. Calculate an improvement index by adding
together the unit cost figures found in each
square containing a plus sign and then
subtracting the unit costs in each square
containing a minus sign
5. Repeat steps 1 to 4 until an improvement index
has been calculated for all unused squares. If all
indices computed are greater than or equal to
zero, an optimal solution has been reached. If
not, it is possible to improve the current solution
and decrease total shipping costs.
© 2009 Prentice-Hall, Inc. 10 – 28
Five Steps to Test Unused Squares
with the Stepping-Stone Method
For the Executive Furniture Corporation data
Steps 1 and 2.
2 Beginning with Des Moines–Boston
route we trace a closed path using only currently
occupied squares, alternately placing plus and
minus signs in the corners of the path
In a closed path,
path only squares currently used for
shipping can be used in turning corners
Only one closed route is possible for each square
we wish to test
© 2009 Prentice-Hall, Inc. 10 – 29
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Step 3.
3 We want to test the cost-effectiveness of the
Des Moines–Boston shipping route so we pretend
we are shipping one desk from Des Moines to
Boston and put a plus in that box
But if we ship one more unit out of Des Moines
we will be sending out 101 units
Since the Des Moines factory capacity is only
100, we must ship fewer desks from Des Moines
to Albuquerque so we place a minus sign in that
box
But that leaves Albuquerque one unit short so we
must increase the shipment from Evansville to
Albuquerque by one unit and so on until we
complete the entire closed path
© 2009 Prentice-Hall, Inc. 10 – 30
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Evaluating the unused Warehouse A Warehouse B
Des Moines–Boston Factory
100
$5 $4
D
shipping route – +
Factory
+ $8 – $4
E 200 100
TO FACTORY
ALBUQUERQUE BOSTON CLEVELAND
FROM CAPACITY
$5 $4 $3
DES MOINES 100 100
$8 $4 $3
EVANSVILLE 200 100 300
$9 $7 $5
FORT LAUDERDALE 100 200 300
WAREHOUSE
REQUIREMENTS 300 200 200 700 Table 10.4
© 2009 Prentice-Hall, Inc. 10 – 31
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Evaluating the unused Warehouse A Warehouse B
Des Moines–Boston Factory
99
100
$5
1
$4
D
shipping route – +
Factory 201 + $8 99 – $4
E 200 100
TO FACTORY
ALBUQUERQUE BOSTON CLEVELAND
FROM CAPACITY
$5 $4 $3
DES MOINES 100 100
$8 $4 $3
EVANSVILLE 200 100 300
$9 $7 $5
FORT LAUDERDALE 100 200 300
WAREHOUSE
REQUIREMENTS 300 200 200 700 Table 10.4
© 2009 Prentice-Hall, Inc. 10 – 32
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Evaluating the unused Warehouse A Warehouse B
Des Moines–Boston Factory
99
100
$5
1
$4
D
shipping route – +
Factory 201 + $8 99 – $4
E 200 100
TO FACTORY
ALBUQUERQUE BOSTON CLEVELAND
FROM Result
CAPACITY of Proposed
$5 $4 $3
Shift in Allocation
DES MOINES 100 100
= 1 x $4
$8 $4 $3 – 1 x $5
EVANSVILLE 200 100 300 + 1 x $8
– 1 x $4 = +$3
$9 $7 $5
FORT LAUDERDALE 100 200 300
WAREHOUSE
REQUIREMENTS 300 200 200 700 Table 10.4
© 2009 Prentice-Hall, Inc. 10 – 33
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Step 4. 4 We can now compute an improvement index
(Iij) for the Des Moines–Boston route
We add the costs in the squares with plus signs
and subtract the costs in the squares with minus
signs
Des Moines–
Boston index = IDB = +$4 – $5 + $5 – $4 = + $3
This means for every desk shipped via the Des
Moines–Boston route, total transportation cost
will increase by $3 over their current level
© 2009 Prentice-Hall, Inc. 10 – 34
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Step 5.
5 We can now examine the Des Moines–
Cleveland unused route which is slightly more
difficult to draw
Again we can only turn corners at squares that
represent existing routes
We must pass through the Evansville–Cleveland
square but we can not turn there or put a + or –
sign
The closed path we will use is
+ DC – DA + EA – EB + FB – FC
© 2009 Prentice-Hall, Inc. 10 – 35
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Evaluating the Des Moines–Cleveland shipping
route
TO FACTORY
ALBUQUERQUE BOSTON CLEVELAND
FROM CAPACITY
$5 $4 Start $3
DES MOINES 100 100
– +
$8 $4 $3
EVANSVILLE 200 100 300
+ –
$9 $7 $5
FORT LAUDERDALE 100 200 300
+ –
WAREHOUSE 300 200 200 700
REQUIREMENTS
Table 10.5
Des Moines–Cleveland
improvement index = IDC = + $3 – $5 + $8 – $4 + $7 – $5 = + $4
© 2009 Prentice-Hall, Inc. 10 – 36
Five Steps to Test Unused Squares
with the Stepping-Stone Method
Opening the Des Moines–Cleveland route will not
lower our total shipping costs
Evaluating the other two routes we find
Evansville-
Cleveland index = IEC = + $3 – $4 + $7 – $5 = + $1
The closed path is
+ EC – EB + FB – FC
Fort Lauderdale–
Albuquerque index = IFA = + $9 – $7 + $4 – $8 = – $2
The closed path is
+ FA – FB + EB – EA
So opening the Fort Lauderdale-Albuquerque
route will lower our total transportation costs
© 2009 Prentice-Hall, Inc. 10 – 37
Obtaining an Improved Solution
In the Executive Furniture problem there is only
one unused route with a negative index (Fort
Lauderdale-Albuquerque)
If there was more than one route with a negative
index, we would choose the one with the largest
improvement
We now want to ship the maximum allowable
number of units on the new route
The quantity to ship is found by referring to the
closed path of plus and minus signs for the new
route and selecting the smallest number found in
those squares containing minus signs
© 2009 Prentice-Hall, Inc. 10 – 38
Obtaining an Improved Solution
To obtain a new solution, that number is added to
all squares on the closed path with plus signs
and subtracted from all squares the closed path
with minus signs
All other squares are unchanged
In this case, the maximum number that can be
shipped is 100 desks as this is the smallest value
in a box with a negative sign (FB route)
We add 100 units to the FA and EB routes and
subtract 100 from FB and EA routes
This leaves balanced rows and columns and an
improved solution
© 2009 Prentice-Hall, Inc. 10 – 39
Obtaining an Improved Solution
Stepping-stone path used to evaluate route FA
TO FACTORY
A B C
FROM CAPACITY
$5 $4 $3
D 100 100
$8 $4 $3
E 200 100 300
– +
$9 $7 $5
F 100 200 300
+ –
WAREHOUSE 300 200 200 700
REQUIREMENTS
Table 10.6
© 2009 Prentice-Hall, Inc. 10 – 40
Obtaining an Improved Solution
Second solution to the Executive Furniture
problem
TO FACTORY
A B C
FROM CAPACITY
$5 $4 $3
D 100 100
$8 $4 $3
E 100 200 300
$9 $7 $5
F 100 200 300
WAREHOUSE 300 200 200 700
REQUIREMENTS
Table 10.7
Total shipping costs have been reduced by (100
units) x ($2 saved per unit) and now equals $4,000
© 2009 Prentice-Hall, Inc. 10 – 41
Obtaining an Improved Solution
This second solution may or may not be optimal
To determine whether further improvement is
possible, we return to the first five steps to test
each square that is now unused
The four new improvement indices are
D to B = IDB = + $4 – $5 + $8 – $4 = + $3
(closed path: + DB – DA + EA – EB)
D to C = IDC = + $3 – $5 + $9 – $5 = + $2
(closed path: + DC – DA + FA – FC)
E to C = IEC = + $3 – $8 + $9 – $5 = – $1
(closed path: + EC – EA + FA – FC)
F to B = IFB = + $7 – $4 + $8 – $9 = + $2
(closed path: + FB – EB + EA – FA)
© 2009 Prentice-Hall, Inc. 10 – 42
Obtaining an Improved Solution
Path to evaluate for the EC route
TO FACTORY
A B C
FROM CAPACITY
$5 $4 $3
D 100 100
E 100
$8
200
$4 Start $3
300
– +
$9 $7 $5
F 100 200 300
+ –
WAREHOUSE 300 200 200 700
REQUIREMENTS
Table 10.8
An improvement can be made by shipping the
maximum allowable number of units from E to C
© 2009 Prentice-Hall, Inc. 10 – 43
Obtaining an Improved Solution
Total cost of third solution
ROUTE
FROM TO DESKS PER UNIT TOTAL
SHIPPED x COST ($) = COST ($)
D A 100 5 500
E B 200 4 800
E C 100 3 300
F A 200 9 1,800
F C 100 5 500
3,900
© 2009 Prentice-Hall, Inc. 10 – 44
Obtaining an Improved Solution
Third and optimal solution
TO FACTORY
A B C
FROM CAPACITY
$5 $4 $3
D 100 100
$8 $4 $3
E 200 100 300
$9 $7 $5
F 200 100 300
WAREHOUSE 300 200 200 700
REQUIREMENTS
Table 10.9
© 2009 Prentice-Hall, Inc. 10 – 45
Obtaining an Improved Solution
This solution is optimal as the improvement
indices that can be computed are all greater than
or equal to zero
D to B = IDB = + $4 – $5 + $9 – $5 + $3 – $4 = + $2
(closed path: + DB – DA + FA – FC + EC – EB)
D to C = IDC = + $3 – $5 + $9 – $5 = + $2
(closed path: + DC – DA + FA – FC)
E to A = IEA = + $8 – $9 + $5 – $3 = + $1
(closed path: + EA – FA + FC – EC)
F to B = IFB = + $7 – $5 + $3 – $4 = + $1
(closed path: + FB – FC + EC – EB)
© 2009 Prentice-Hall, Inc. 10 – 46
Summary of Steps in Transportation
Algorithm (Minimization)
1. Set up a balanced transportation table
2. Develop initial solution using either the
northwest corner method or Vogel’s
approximation method
3. Calculate an improvement index for each empty
cell using either the stepping-stone method or
the MODI method. If improvement indices are all
nonnegative, stop as the optimal solution has
been found. If any index is negative, continue to
step 4.
4. Select the cell with the improvement index
indicating the greatest decrease in cost. Fill this
cell using the stepping-stone path and go to step
3.
© 2009 Prentice-Hall, Inc. 10 – 47
Using Excel QM to Solve
Transportation Problems
Excel QM input screen and formulas
Program 10.1A
© 2009 Prentice-Hall, Inc. 10 – 48
Using Excel QM to Solve
Transportation Problems
Output from Excel QM with optimal solution
Program 10.1B
© 2009 Prentice-Hall, Inc. 10 – 49
MODI Method
The MODI (modified distribution)
distribution method allows
us to compute improvement indices quickly for
each unused square without drawing all of the
closed paths
Because of this, it can often provide considerable
time savings over the stepping-stone method for
solving transportation problems
If there is a negative improvement index, then only
one stepping-stone path must be found
This is used in the same manner as before to
obtain an improved solution
© 2009 Prentice-Hall, Inc. 10 – 50
How to Use the MODI Approach
In applying the MODI method, we begin with an
initial solution obtained by using the northwest
corner rule
We now compute a value for each row (call the
values R1, R2, R3 if there are three rows) and for
each column (K1, K2, K3) in the transportation table
In general we let
Ri = value for assigned row i
Kj = value for assigned column j
Cij = cost in square ij (cost of
shipping from source i to destination j)
© 2009 Prentice-Hall, Inc. 10 – 51
Five Steps in the MODI Method to
Test Unused Squares
1. Compute the values for each row and column, set
Ri + Kj = Cij
but only for those squares that are currently used
or occupied
2. After all equations have been written, set R1 = 0
3. Solve the system of equations for R and K values
4. Compute the improvement index for each unused
square by the formula
Improvement Index (Iij) = Cij – Ri – Kj
5. Select the best negative index and proceed to
solve the problem as you did using the stepping-
stone method
© 2009 Prentice-Hall, Inc. 10 – 52
Solving the Executive Furniture
Corporation Problem with MODI
The initial northwest corner solution is repeated
in Table 10.10
Note that to use the MODI method we have added
the Ris (rows) and Kjs (columns)
Kj K1 K2 K3
TO FACTORY
Ri A B C
FROM CAPACITY
$5 $4 $3
R1 D 100 100
$8 $4 $3
R2 E 200 100 300
$9 $7 $5
R3 F 100 200 300
WAREHOUSE
300 200 200 700
REQUIREMENTS
Table 10.10
© 2009 Prentice-Hall, Inc. 10 – 53
Solving the Executive Furniture
Corporation Problem with MODI
The first step is to set up an equation for each
occupied square
By setting R1 = 0 we can easily solve for K1, R2, K2,
R3, and K3
(1) R1 + K1 = 5 0 + K1 = 5 K1 = 5
(2) R2 + K1 = 8 R2 + 5 = 8 R2 = 3
(3) R2 + K2 = 4 3 + K2 = 4 K2 = 1
(4) R3 + K2 = 7 R3 + 1 = 7 R3 = 6
(5) R3 + K3 = 5 6 + K3 = 5 K3 = –1
© 2009 Prentice-Hall, Inc. 10 – 54
Solving the Executive Furniture
Corporation Problem with MODI
The next step is to compute the improvement
index for each unused cell using the formula
Improvement index (Iij) = Cij – Ri – Kj
We have
Des Moines- IDB = C12 – R1 – K2 = 4 – 0 – 1
Boston index = +$3
Des Moines- IDC = C13 – R1 – K3 = 3 – 0 – (–1)
Cleveland index = +$4
Evansville- IEC = C23 – R2 – K3 = 3 – 3 – (–1)
Cleveland index = +$1
Fort Lauderdale- I = C31 – R3 – K1 = 9 – 6 – 5
Albuquerque indexFA
= –$2
© 2009 Prentice-Hall, Inc. 10 – 55
Solving the Executive Furniture
Corporation Problem with MODI
The steps we follow to develop an improved
solution after the improvement indices have been
computed are
1. Beginning at the square with the best
improvement index, trace a closed path back
to the original square via squares that are
currently being used
2. Beginning with a plus sign at the unused
square, place alternate minus signs and plus
signs on each corner square of the closed
path just traced
© 2009 Prentice-Hall, Inc. 10 – 56
Solving the Executive Furniture
Corporation Problem with MODI
3. Select the smallest quantity found in those
squares containing the minus signs and add
that number to all squares on the closed path
with plus signs; subtract the number from
squares with minus signs
4. Compute new improvement indices for this
new solution using the MODI method
Note that new Ri and Kj values must be
calculated
Follow this procedure for the second and third
solutions
© 2009 Prentice-Hall, Inc. 10 – 57
Vogel’s Approximation Method:
Another Way To Find An Initial Solution
Vogel’s Approximation Method (VAM)
VAM is not as
simple as the northwest corner method, but it
provides a very good initial solution, often one
that is the optimal solution
VAM tackles the problem of finding a good initial
solution by taking into account the costs
associated with each route alternative
This is something that the northwest corner rule
does not do
To apply VAM, we first compute for each row and
column the penalty faced if we should ship over
the second-best route instead of the least-cost
route
© 2009 Prentice-Hall, Inc. 10 – 58
Vogel’s Approximation Method
The six steps involved in determining an initial
VAM solution are illustrated below beginning with
the same layout originally shown in Table 10.2
VAM Step 1.1 For each row and column of the
transportation table, find the difference between
the distribution cost on the best route in the row
or column and the second best route in the row or
column
This is the opportunity cost of not using the
best route
Step 1 has been done in Table 10.11
© 2009 Prentice-Hall, Inc. 10 – 59
Vogel’s Approximation Method
Transportation table with VAM row and column
differences shown
OPPORTUNITY
3 0 0 COSTS
TO TOTAL
FROM
A B C AVAILABLE
$5 $4 $3
D 100 100 1
$8 $4 $3
E 200 100 300 1
$9 $7 $5
F 100 200 300 2
TOTAL REQUIRED 300 200 200 700
Table 10.11
© 2009 Prentice-Hall, Inc. 10 – 60
Vogel’s Approximation Method
VAM Step 2.
2 identify the row or column with the
greatest opportunity cost, or difference (column A in
this example)
VAM Step 3.Assign
3 as many units as possible to the
lowest-cost square in the row or column selected
VAM Step 4.
4 Eliminate any row or column that has
been completely satisfied by the assignment just
made by placing Xs in each appropriate square
VAM Step 5.
5 Recompute the cost differences for the
transportation table, omitting rows or columns
eliminated in the previous step
© 2009 Prentice-Hall, Inc. 10 – 61
Vogel’s Approximation Method
VAM assignment with D’s requirements satisfied
OPPORTUNITY
31 03 02 COSTS
TO TOTAL
FROM
A B C AVAILABLE
$5 $4 $3
D 100 X X 100 1
$8 $4 $3
E 300 1
$9 $7 $5
F 300 2
TOTAL REQUIRED 300 200 200 700
Table 10.12
© 2009 Prentice-Hall, Inc. 10 – 62
Vogel’s Approximation Method
VAM Step 6. 6 Return to step 2 for the rows and
columns remaining and repeat the steps until an
initial feasible solution has been obtained
In this case column B now has the greatest
difference, 3
We assign 200 units to the lowest-cost square in
the column, EB
We recompute the differences and find the
greatest difference is now in row E
We assign 100 units to the lowest-cost square in
the column, EC
© 2009 Prentice-Hall, Inc. 10 – 63
Vogel’s Approximation Method
Second VAM assignment with B’s requirements
satisfied
OPPORTUNITY
31 03 02 COSTS
TO TOTAL
FROM
A B C AVAILABLE
$5 $4 $3
D 100 X X 100 1
$8 $4 $3
E 200 300 1
$9 $7 $5
F X 300 2
TOTAL REQUIRED 300 200 200 700
Table 10.13
© 2009 Prentice-Hall, Inc. 10 – 64
Vogel’s Approximation Method
Third VAM assignment with E’s requirements
satisfied
TO TOTAL
FROM
A B C AVAILABLE
$5 $4 $3
D 100 X X 100
$8 $4 $3
E X 200 100 300
$9 $7 $5
F X 300
TOTAL REQUIRED 300 200 200 700
Table 10.14
© 2009 Prentice-Hall, Inc. 10 – 65
Vogel’s Approximation Method
Final assignments to balance column and row
requirements
TO TOTAL
FROM
A B C AVAILABLE
$5 $4 $3
D 100 X X 100
$8 $4 $3
E X 200 100 300
$9 $7 $5
F 200 X 100 300
TOTAL REQUIRED 300 200 200 700
Table 10.15
© 2009 Prentice-Hall, Inc. 10 – 66
Unbalanced Transportation Problems
In real-life problems, total demand is frequently
not equal to total supply
These unbalanced problems can be handled
easily by introducing dummy sources or dummy
destinations
If total supply is greater than total demand, a
dummy destination (warehouse), with demand
exactly equal to the surplus, is created
If total demand is greater than total supply, we
introduce a dummy source (factory) with a supply
equal to the excess of demand over supply
© 2009 Prentice-Hall, Inc. 10 – 67
Unbalanced Transportation Problems
In either case, shipping cost coefficients of zero
are assigned to each dummy location or route as
no goods will actually be shipped
Any units assigned to a dummy destination
represent excess capacity
Any units assigned to a dummy source represent
unmet demand
© 2009 Prentice-Hall, Inc. 10 – 68
Demand Less Than Supply
Suppose that the Des Moines factory increases its
rate of production from 100 to 250 desks
The firm is now able to supply a total of 850 desks
each period
Warehouse requirements remain the same (700) so
the row and column totals do not balance
We add a dummy column that will represent a fake
warehouse requiring 150 desks
This is somewhat analogous to adding a slack
variable
We use the northwest corner rule and either
stepping-stone or MODI to find the optimal
solution
© 2009 Prentice-Hall, Inc. 10 – 69
Demand Less Than Supply
Initial solution to an unbalanced problem where
demand is less than supply
TO DUMMY TOTAL
FROM
A B C WAREHOUSE AVAILABLE
$5 $4 $3 0
D 250 250
$8 $4 $3 0
E 50 200 50 300
$9 $7 $5 0
F 150 150 300
WAREHOUSE 300 200 200 150 850
REQUIREMENTS
Total cost = 250($5) + 50($8) + 200($4) + 50($3) + 150($5) + 150(0) = $3,350
New Des Moines
Table 10.16
capacity
© 2009 Prentice-Hall, Inc. 10 – 70
Demand Greater than Supply
The second type of unbalanced condition occurs
when total demand is greater than total supply
In this case we need to add a dummy row
representing a fake factory
The new factory will have a supply exactly equal
to the difference between total demand and total
real supply
The shipping costs from the dummy factory to
each destination will be zero
© 2009 Prentice-Hall, Inc. 10 – 71
Demand Greater than Supply
Unbalanced transportation table for Happy
Sound Stereo Company
TO WAREHOUSE WAREHOUSE WAREHOUSE PLANT
FROM A B C SUPPLY
$6 $4 $9
PLANT W 200
$10 $5 $8
PLANT X 175
$12 $7 $6
PLANT Y 75
Totals
WAREHOUSE 450
DEMAND
250 100 150
500
do not
balance
Table 10.17
© 2009 Prentice-Hall, Inc. 10 – 72
Demand Greater than Supply
Initial solution to an unbalanced problem in
which demand is greater than supply
TO WAREHOUSE WAREHOUSE WAREHOUSE PLANT
FROM A B C SUPPLY
$6 $4 $9
PLANT W 200 200
$10 $5 $8
PLANT X 50 100 25 175
$12 $7 $6
PLANT Y 75 75
0 0 0
PLANT Y 50 50
WAREHOUSE 250 100 150 500
DEMAND
Total cost of initial solution = 200($6) + 50($10) + 100($5) + 25($8) + 75($6)
+ $50(0) = $2,850
Table 10.18
© 2009 Prentice-Hall, Inc. 10 – 73
Degeneracy in Transportation
Problems
Degeneracy occurs when the number of occupied
squares or routes in a transportation table
solution is less than the number of rows plus the
number of columns minus 1
Such a situation may arise in the initial solution
or in any subsequent solution
Degeneracy requires a special procedure to
correct the problem since there are not enough
occupied squares to trace a closed path for each
unused route and it would be impossible to apply
the stepping-stone method or to calculate the R
and K values needed for the MODI technique
© 2009 Prentice-Hall, Inc. 10 – 74
Degeneracy in Transportation
Problems
To handle degenerate problems, create an
artificially occupied cell
That is, place a zero (representing a fake
shipment) in one of the unused squares and then
treat that square as if it were occupied
The square chosen must be in such a position as
to allow all stepping-stone paths to be closed
There is usually a good deal of flexibility in
selecting the unused square that will receive the
zero
© 2009 Prentice-Hall, Inc. 10 – 75
Degeneracy in an Initial Solution
The Martin Shipping Company example illustrates
degeneracy in an initial solution
They have three warehouses which supply three
major retail customers
Applying the northwest corner rule the initial
solution has only four occupied squares
This is less than the amount required to use
either the stepping-stone or MODI method to
improve the solution (3 rows + 3 columns – 1 = 5)
To correct this problem, place a zero in an
unused square, typically one adjacent to the last
filled cell
© 2009 Prentice-Hall, Inc. 10 – 76
Degeneracy in an Initial Solution
Initial solution of a degenerate problem
TO WAREHOUSE
CUSTOMER 1 CUSTOMER 2 CUSTOMER 3
FROM SUPPLY
$8 $2 $6
WAREHOUSE 1 100 0 100
$10 $9 $9
WAREHOUSE 2 0 100 20 120
$7 $10 $7
WAREHOUSE 3 80 80
CUSTOMER 100 100 100 300
DEMAND
Table 10.19
Possible choices of
cells to address the
degenerate solution
© 2009 Prentice-Hall, Inc. 10 – 77
Degeneracy During
Later Solution Stages
A transportation problem can become degenerate
after the initial solution stage if the filling of an
empty square results in two or more cells
becoming empty simultaneously
This problem can occur when two or more cells
with minus signs tie for the lowest quantity
To correct this problem, place a zero in one of the
previously filled cells so that only one cell
becomes empty
© 2009 Prentice-Hall, Inc. 10 – 78
Degeneracy During
Later Solution Stages
Bagwell Paint Example
After one iteration, the cost analysis at
Bagwell Paint produced a transportation table
that was not degenerate but was not optimal
The improvement indices are
factory A – warehouse 2 index = +2
factory A – warehouse 3 index = +1
factory B – warehouse 3 index = –15
factory C – warehouse 2 index = +11
Only route with
a negative index
© 2009 Prentice-Hall, Inc. 10 – 79
Degeneracy During
Later Solution Stages
Bagwell Paint transportation table
TO WAREHOUSE WAREHOUSE WAREHOUSE
1 2 3 FACTORY
FROM CAPACITY
$8 $5 $16
FACTORY A 70 70
$15 $10 $7
FACTORY B 50 80 130
$3 $9 $10
FACTORY C 30 50 80
WAREHOUSE
150 80 50 280
REQUIREMENT
Table 10.20
© 2009 Prentice-Hall, Inc. 10 – 80
Degeneracy During
Later Solution Stages
Tracing a closed path for the factory B –
warehouse 3 route
TO
WAREHOUSE 1 WAREHOUSE 3
FROM
$15 $7
FACTORY B 50
– +
$3 $10
FACTORY C 30 50
+ –
Table 10.21
This would cause two cells to drop to zero
We need to place an artificial zero in one of these
cells to avoid degeneracy
© 2009 Prentice-Hall, Inc. 10 – 81
More Than One Optimal Solution
It is possible for a transportation problem to have
multiple optimal solutions
This happens when one or more of the
improvement indices zero in the optimal solution
This means that it is possible to design
alternative shipping routes with the same total
shipping cost
The alternate optimal solution can be found by
shipping the most to this unused square using a
stepping-stone path
In the real world, alternate optimal solutions
provide management with greater flexibility in
selecting and using resources
© 2009 Prentice-Hall, Inc. 10 – 82
Maximization Transportation Problems
If the objective in a transportation problem is to
maximize profit, a minor change is required in the
transportation algorithm
Now the optimal solution is reached when all the
improvement indices are negative or zero
The cell with the largest positive improvement
index is selected to be filled using a stepping-
stone path
This new solution is evaluated and the process
continues until there are no positive improvement
indices
© 2009 Prentice-Hall, Inc. 10 – 83
Unacceptable Or Prohibited Routes
At times there are transportation problems in
which one of the sources is unable to ship to one
or more of the destinations
When this occurs, the problem is said to have an
unacceptable or prohibited route
In a minimization problem, such a prohibited
route is assigned a very high cost to prevent this
route from ever being used in the optimal
solution
In a maximization problem, the very high cost
used in minimization problems is given a
negative sign, turning it into a very bad profit
© 2009 Prentice-Hall, Inc. 10 – 84
Facility Location Analysis
The transportation method is especially useful in
helping a firm to decide where to locate a new
factory or warehouse
Each alternative location should be analyzed
within the framework of one overall distribution
system
The new location that yields the minimum cost
for the entire system is the one that should be
chosen
© 2009 Prentice-Hall, Inc. 10 – 85
Locating a New Factory for
Hardgrave Machine Company
Hardgrave Machine produces computer
components at three plants and they ship to four
warehouses
The plants have not been able to keep up with
demand so the firm wants to build a new plant
Two sites are being considered, Seattle and
Birmingham
Data has been collected for each possible
location
Which new location will yield the lowest cost for
the firm in combination with the existing plants
and warehouses
© 2009 Prentice-Hall, Inc. 10 – 86
Locating a New Factory for
Hardgrave Machine Company
Hardgrave’s demand and supply data
MONTHLY
DEMAND PRODUCTION MONTHLY COST TO PRODUCE
WAREHOUSE (UNITS) PLANT SUPPLY ONE UNIT ($)
Detroit 10,000 Cincinnati 15,000 48
Dallas 12,000 Salt Lake 6,000 50
New York 15,000 Pittsburgh 14,000 52
Los Angeles 9,000 35,000
46,000
Supply needed from new plant = 46,000 – 35,000 = 11,000 units per month
Table 10.22
ESTIMATED PRODUCTION COST
PER UNIT AT PROPOSED PLANTS
Seattle $53
Birmingham $49
© 2009 Prentice-Hall, Inc. 10 – 87
Locating a New Factory for
Hardgrave Machine Company
Hardgrave’s shipping costs
TO LOS
FROM DETROIT DALLAS NEW YORK ANGELES
CINCINNATI $25 $55 $40 $60
SALT LAKE 35 30 50 40
PITTSBURGH 36 45 26 66
SEATTLE 60 38 65 27
BIRMINGHAM 35 30 41 50
Table 10.23
© 2009 Prentice-Hall, Inc. 10 – 88
Locating a New Factory for
Hardgrave Machine Company
Optimal solution for the Birmingham location
TO LOS FACTORY
FROM DETROIT DALLAS NEW YORK ANGELES CAPACITY
73 103 88 108
CINCINNATI 10,000 1,000 4,000 15,000
85 80 100 90
SALT LAKE 1,000 5,000 6,000
88 97 78 118
PITTSBURGH 14,000 14,000
84 79 90 99
BIRMINGHAM 11,000 11,000
WAREHOUSE 10,000 12,000 15,000 9,000 46,000
REQUIREMENT
Table 10.24
© 2009 Prentice-Hall, Inc. 10 – 89
Locating a New Factory for
Hardgrave Machine Company
Optimal solution for the Seattle location
TO LOS FACTORY
FROM DETROIT DALLAS NEW YORK ANGELES CAPACITY
73 103 88 108
CINCINNATI 10,000 4,000 1,000 15,000
85 80 100 90
SALT LAKE 6,000 6,000
88 97 78 118
PITTSBURGH 14,000 14,000
113 91 118 80
SEATTLE 2,000 9,000 11,000
WAREHOUSE 10,000 12,000 15,000 9,000 46,000
REQUIREMENT
Table 10.25
© 2009 Prentice-Hall, Inc. 10 – 90
Locating a New Factory for
Hardgrave Machine Company
By comparing the total system costs of the two
alternatives, Hardgrave can select the lowest cost
option
The Birmingham location yields a total system
cost of $3,741,000
The Seattle location yields a total system cost of
$3,704,000
With the lower total system cost, the Seattle
location is favored
Excel QM can also be used as a solution tool
© 2009 Prentice-Hall, Inc. 10 – 91
Locating a New Factory for
Hardgrave Machine Company
Excel input screen
Program 10.2A
© 2009 Prentice-Hall, Inc. 10 – 92
Locating a New Factory for
Hardgrave Machine Company
Output from Excel QM analysis
Program 10.2A
© 2009 Prentice-Hall, Inc. 10 – 93
Assignment Model Approach
The second special-purpose LP algorithm is the
assignment method
Each assignment problem has associated with it
a table, or matrix
Generally, the rows contain the objects or people
we wish to assign, and the columns comprise the
tasks or things we want them assigned to
The numbers in the table are the costs associated
with each particular assignment
An assignment problem can be viewed as a
transportation problem in which the capacity
from each source is 1 and the demand at each
destination is 1
© 2009 Prentice-Hall, Inc. 10 – 94
Assignment Model Approach
The Fix-It Shop has three rush projects to repair
They have three repair persons with different
talents and abilities
The owner has estimates of wage costs for each
worker for each project
The owner’s objective is to assign the three
project to the workers in a way that will result in
the lowest cost to the shop
Each project will be assigned exclusively to one
worker
© 2009 Prentice-Hall, Inc. 10 – 95
Assignment Model Approach
Estimated project repair costs for the Fix-It shop
assignment problem
PROJECT
PERSON 1 2 3
Adams $11 $14 $6
Brown 8 10 11
Cooper 9 12 7
Table 10.26
© 2009 Prentice-Hall, Inc. 10 – 96
Assignment Model Approach
Summary of Fix-It Shop assignment alternatives
and costs
PRODUCT ASSIGNMENT
LABOR TOTAL
1 2 3 COSTS ($) COSTS ($)
Adams Brown Cooper 11 + 10 + 7 28
Adams Cooper Brown 11 + 12 + 11 34
Brown Adams Cooper 8 + 14 + 7 29
Brown Cooper Adams 8 + 12 + 6 26
Cooper Adams Brown 9 + 14 + 11 34
Cooper Brown Adams 9 + 10 + 6 25
Table 10.27
© 2009 Prentice-Hall, Inc. 10 – 97
The Hungarian Method
(Flood’s Technique)
The Hungarian method is an efficient method of
finding the optimal solution to an assignment
problem without having to make direct
comparisons of every option
It operates on the principle of matrix reduction
By subtracting and adding appropriate numbers
in the cost table or matrix, we can reduce the
problem to a matrix of opportunity costs
Opportunity costs show the relative penalty
associated with assigning any person to a project
as opposed to making the best assignment
We want to make assignment so that the
opportunity cost for each assignment is zero
© 2009 Prentice-Hall, Inc. 10 – 98
Three Steps of the Assignment Method
1. Find the opportunity cost table by:
by
(a) Subtracting the smallest number in each row
of the original cost table or matrix from every
number in that row
(b) Then subtracting the smallest number in
each column of the table obtained in part (a)
from every number in that column
2. Test the table resulting from step 1 to see
whether an optimal assignment can be made by
drawing the minimum number of vertical and
horizontal straight lines necessary to cover all
the zeros in the table. If the number of lines is
less than the number of rows or columns,
proceed to step 3.
© 2009 Prentice-Hall, Inc. 10 – 99
Three Steps of the Assignment Method
3. Revise the present opportunity cost table by
subtracting the smallest number not covered by
a line from every other uncovered number. This
same number is also added to any number(s)
lying at the intersection of horizontal and
vertical lines. Return to step 2 and continue the
cycle until an optimal assignment is possible.
© 2009 Prentice-Hall, Inc. 10 – 100
Steps in the Assignment Method
Not
Set up cost table for problem Revise opportunity cost table
optimal
in two steps:
Step 1 (a) Subtract the smallest
number not covered by a line
from itself and every other
Find opportunity cost
uncovered number
(a) Subtract smallest number in
(b) add this number at every
each row from every number
intersection of any two lines
in that row, then
(b) subtract smallest number in
each column from every
number in that column Optimal solution at zero
locations. Systematically make
final assignments.
Step 2
(a) Check each row and
column for a unique zero and
Test opportunity cost table to make the first assignment in
see if optimal assignments are that row or column
possible by drawing the
minimum possible lines on (b) Eliminate that row and
columns and/or rows such that Optimal column and search for another
all zeros are covered unique zero. Make that
assignment and proceed in a
like manner.
Figure 10.3
© 2009 Prentice-Hall, Inc. 10 – 101
The Hungarian Method
(Flood’s Technique)
Step 1: Find the opportunity cost table
We can compute row opportunity costs and
column opportunity costs
What we need is the total opportunity cost
We derive this by taking the row opportunity
costs and subtract the smallest number in that
column from each number in that column
© 2009 Prentice-Hall, Inc. 10 – 102
The Hungarian Method
(Flood’s Technique)
Cost of each person- Row opportunity
project assignment cost table
PROJECT PROJECT
PERSON 1 2 3 PERSON 1 2 3
Adams $11 $14 $6 Adams $5 $8 $0
Brown 8 10 11 Brown 0 2 3
Cooper 9 12 7 Cooper 2 5 0
Table 10.28 Table 10.29
The opportunity cost of assigning Cooper to
project 2 is $12 – $7 = $5
© 2009 Prentice-Hall, Inc. 10 – 103
The Hungarian Method
(Flood’s Technique)
We derive the total opportunity costs by taking
the costs in Table 29 and subtract the smallest
number in each column from each number in that
column
Row opportunity Total opportunity
cost table cost table
PROJECT PROJECT
PERSON 1 2 3 PERSON 1 2 3
Adams $5 $8 $0 Adams $5 $6 $0
Brown 0 2 3 Brown 0 0 3
Cooper 2 5 0 Cooper 2 3 0
Table 10.29 Table 10.30
© 2009 Prentice-Hall, Inc. 10 – 104
The Hungarian Method
(Flood’s Technique)
Step 2: Test for the optimal assignment
We want to assign workers to projects in such
a way that the total labor costs are at a
minimum
We would like to have a total assigned
opportunity cost of zero
The test to determine if we have reached an
optimal solution is simple
We find the minimum number of straight lines
necessary to cover all the zeros in the table
If the number of lines equals the number of
rows or columns, an optimal solution has been
reached
© 2009 Prentice-Hall, Inc. 10 – 105
The Hungarian Method
(Flood’s Technique)
Test for optimal solution
PROJECT
PERSON 1 2 3
Adams $5 $6 $0
Brown 0 0 3 Covering line 1
Cooper 2 3 0
Table 10.31 Covering line 2
This requires only two lines to cover the zeros so
the solution is not optimal
© 2009 Prentice-Hall, Inc. 10 – 106
The Hungarian Method
(Flood’s Technique)
Step 3: Revise the opportunity-cost table
We subtract the smallest number not covered
by a line from all numbers not covered by a
straight line
The same number is added to every number
lying at the intersection of any two lines
We then return to step 2 to test this new table
© 2009 Prentice-Hall, Inc. 10 – 107
The Hungarian Method
(Flood’s Technique)
Revised opportunity cost table (derived by
subtracting 2 from each cell not covered by a line
and adding 2 to the cell at the intersection of the
lines)
PROJECT
PERSON 1 2 3
Adams $3 $4 $0
Brown 0 0 5
Cooper 0 1 0
Table 10.32
© 2009 Prentice-Hall, Inc. 10 – 108
The Hungarian Method
(Flood’s Technique)
Optimality test on the revised opportunity cost
table
PROJECT
PERSON 1 2 3
Adams $3 $4 $0
Brown 0 0 5 Covering line 2
Cooper 0 1 0
Table 10.33 Covering line 1 Covering line 3
This requires three lines to cover the zeros so the
solution is optimal
© 2009 Prentice-Hall, Inc. 10 – 109
Making the Final Assignment
The optimal assignment is Adams to project 3,
Brown to project 2, and Cooper to project 1
But this is a simple problem
For larger problems one approach to making the
final assignment is to select a row or column that
contains only one zero
Make the assignment to that cell and rule out its
row and column
Follow this same approach for all the remaining
cells
© 2009 Prentice-Hall, Inc. 10 – 110
Making the Final Assignment
Total labor costs of this assignment are
ASSIGNMENT COST ($)
Adams to project 3 6
Brown to project 2 10
Cooper to project 1 9
Total cost 25
© 2009 Prentice-Hall, Inc. 10 – 111
Making the Final Assignment
Making the final assignments
(A) FIRST (B) SECOND (C) THIRD
ASSIGNMENT ASSIGNMENT ASSIGNMENT
1 2 3 1 2 3 1 2 3
Adams 3 4 0 Adams 3 4 0 Adams 3 4 0
Brown 0 0 5 Brown 0 0 5 Brown 0 0 5
Cooper 0 1 0 Cooper 0 1 0 Cooper 0 1 0
Table 10.34
© 2009 Prentice-Hall, Inc. 10 – 112
Using Excel QM for the Fix-It Shop
Assignment Problem
Excel QM assignment module
Program 10.3A
© 2009 Prentice-Hall, Inc. 10 – 113
Using Excel QM for the Fix-It Shop
Assignment Problem
Excel QM output screen
Program 10.3A
© 2009 Prentice-Hall, Inc. 10 – 114
Unbalanced Assignment Problems
Often the number of people or objects to be
assigned does not equal the number of tasks or
clients or machines listed in the columns, and the
problem is unbalanced
When this occurs, and there are more rows than
columns, simply add a dummy column or task
If the number of tasks exceeds the number of
people available, we add a dummy row
Since the dummy task or person is nonexistent,
we enter zeros in its row or column as the cost or
time estimate
© 2009 Prentice-Hall, Inc. 10 – 115
Unbalanced Assignment Problems
The Fix-It Shop has another worker available
The shop owner still has the same basic problem
of assigning workers to projects
But the problem now needs a dummy column to
balance the four workers and three projects
PROJECT
PERSON 1 2 3 DUMMY
Adams $11 $14 $6 $0
Brown 8 10 11 0
Cooper 9 12 7 0
Davis 10 13 8 0
Table 10.35
© 2009 Prentice-Hall, Inc. 10 – 116
Maximization Assignment Problems
Some assignment problems are phrased in terms
of maximizing the payoff, profit, or effectiveness
It is easy to obtain an equivalent minimization
problem by converting all numbers in the table to
opportunity costs
This is brought about by subtracting every
number in the original payoff table from the
largest single number in that table
Transformed entries represent opportunity costs
Once the optimal assignment has been found, the
total payoff is found by adding the original
payoffs of those cells that are in the optimal
assignment
© 2009 Prentice-Hall, Inc. 10 – 117
Maximization Assignment Problems
The British navy wishes to assign four ships to
patrol four sectors of the North Sea
Ships are rated for their probable efficiency in
each sector
The commander wants to determine patrol
assignments producing the greatest overall
efficiencies
© 2009 Prentice-Hall, Inc. 10 – 118
Maximization Assignment Problems
Efficiencies of British ships in patrol sectors
SECTOR
SHIP A B C D
1 20 60 50 55
2 60 30 80 75
3 80 100 90 80
4 65 80 75 70
Table 10.36
© 2009 Prentice-Hall, Inc. 10 – 119
Maximization Assignment Problems
Opportunity cost of British ships
SECTOR
SHIP A B C D
1 80 40 50 45
2 40 70 20 25
3 20 0 10 20
4 35 20 25 30
Table 10.37
© 2009 Prentice-Hall, Inc. 10 – 120
Maximization Assignment Problems
First convert the maximization efficiency table
into a minimizing opportunity cost table by
subtracting each rating from 100, the largest
rating in the whole table
The smallest number in each row is subtracted
from every number in that row and the smallest
number in each column is subtracted from every
number in that column
The minimum number of lines needed to cover
the zeros in the table is four, so this represents
an optimal solution
© 2009 Prentice-Hall, Inc. 10 – 121
Maximization Assignment Problems
The overall efficiency
ASSIGNMENT EFFICIENCY
Ship 1 to sector D 55
Ship 2 to sector C 80
Ship 3 to sector B 100
Ship 4 to sector A 65
Total efficiency 300
© 2009 Prentice-Hall, Inc. 10 – 122
Chapter 3
Decision Analysis
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Chapter Outline
3.1 Introduction
3.2 The Six Steps in Decision Making
3.3 Types of Decision-Making
Environments
3.4 Decision Making under Uncertainty
3.5 Decision Making under Risk
3.5.1 EMV
3.5.2 Sensitivity Analysis
© 2009 Prentice-Hall, Inc. 3–2
Introduction
What is involved in making a good
decision?
Decision theory is an analytic and
systematic approach to the study of
decision making
A good decision is one that is based
on logic, considers all available data
and possible alternatives, and the
quantitative approach described here
© 2009 Prentice-Hall, Inc. 3–3
The Six Steps in Decision Making
1. Clearly define the problem at hand
2. List the possible alternatives
3. Identify the possible outcomes or states
of nature
4. List the payoff or profit of each
combination of alternatives and
outcomes
5. Select one of the mathematical decision
theory models
6. Apply the model and make your decision
© 2009 Prentice-Hall, Inc. 3–4
Thompson Lumber Company
Step 1 – Define the problem
Expand by manufacturing and
marketing a new product, backyard
storage sheds
Step 2 – List alternatives
Construct a large new plant
A small plant
No plant at all
Step 3 – Identify possible outcomes
The market could be favorable or
unfavorable
© 2009 Prentice-Hall, Inc. 3–5
Thompson Lumber Company
Step 4 – List the payoffs
Identify conditional values for the
profits for large, small, and no plants
for the two possible market
conditions
STATE OF NATURE
FAVORABLE UNFAVORABLE
ALTERNATIVE MARKET ($) MARKET ($)
Construct a large plant 200,000 –180,000
Construct a small plant 100,000 –20,000
Do nothing 0 0
© 2009 Prentice-Hall, Inc. 3–6
Thompson Lumber Company
Step 5 – Select the decision model
Depends on the environment and amount
of risk and uncertainty
Step 6 – Apply the model to the data
Solution and analysis used to help the
decision making
© 2009 Prentice-Hall, Inc. 3–7
Types of Decision-Making
Environments
Type 1: Decision making under certainty
Decision maker knows with certainty the
consequences of every alternative or
decision choice
Type 2: Decision making under uncertainty
The decision maker does not know the
probabilities of the various outcomes
Type 3: Decision making under risk
The decision maker knows the
probabilities of the various outcomes
© 2009 Prentice-Hall, Inc. 3–8
Decision Making Under
Uncertainty
There are several criteria for making decisions
under uncertainty
1. Maximax (optimistic)
2. Maximin (pessimistic)
3. Criterion of realism (Hurwicz)
4. Equally likely (Laplace)
5. Minimax regret
© 2009 Prentice-Hall, Inc. 3–9
Maximax
Used to find the alternative that maximizes
the maximum payoff
Locate the maximum payoff for each alternative
Select the alternative with the maximum
number
STATE OF NATURE
FAVORABLE UNFAVORABLE MAXIMUM IN
ALTERNATIVE MARKET ($) MARKET ($) A ROW ($)
Construct a large 200,000 –180,000 200,000
plant
Maximax
Construct a small 100,000 –20,000 100,000
plant
Do nothing 0 0 0
Table 3.2
© 2009 Prentice-Hall, Inc. 3 – 10
Maximin
Used to find the alternative that maximizes
the minimum payoff
Locate the minimum payoff for each alternative
Select the alternative with the maximum
number
STATE OF NATURE
FAVORABLE UNFAVORABLE MINIMUM IN
ALTERNATIVE MARKET ($) MARKET ($) A ROW ($)
Construct a large 200,000 –180,000 –180,000
plant
Construct a small 100,000 –20,000 –20,000
plant
Do nothing 0 0 0
Table 3.3
Maximin
© 2009 Prentice-Hall, Inc. 3 – 11
Criterion of Realism (Hurwicz)
A weighted average compromise between
optimistic and pessimistic
Select a coefficient of realism α
Coefficient is between 0 and 1
A value of 1 is 100% optimistic
Compute the weighted averages for each
alternative
Select the alternative with the highest value
Weighted average = α (maximum in row)
+ (1 – α )(minimum in row)
© 2009 Prentice-Hall, Inc. 3 – 12
Criterion of Realism (Hurwicz)
For the large plant alternative using α = 0.8
(0.8)(200,000) + (1 – 0.8)(–180,000) = 124,000
For the small plant alternative using α = 0.8
(0.8)(100,000) + (1 – 0.8)(–20,000) = 76,000
STATE OF NATURE
CRITERION
FAVORABLE UNFAVORABLE OF REALISM
ALTERNATIVE MARKET ($) MARKET ($) (α = 0.8)$
Construct a large
200,000 –180,000 124,000
plant
Realism
Construct a small
100,000 –20,000 76,000
plant
Do nothing 0 0 0
Table 3.4
© 2009 Prentice-Hall, Inc. 3 – 13
Equally Likely (Laplace)
Considers all the payoffs for each alternative
Find the average payoff for each alternative
Select the alternative with the highest average
STATE OF NATURE
FAVORABLE UNFAVORABLE ROW
ALTERNATIVE MARKET ($) MARKET ($) AVERAGE ($)
Construct a large 200,000 –180,000 10,000
plant
Construct a small
100,000 –20,000 40,000
plant
Equally likely
Do nothing 0 0 0
Table 3.5
© 2009 Prentice-Hall, Inc. 3 – 14
Minimax Regret
Based on opportunity loss or regret,
regret the
difference between the optimal profit and
actual payoff for a decision.
Create an opportunity loss table by determining
the opportunity loss for not choosing the best
alternative
© 2009 Prentice-Hall, Inc. 3 – 15
Minimax Regret
Opportunity loss is calculated by subtracting
each payoff in the column from the best payoff in
the column STATE OF NATURE
FAVORABLE UNFAVORABLE
MARKET ($) MARKET ($)
Table 3.6 200,000 – 200,000 0 – (–180,000)
200,000 – 100,000 0 – (–20,000)
200,000 – 0 0–0
STATE OF NATURE
FAVORABLE UNFAVORABLE
ALTERNATIVE MARKET ($) MARKET ($)
Construct a large plant 0 180,000 Table 3.7
Construct a small plant 100,000 20,000
Do nothing 200,000 0© 2009 Prentice-Hall, Inc. 3 – 16
Minimax Regret
Find the maximum opportunity loss for each
alternative and pick the alternative with the
minimum number
STATE OF NATURE
FAVORABLE UNFAVORABLE MAXIMUM IN
ALTERNATIVE MARKET ($) MARKET ($) A ROW ($)
Construct a large
0 180,000 180,000
plant
Construct a small 100,000 20,000 100,000
plant
Minimax
Do nothing 200,000 0 200,000
Table 3.8
© 2009 Prentice-Hall, Inc. 3 – 17
Decision Making Under Risk
Decision making when there are several possible
states of nature and we know the probabilities
associated with each possible state
Most popular method is to choose the alternative
with the highest expected monetary value (EMV)
native i) = (payoff of first state of nature)
x (probability of first state of nature)
+ (payoff of second state of nature)
x (probability of second state of nature)
+ … + (payoff of last state of nature)
x (probability of last state of nature)
© 2009 Prentice-Hall, Inc. 3 – 18
EMV for Thompson Lumber
Each market has a probability of 0.50
Which alternative would give the highest EMV?
The calculations are
rge plant) = (0.50)($200,000) + (0.50)(–$180,000)
= $10,000
mall plant) = (0.50)($100,000) + (0.50)(–$20,000)
= $40,000
o nothing) = (0.50)($0) + (0.50)($0)
= $0
© 2009 Prentice-Hall, Inc. 3 – 19
EMV for Thompson Lumber
STATE OF NATURE
FAVORABLE UNFAVORABLE
ALTERNATIVE MARKET ($) MARKET ($) EMV ($)
Construct a large 200,000 –180,000 10,000
plant
Construct a small 100,000 –20,000 40,000
plant
Do nothing 0 0 0
Probabilities 0.50 0.50
Table 3.9 Largest EMV
© 2009 Prentice-Hall, Inc. 3 – 20
Expected Value of Perfect
Information (EVPI)
EVPI places an upper bound on what you should
pay for additional information
EVPI = EVwPI – Maximum EMV
EVwPI is the long run average return if we have
perfect information before a decision is made
EVwPI = (best payoff for first state of nature)
x (probability of first state of nature)
+ (best payoff for second state of nature)
x (probability of second state of nature)
+ … + (best payoff for last state of nature)
x (probability of last state of nature)
© 2009 Prentice-Hall, Inc. 3 – 21
Expected Value of Perfect
Information (EVPI)
Scientific Marketing, Inc. offers analysis
that will provide certainty about market
conditions (favorable)
Additional information will cost $65,000
Is it worth purchasing the information?
© 2009 Prentice-Hall, Inc. 3 – 22
Expected Value of Perfect
Information (EVPI)
1. Best alternative for favorable state of nature is
build a large plant with a payoff of $200,000
Best alternative for unfavorable state of nature is
to do nothing with a payoff of $0
EVwPI = ($200,000)(0.50) + ($0)(0.50) = $100,000
2. The maximum EMV without additional
information is $40,000
EVPI = EVwPI – Maximum EMV
= $100,000 - $40,000
= $60,000
© 2009 Prentice-Hall, Inc. 3 – 23
Expected Value of Perfect
Information (EVPI)
1. Best alternative for favorable state of nature is
build a large plant with a payoff of $200,000
So the maximum Thompson
Best alternative for unfavorable
should pay for the state of nature is
additional
to do nothinginformation
with a payoff
is of $0
$60,000
EVwPI = ($200,000)(0.50) + ($0)(0.50) = $100,000
2. The maximum EMV without additional
information is $40,000
EVPI = EVwPI – Maximum EMV
= $100,000 - $40,000
= $60,000
© 2009 Prentice-Hall, Inc. 3 – 24
Expected Opportunity Loss
Expected opportunity loss (EOL) is the
cost of not picking the best solution
First construct an opportunity loss table
For each alternative, multiply the
opportunity loss by the probability of that
loss for each possible outcome and add
these together
Minimum EOL will always result in the
same decision as maximum EMV
Minimum EOL will always equal EVPI
© 2009 Prentice-Hall, Inc. 3 – 25
Expected Opportunity Loss
STATE OF NATURE
FAVORABLE UNFAVORABLE
ALTERNATIVE MARKET ($) MARKET ($) EOL
Construct a large plant 0 180,000 90,000
Construct a small 100,000 20,000 60,000
plant
Do nothing 200,000 0 100,000
Probabilities 0.50 0.50
Table 3.10 Minimum EOL
arge plant) = (0.50)($0) + (0.50)($180,000)
= $90,000
mall plant) = (0.50)($100,000) + (0.50)($20,000)
= $60,000
o nothing) = (0.50)($200,000) + (0.50)($0)
= $100,000
© 2009 Prentice-Hall, Inc. 3 – 26
Sensitivity Analysis
Sensitivity analysis examines how our decision
might change with different input data
For the Thompson Lumber example
P = probability of a favorable market
(1 – P) = probability of an unfavorable market
© 2009 Prentice-Hall, Inc. 3 – 27
Sensitivity Analysis
EMV(Large Plant) = $200,000P – $180,000)(1 – P)
= $200,000P – $180,000 + $180,000P
= $380,000P – $180,000
EMV(Small Plant) = $100,000P – $20,000)(1 – P)
= $100,000P – $20,000 + $20,000P
= $120,000P – $20,000
EMV(Do Nothing) = $0P + 0(1 – P)
= $0
© 2009 Prentice-Hall, Inc. 3 – 28
Sensitivity Analysis
EMV Values
$300,000
$200,000
Point 2 EMV (large plant)
$100,000
Point 1 EMV (small plant)
0
EMV (do nothing)
–$100,000
.167 .615 1
–$200,000 Values of P
Figure 3.1
© 2009 Prentice-Hall, Inc. 3 – 29
Sensitivity Analysis
Point 1:
EMV(do nothing) = EMV(small plant)
20,000
0 = $120,000 P − $20,000 P= = 0.167
120,000
Point 2:
EMV(small plant) = EMV(large plant)
$120,000 P − $20,000 = $380,000 P − $180,000
160,000
P= = 0.615
260,000
© 2009 Prentice-Hall, Inc. 3 – 30
Sensitivity Analysis
BEST RANGE OF P
ALTERNATIVE VALUES
Do nothing Less than 0.167
EMV Values Construct a small plant 0.167 – 0.615
$300,000 Construct a large plant Greater than 0.615
$200,000
Point 2 EMV (large plant)
$100,000
Point 1 EMV (small plant)
0
EMV (do nothing)
–$100,000
.167 .615 1
–$200,000 Values of P
Figure 3.1
© 2009 Prentice-Hall, Inc. 3 – 31
Using Excel QM to Solve
Decision Theory Problems
Program 3.1A
© 2009 Prentice-Hall, Inc. 3 – 32
Using Excel QM to Solve
Decision Theory Problems
Program 3.1B
© 2009 Prentice-Hall, Inc. 3 – 33
Decision Trees
Any problem that can be presented in a
decision table can also be graphically
represented in a decision tree
Decision trees are most beneficial when a
sequence of decisions must be made
All decision trees contain decision points
or nodes and state-of-nature points or
nodes
A decision node from which one of several
alternatives may be chosen
A state-of-nature node out of which one state
of nature will occur
© 2009 Prentice-Hall, Inc. 3 – 34
Five Steps to
Decision Tree Analysis
1. Define the problem
2. Structure or draw the decision tree
3. Assign probabilities to the states of
nature
4. Estimate payoffs for each possible
combination of alternatives and states of
nature
5. Solve the problem by computing
expected monetary values (EMVs) for
each state of nature node
© 2009 Prentice-Hall, Inc. 3 – 35
Structure of Decision Trees
Trees start from left to right
Represent decisions and outcomes in
sequential order
Squares represent decision nodes
Circles represent states of nature nodes
Lines or branches connect the decisions
nodes and the states of nature
© 2009 Prentice-Hall, Inc. 3 – 36
Thompson’s Decision Tree
A State-of-Nature Node
Favorable Market
A Decision Node
1
Unfavorable Market
uct nt
r
n st Pla
e
Co arg
L Favorable Market
Construct
2
Small Plant Unfavorable Market
Do
No
th
in
g
Figure 3.2
© 2009 Prentice-Hall, Inc. 3 – 37
Thompson’s Decision Tree
EMV for Node = (0.5)($200,000) + (0.5)(–$180,000)
1 = $10,000
Payoffs
Favorable Market (0.5)
$200,000
Alternative with best
EMV is selected 1
Unfavorable Market (0.5)
ct nt –$180,000
r u
n st Pla
e
Co arg
L Favorable Market (0.5)
$100,000
Construct
2
Small Plant Unfavorable Market (0.5)
–$20,000
Do
No
th EMV for Node = (0.5)($100,000)
in
g 2 = $40,000 + (0.5)(–$20,000)
Figure 3.3
$0
© 2009 Prentice-Hall, Inc. 3 – 38
Thompson’s Complex Decision Tree
First Decision Second Decision Payoffs
Point Point
Favorable Market (0.78)
$190,000
nt 2 Unfavorable Market (0.22)
P la –$190,000
ge
Lar Favorable Market (0.78)
$90,000
) Small
. 45 Plant
3 Unfavorable Market (0.22)
(0 –$30,000
ey ts e
urv sul abl No Plant
–$10,000
S
R e vor
a
1 Surv F Favorable Market (0.27)
e $190,000
y
Re y (
ve
nt 4 Unfavorable Market (0.73)
Ne su 0.5 P la
ur
–$190,000
5) ge
tS
ga lts
tiv Lar Favorable Market (0.27)
$90,000
Small
ke
e 5 Unfavorable Market (0.73)
ar
Plant –$30,000
M
t
uc
No Plant
–$10,000
nd
Co
Do Favorable Market (0.50)
Not $200,000
Con t 6 Unfavorable Market (0.50)
duc la n
P –$180,000
tS ge
urv Lar Favorable Market (0.50)
$100,000
e y Small
Plant
7 Unfavorable Market (0.50)
–$20,000
No Plant
$0
Figure 3.4
© 2009 Prentice-Hall, Inc. 3 – 39
Thompson’s Complex Decision Tree
Given favorable survey results,
EMV(node 2) = EMV(large plant | positive survey)
= (0.78)($190,000) + (0.22)(–$190,000) = $106,400
EMV(node 3) = EMV(small plant | positive survey)
= (0.78)($90,000) + (0.22)(–$30,000) = $63,600
EMV for no plant = –$10,000
Given negative survey results,
EMV(node 4) = EMV(large plant | negative survey)
= (0.27)($190,000) + (0.73)(–$190,000) = –$87,400
EMV(node 5) = EMV(small plant | negative survey)
= (0.27)($90,000) + (0.73)(–$30,000) = $2,400
EMV for no plant = –$10,000
© 2009 Prentice-Hall, Inc. 3 – 40
Thompson’s Complex Decision Tree
Compute the expected value of the market survey,
EMV(node 1) = EMV(conduct survey)
= (0.45)($106,400) + (0.55)($2,400)
= $47,880 + $1,320 = $49,200
f the market survey is not conducted,
EMV(node 6) = EMV(large plant)
= (0.50)($200,000) + (0.50)(–$180,000) = $10,000
EMV(node 7) = EMV(small plant)
= (0.50)($100,000) + (0.50)(–$20,000) = $40,000
EMV for no plant = $0
Best choice is to seek marketing information
© 2009 Prentice-Hall, Inc. 3 – 41
Thompson’s Complex Decision Tree
First Decision Second Decision Payoffs
Point Point
$106,400 Favorable Market (0.78)
$190,000
l a nt Unfavorable Market (0.22)
$106,400
P –$190,000
a rge $63,600 Favorable Market (0.78)
L $90,000
) Small
. 45 Plant
Unfavorable Market (0.22)
–$30,000
(0
ey ts e
urv sul abl No Plant
–$10,000
S
R e vor
Su Fa –$87,400 Favorable Market (0.27)
rv $190,000
e
y
Re y (
ve
Ne su 0.5 ant
Pl
Unfavorable Market (0.73)
ur
–$190,000
5) ge
tS
ga lts r $2,400
$2,400
La Favorable Market (0.27)
tiv Small $90,000
ke
e Unfavorable Market (0.73)
ar
Plant –$30,000
Mt
uc
No Plant
–$10,000
d
on
$49,200
C
Do $10,000 Favorable Market (0.50)
Not $200,000
Con t Unfavorable Market (0.50)
duc P la n –$180,000
ge
$40,000
tS urv Lar $40,000 Favorable Market (0.50)
e y Small $100,000
Unfavorable Market (0.50)
Plant –$20,000
No Plant
$0
Figure 3.4
© 2009 Prentice-Hall, Inc. 3 – 42
Expected Value of Sample Information
Thompson wants to know the actual value
of doing the survey
Expected value Expected value
with sample of best decision
EVSI = information, assuming – without sample
no cost to gather it information
= (EV with sample information + cost)
– (EV without sample information)
EVSI = ($49,200 + $10,000) – $40,000 = $19,200
© 2009 Prentice-Hall, Inc. 3 – 43
Sensitivity Analysis
How sensitive are the decisions to
changes in the probabilities?
How sensitive is our decision to the
probability of a favorable survey result?
That is, if the probability of a favorable
result (p = .45) where to change, would we
make the same decision?
How much could it change before we would
make a different decision?
© 2009 Prentice-Hall, Inc. 3 – 44
Sensitivity Analysis
p = probability of a favorable survey
result
(1 – p) = probability of a negative survey
result 1) = ($106,400)p +($2,400)(1 – p)
EMV(node
= $104,000p + $2,400
We are indifferent when the EMV of node 1 is the
same as the EMV of not conducting the survey,
$40,000
$104,000p + $2,400= $40,000
$104,000p = $37,600
p = $37,600/$104,000 = 0.36
© 2009 Prentice-Hall, Inc. 3 – 45
Bayesian Analysis
Many ways of getting probability data
It can be based on
Management’s experience and intuition
Historical data
Computed from other data using Bayes’
theorem
Bayes’ theorem incorporates initial
estimates and information about the
accuracy of the sources
Allows the revision of initial estimates
based on new information
© 2009 Prentice-Hall, Inc. 3 – 46
Calculating Revised Probabilities
In the Thompson Lumber case we used these four
conditional probabilities
P (favorable market(FM) | survey results positive) = 0.78
P (unfavorable market(UM) | survey results positive) = 0.22
P (favorable market(FM) | survey results negative) = 0.27
P (unfavorable market(UM) | survey results negative) = 0.73
The prior probabilities of these markets are
P (FM) = 0.50
P (UM) = 0.50
© 2009 Prentice-Hall, Inc. 3 – 47
Calculating Revised Probabilities
Through discussions with experts Thompson has
learned the following
He can use this information and Bayes’ theorem
to calculate posterior probabilities
STATE OF NATURE
RESULT OF FAVORABLE MARKET UNFAVORABLE MARKET
SURVEY (FM) (UM)
Positive (predicts P (survey positive | FM) P (survey positive | UM)
favorable market
for product) = 0.70 = 0.20
Negative (predicts
unfavorable P (survey negative | FM) P (survey negative | UM)
market for = 0.30 = 0.80
product)
Table 3.11
© 2009 Prentice-Hall, Inc. 3 – 48
Calculating Revised Probabilities
Recall Bayes’ theorem is
P ( B | A ) × P ( A)
P( A | B) =
P ( B | A) × P ( A) + P ( B | A′ ) × P ( A′ )
where
A, B = any two events
A′ = complement of A
For this example, A will represent a favorable
market and B will represent a positive survey
© 2009 Prentice-Hall, Inc. 3 – 49
Calculating Revised Probabilities
P (FM | survey positive)
P ( survey positive | FM ) × P ( FM )
=
P(survey positive |FM) × P(FM) + P(survey positive |UM) × P(UM)
(0.70 )(0.50 ) 0.35
= = = 0.78
(0.70)(0.50) + (0.20)(0.50 ) 0.45
P (UM | survey positive)
P ( survey positive | UM ) × P (UM )
=
P(survey positive |UM) × P(UM) + P(survey positive |FM) × P(FM)
(0.20 )(0.50 ) 0.10
= = = 0.22
(0.20)(0.50) + (0.70)(0.50 ) 0.45
© 2009 Prentice-Hall, Inc. 3 – 50
Calculating Revised Probabilities
POSTERIOR PROBABILITY
CONDITIONAL
PROBABILITY P(STATE OF
P(SURVEY NATURE |
STATE OF POSITIVE | STATE PRIOR JOINT SURVEY
NATURE OF NATURE) PROBABILITY PROBABILITY POSITIVE)
FM 0.70 X 0.50 = 0.35 0.35/0.45 = 0.78
UM 0.20 X 0.50 = 0.10 0.10/0.45 = 0.22
P(survey results positive) = 0.45 1.00
Table 3.12
© 2009 Prentice-Hall, Inc. 3 – 51
Calculating Revised Probabilities
P (FM | survey negative)
P ( survey negative | FM ) × P ( FM )
=
P(survey negative |FM) × P(FM) + P(survey negative |UM) × P(UM)
(0.30 )(0.50 ) 0.15
= = = 0.27
(0.30)(0.50) + (0.80)(0.50 ) 0.55
P (UM | survey negative)
P ( survey negative | UM ) × P (UM )
=
P(survey negative |UM) × P(UM) + P(survey negative |FM) × P(FM)
(0.80 )(0.50 ) 0.40
= = = 0.73
(0.80)(0.50) + (0.30)(0.50 ) 0.55
© 2009 Prentice-Hall, Inc. 3 – 52
Calculating Revised Probabilities
POSTERIOR PROBABILITY
CONDITIONAL
PROBABILITY P(STATE OF
P(SURVEY NATURE |
STATE OF NEGATIVE | STATE PRIOR JOINT SURVEY
NATURE OF NATURE) PROBABILITY PROBABILITY NEGATIVE)
FM 0.30 X 0.50 = 0.15 0.15/0.55 = 0.27
UM 0.80 X 0.50 = 0.40 0.40/0.55 = 0.73
P(survey results positive) = 0.55 1.00
Table 3.13
© 2009 Prentice-Hall, Inc. 3 – 53
Potential Problems Using
Survey Results
We can not always get the necessary
data for analysis
Survey results may be based on cases
where an action was taken
Conditional probability information
may not be as accurate as we would
like
© 2009 Prentice-Hall, Inc. 3 – 54
Utility Theory
Monetary value is not always a true
indicator of the overall value of the
result of a decision
The overall value of a decision is called
utility
Rational people make decisions to
maximize their utility
© 2009 Prentice-Hall, Inc. 3 – 55
Utility Theory
$2,000,000
Accept
Offer
$0
Heads
(0.5)
Reject
Offer
Tails
(0.5)
EMV = $2,500,000 $5,000,000
Figure 3.6
© 2009 Prentice-Hall, Inc. 3 – 56
Utility Theory
Utility assessment assigns the worst outcome
a utility of 0, and the best outcome, a utility of 1
A standard gamble is used to determine utility
values
When you are indifferent, the utility values are
equal
Expected utility of alternative 2 =
Expected utility of alternative 1
Utility of other outcome = (p)
(utility of best outcome, which is 1)
+ (1 – p)(utility of the worst
outcome, which is 0)
Utility of other outcome = (p)
(1) + (1 – p)(0) = p © 2009 Prentice-Hall, Inc. 3 – 57
Standard Gamble
(p)
Best Outcome
Utility = 1
1 (1 – p) Worst Outcome
ve
ati Utility = 0
e rn
Alt
Alt
ern
a ti
ve
2
Other Outcome
Utility = ?
Figure 3.7
© 2009 Prentice-Hall, Inc. 3 – 58
Investment Example
Jane Dickson wants to construct a utility curve
revealing her preference for money between $0
and $10,000
A utility curve plots the utility value versus the
monetary value
An investment in a bank will result in $5,000
An investment in real estate will result in $0 or
$10,000
Unless there is an 80% chance of getting $10,000
from the real estate deal, Jane would prefer to
have her money in the bank
So if p = 0.80, Jane is indifferent between the bank
or the real estate investment
© 2009 Prentice-Hall, Inc. 3 – 59
Investment Example
p = 0.80 $10,000
U($10,000) = 1.0
(1 – p) = 0.20 $0
t in U($0.00) = 0.0
ves tate
In s
alE
Re
Inv
es
t in
Ba
n k
$5,000
Figure 3.8
U($5,000) = p = 1.0
Utility for $5,000 = U($5,000) = pU($10,000) + (1 – p)U($0)
= (0.8)(1) + (0.2)(0) = 0.8
© 2009 Prentice-Hall, Inc. 3 – 60
Investment Example
We can assess other utility values in the same way
For Jane these are
Utility for $7,000 = 0.90
Utility for $3,000 = 0.50
Using the three utilities for different dollar amounts,
she can construct a utility curve
© 2009 Prentice-Hall, Inc. 3 – 61
Utility Curve
1.0 –
U ($10,000) = 1.0
0.9 – U ($7,000) = 0.90
0.8 – U ($5,000) = 0.80
0.7 –
0.6 –
0.5 –
Utility
0.4 –
U ($3,000) = 0.50
0.3 –
0.2 –
0.1 –
U ($0) = 0
| | | | | | | | | | |
$0 $1,000 $3,000 $5,000 $7,000 $10,000
Monetary Value
Figure 3.9
© 2009 Prentice-Hall, Inc. 3 – 62
Utility Curve
Jane’s utility curve is typical of a risk
avoider
A risk avoider gets less utility from greater risk
Avoids situations where high losses might
occur
As monetary value increases, the utility curve
increases at a slower rate
A risk seeker gets more utility from greater risk
As monetary value increases, the utility curve
increases at a faster rate
Someone who is indifferent will have a linear
utility curve
© 2009 Prentice-Hall, Inc. 3 – 63
Utility Curve
Risk
Avoider
e
nc
re
Utility
ffe
di
In
k
is
R
Risk
Seeker
Monetary Outcome
Figure 3.10
© 2009 Prentice-Hall, Inc. 3 – 64
Utility as a
Decision-Making Criteria
Once a utility curve has been developed
it can be used in making decisions
Replace monetary outcomes with utility
values
The expected utility is computed instead
of the EMV
© 2009 Prentice-Hall, Inc. 3 – 65
Utility as a
Decision-Making Criteria
Mark Simkin loves to gamble
He plays a game tossing thumbtacks in
the air
If the thumbtack lands point up, Mark wins
$10,000
If the thumbtack lands point down, Mark
loses $10,000
Should Mark play the game (alternative 1)?
© 2009 Prentice-Hall, Inc. 3 – 66
Utility as a
Decision-Making Criteria
Tack Lands
Point Up (0.45)
$10,000
Tack Lands
1 ame Point Down (0.55)
t ive the G –$10,000
er na ys
a
Alt rk Pl
Ma
Alt
er n
ati
ve
2
Mark Does Not Play the Game
$0
Figure 3.11
© 2009 Prentice-Hall, Inc. 3 – 67
Utility as a
Decision-Making Criteria
Step 1– Define Mark’s utilities
U (–$10,000) = 0.05
U ($0) = 0.15
U ($10,000) = 0.30
Step 2 – Replace monetary values with
utility values
E(alternative 1: play the game) = (0.45)(0.30) + (0.55)(0.05)
= 0.135 + 0.027 = 0.162
E(alternative 2: don’t play the game) = 0.15
© 2009 Prentice-Hall, Inc. 3 – 68
Utility as a
Decision-Making Criteria
1.00 –
0.75 –
Utility
0.50 –
0.30 –
0.25 –
0.15 –
0.05 –
0 |– | | | |
–$20,000 –$10,000 $0 $10,000 $20,000
Monetary Outcome Figure 3.12
© 2009 Prentice-Hall, Inc. 3 – 69
Utility as a
Decision-Making Criteria
Tack Lands Utility
E = 0.162 Point Up (0.45)
0.30
Tack Lands
1 ame Point Down (0.55)
t ive the G 0.05
er na ys
a
Alt rk Pl
Ma
Alt
er n
ati
ve
2
Don’t Play
0.15
Figure 3.13
© 2009 Prentice-Hall, Inc. 3 – 70
Queuing Theory Models
or
Waiting Lines
Learning Objectives
After completing this chapter, we will be able to:
1. Describe the trade-off curves for cost-of-
waiting time and cost of service.
2. Understand the three parts of a queuing
system: the calling population, the queue
itself, and the service facility.
3. Describe the basic queuing system
configurations.
4. Understand the assumptions of the common
models dealt with in this chapter.
5. Analyze a variety of operating characteristics
of waiting lines.
Copyright © 2012 Pearson Education 13-2
Chapter Outline
1 Introduction
2 Waiting Line Costs
3 Characteristics of a Queuing System
4 Single-Channel Queuing Model with
Poisson Arrivals and Exponential Service
Times (M/M/1)
Copyright © 2012 Pearson Education 13-3
Introduction
Queuing theory is the study of waiting lines.
It is one of the oldest and most widely used
quantitative analysis techniques.
The three basic components of a queuing
process are arrivals, service facilities, and the
actual waiting line.
Analytical models of waiting lines can help
managers evaluate the cost and effectiveness
of service systems.
Copyright © 2012 Pearson Education 13-4
Waiting Line Costs
Most waiting line problems are focused on
finding the ideal level of service a firm should
provide.
In most cases, this service level is something
management can control.
When an organization does have control, they
often try to find the balance between two
extremes.
Copyright © 2012 Pearson Education 13-5
Waiting Line Costs
There is generally a trade-off between cost of
providing service and cost of waiting time.
A large staff and many service facilities generally results
in high levels of service but have high costs.
Having the minimum number of service facilities keeps
service cost down but may result in dissatisfied
customers.
Service facilities are evaluated on their total
expected cost which is the sum of service costs
and waiting costs.
Organizations typically want to find the service
level that minimizes the total expected cost.
Copyright © 2012 Pearson Education 13-6
Queuing Costs and Service Levels
Figure 1
Copyright © 2012 Pearson Education 13-7
Characteristics of a Queuing System
There are three parts to a queuing system:
1. The arrivals or inputs to the system
(sometimes referred to as the calling
population).
2. The queue or waiting line itself.
3. The service facility.
These components have their own
characteristics that must be examined
before mathematical models can be
developed.
Copyright © 2012 Pearson Education 13-8
Characteristics of a Queuing System
Arrival Characteristics have three major
characteristics: size, pattern, and behavior.
The size of the calling population can be either
unlimited (essentially infinite) or limited (finite).
The pattern of arrivals can arrive according to
a known pattern or can arrive randomly.
Random arrivals generally follow a Poisson
distribution.
Copyright © 2012 Pearson Education 13-9
Characteristics of a Queuing System
Behavior of arrivals
Most queuing models assume customers are
patient and will wait in the queue until they are
served and do not switch lines.
Balking refers to customers who refuse to join
the queue.
Reneging customers enter the queue but
become impatient and leave without receiving
their service.
That these behaviors exist is a strong
argument for the use of queuing theory to
managing waiting lines.
Copyright © 2012 Pearson Education 13-10
Characteristics of a Queuing System
Waiting Line Characteristics
Waiting lines can be either limited or unlimited.
Queue discipline refers to the rule by which
customers in the line receive service.
The most common rule is first-in, first-out (FIFO).
Other rules are possible and may be based on other
important characteristics.
Other rules can be applied to select which
customers enter which queue, but may apply
FIFO once they are in the queue.
Copyright © 2012 Pearson Education 13-11
Characteristics of a Queuing System
Service Facility Characteristics
Basic queuing system configurations:
Service systems are classified in terms of
the number of channels, or servers, and the
number of phases, or service stops.
A single-channel system with one server is
quite common.
Multichannel systems exist when multiple
servers are fed by one common waiting line.
In a single-phase system, the customer
receives service form just one server.
In a multiphase system, the customer has to
go through more than one server.
Copyright © 2012 Pearson Education 13-12
Four basic
queuing
system
configurations
Figure 2
Copyright © 2012 Pearson Education 13-13
Characteristics of a Queuing System
Service time distribution
Service patterns can be either constant or
random.
Constant service times are often machine
controlled.
More often, service times are randomly
distributed according to a negative
exponential probability distribution.
Analysts should observe, collect, and plot
service time data to ensure that the
observations fit the assumed distributions
when applying these models.
Copyright © 2012 Pearson Education 13-14
Identifying Models Using
Kendall Notation
D. G. Kendall developed a notation for queuing
models that specifies the pattern of arrival, the
service time distribution, and the number of
channels.
Notation takes the form:
Arrival Service time Number of service
distribution distribution channels open
Specific letters are used to represent probability
distributions.
M = Poisson distribution for number of occurrences
D = constant (deterministic) rate
G = general distribution with known mean and variance
Copyright © 2012 Pearson Education 13-15
Identifying Models Using
Kendall Notation
A single-channel model with Poisson arrivals and
exponential service times would be represented
by:
M/M/1
If a second channel is added the notation would
read:
M/M/2
A three-channel system with Poisson arrivals and
constant service time would be
M/D/3
A four-channel system with Poisson arrivals and
normally distributed service times would be
M/G/4
Copyright © 2012 Pearson Education 13-16
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1)
Assumptions of the model:
Arrivals are served on a FIFO basis.
There is no balking or reneging.
Arrivals are independent of each other but the
arrival rate is constant over time.
Arrivals follow a Poisson distribution.
Service times are variable and independent but
the average is known.
Service times follow a negative exponential
distribution.
Average service rate is greater than the
average arrival rate.
Copyright © 2012 Pearson Education 13-17
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1)
When these assumptions are met, we
can develop a series of equations that
define the queue’s operating
characteristics.
Queuing Equations:
Let
= mean number of arrivals per time period
= mean number of customers or units
served per time period
The arrival rate and the service rate must be
defined for the same time period.
Copyright © 2012 Pearson Education 13-18
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1)
1. The average number of customers or units in the
system, LS :
LS
2. The average time a customer spends in the
system, WS :
1
WS
3. The average number of customers in the queue, Lq:
2
Lq
( )
Copyright © 2012 Pearson Education 13-19
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1)
4. The average time a customer spends waiting in
the queue, Wq :
Wq
( )
5. The utilization factor for the system, , the
probability the service facility is being used:
Copyright © 2012 Pearson Education 13-20
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1)
6. The percent idle time, P0, or the probability no one
is in the system:
P0 1
7. The probability that there are n number of
customers in the system :
Pn n (1 )
8. The probability that the number of customers in
the system is greater than k, Pn>k:
k 1
Pn k
Copyright © 2012 Pearson Education 13-21
Arnold’s Muffler Shop
Arnold’s mechanic can install mufflers at a rate of
3 per hour.
Customers arrive at a rate of 2 per hour.
So:
= 2 cars arriving per hour
= 3 cars serviced per hour
2 2
L 2 cars in the system
32 1 on average
1 1
W 1 hour that an average car
32 spends in the system
Copyright © 2012 Pearson Education 13-22
Arnold’s Muffler Shop
2 22 4
Lq 1.33 cars waiting in line
( ) 3(3 2) 3(1) on average
2
Wq hour 40 minutes average
( ) 3 waiting time per car
2
0.67 percentage of time
3 mechanic is busy
2
P0 1 1 0.33 probability that there
3 are 0 cars in the system
Copyright © 2012 Pearson Education 13-23
Arnold’s Muffler Shop
Probability of more than k cars in the system
k Pn>k = (2/3)k+1
0 0.667 Note that this is equal to 1 – P0 = 1 – 0.33 = 0.667
1 0.444
2 0.296
3 0.198 Implies that there is a 19.8% chance that more
than 3 cars are in the system
4 0.132
5 0.088
6 0.058
7 0.039
Copyright © 2012 Pearson Education 13-24
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1:N)
Assumptions of the model:
Arrivals are served on a FIFO basis.
There is no balking or reneging.
Arrivals are independent of each other but the
arrival rate is constant over time.
Arrivals follow a Poisson distribution.
Service times are variable and independent but
the average is known.
Service times follow a negative exponential
distribution.
Average service rate is greater than the
average arrival rate.
Capacity of system is limited to N customers
Copyright © 2012 Pearson Education 13-25
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1:N)
1. The average number of customers or units in the
system, LS :
N
( N 1) N 1
LS nPn
n 1 1 1 N 1
2. The average time a customer spends in the
system, WS :
LS
WS
(1 PN )
3. The average number of customers in the queue, Lq:
Lq LS
Copyright © 2012 Pearson Education 13-26
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1:N)
4. The average time a customer spends waiting in
the queue, Wq :
Lq
Wq
(1 PN )
5. The utilization factor for the system, , the
probability the service facility is being used:
Copyright © 2012 Pearson Education 13-27
Single-Channel Model, Poisson Arrivals,
Exponential Service Times (M/M/1:N)
6. The percent idle time, P0, or the probability no one
is in the system:
1
P0
1 N 1
7. The probability that there are n number of
customers in the system :
(1 ) n
Pn P n
1 N 1
0
7. The probability that the number of customers in
the system is greater than k, Pn>k:
k 1
Pn k
Copyright © 2012 Pearson Education 13-28
Chapter 15
Simulation Modeling
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Learning Objectives
After completing this chapter, students will be able to:
1. Tackle a wide variety of problems by
simulation
2. Understand the seven steps of conducting a
simulation
3. Explain the advantages and disadvantages of
simulation
4. Develop random number intervals and use
them to generate outcomes
5. Understand alternative simulation packages
available
© 2009 Prentice-Hall, Inc. 15 – 2
Chapter Outline
15.1 Introduction
15.2 Advantages and Disadvantages of Simulation
15.3 Monte Carlo Simulation
15.4 Simulation and Inventory Analysis
15.5 Simulation of a Queuing Problem
15.6 Fixed Time Increment and Next Event
Increment Simulation Models
15.7 Simulation Model for a Maintenance Policy
15.8 Two Other Types of Simulation
15.9 Verification and Validation
15.10 Role of Computers in Simulation
© 2009 Prentice-Hall, Inc. 15 – 3
Introduction
◼ Simulation is one of the most widely used
quantitative analysis tools
◼ It is used by over half of the largest US
corporations in corporate planning
◼ To simulate is to try to duplicate the features,
appearance, and characteristics of a real system
◼ We will build a mathematical model that comes
as close as possible to representing the reality
of the system
◼ You can also build physical models to test
systems
© 2009 Prentice-Hall, Inc. 15 – 4
Introduction
◼ The idea behind simulation is to imitate a real-
world situation mathematically
◼ Study its properties and operating characteristics
◼ Draw conclusions and make action decisions
© 2009 Prentice-Hall, Inc. 15 – 5
Introduction
◼ Using simulation, a manager should
1. Define a problem
2. Introduce the variables associated with the
problem
3. Construct a numerical model
4. Set up possible courses of action for testing
5. Run the experiment
6. Consider the results
7. Decide what courses of action to take
© 2009 Prentice-Hall, Inc. 15 – 6
Process of Simulation
Define Problem
Introduce Important
Variables
Construct Simulation
Model
Specify Values of
Variables to Be Tested
Conduct the
Simulation
Examine the
Results
Select Best Course
Figure 15.1 of Action
© 2009 Prentice-Hall, Inc. 15 – 7
Advantages and Disadvantages
of Simulation
◼ Simulation is useful because
1. It is relatively straightforward and flexible
2. Recent advances in computer software make
simulation models very easy to develop
3. Can be used to analyze large and complex
real-world situations
4. Allows “what-if?” type questions
5. Does not interfere with the real-world system
6. Enables study of interactions between
components
7. Enables time compression
8. Enables the inclusion of real-world
complications
© 2009 Prentice-Hall, Inc. 15 – 8
Advantages and Disadvantages
of Simulation
◼ The main disadvantages of simulation are
1. It is often expensive as it may require a long,
complicated process to develop the model
2. Does not generate optimal solutions, it is a
trial-and-error approach
3. Requires managers to generate all conditions
and constraints of real-world problem
4. Each model is unique and the solutions and
inferences are not usually transferable to
other problems
© 2009 Prentice-Hall, Inc. 15 – 9
Monte Carlo Simulation
◼ When systems contain elements that exhibit
chance in their behavior, the Monte Carlo method
of simulation can be applied
◼ Some examples are
1. Inventory demand
2. Lead time for inventory
3. Times between machine breakdowns
4. Times between arrivals
5. Service times
6. Times to complete project activities
7. Number of employees absent
© 2009 Prentice-Hall, Inc. 15 – 10
Monte Carlo Simulation
◼ The basis of the Monte Carlo simulation is
experimentation on the probabilistic elements
through random sampling
◼ It is based on the following five steps
1. Setting up a probability distribution for
important variables
2. Building a cumulative probability distribution
for each variable
3. Establishing an interval of random numbers
for each variable
4. Generating random numbers
5. Actually simulating a series of trials
© 2009 Prentice-Hall, Inc. 15 – 11
Harry’s Auto Tire Example
◼ A popular radial tire accounts for a large portion
of the sales at Harry’s Auto Tire
◼ Harry wishes to determine a policy for managing
this inventory
◼ He wants to simulate the daily demand for a
number of days
Step 1: Establishing probability distributions
◼ One way to establish a probability distribution for
a given variable is to examine historical outcomes
◼ Managerial estimates based on judgment and
experience can also be used
© 2009 Prentice-Hall, Inc. 15 – 12
Harry’s Auto Tire Example
◼ Historical daily demand for radial tires
DEMAND FOR TIRES FREQUENCY (DAYS)
0 10
1 20
2 40
3 60
4 40
5 30
200
Table 15.1
© 2009 Prentice-Hall, Inc. 15 – 13
Harry’s Auto Tire Example
Step 2: Building a cumulative probability distribution
for each variable
◼ Converting from a regular probability to a
cumulative distribution is an easy job
◼ A cumulative probability is the probability that a
variable will be less than or equal to a particular
value
◼ A cumulative distribution lists all of the possible
values and the probabilities
◼ Tables 15.2 and 15.3 show these distributions
© 2009 Prentice-Hall, Inc. 15 – 14
Harry’s Auto Tire Example
◼ Probability of demand for radial tires
DEMAND VARIABLE PROBABILITY OF OCCURRENCE
0 10/200 = 0.05
1 20/200 = 0.10
2 40/200 = 0.20
3 60/200 = 0.30
4 40/200 = 0.20
5 30/200 = 0.15
200/200 = 1.00
Table 15.2
© 2009 Prentice-Hall, Inc. 15 – 15
Harry’s Auto Tire Example
◼ Cumulative probability for radial tires
DAILY DEMAND PROBABILITY CUMULATIVE PROBABILITY
0 0.05 0.05
1 0.10 0.15
2 0.20 0.35
3 0.30 0.65
4 0.20 0.85
5 0.15 1.00
Table 15.3
© 2009 Prentice-Hall, Inc. 15 – 16
Harry’s Auto Tire Example
Step 3: Setting random number intervals
◼ We assign a set of numbers to represent each
possible value or outcome
◼ These are random number intervals
◼ A random number is a series of digits that have
been selected by a totally random process
◼ The range of the random number intervals
corresponds exactly to the probability of the
outcomes as shown in Figure 15.2
© 2009 Prentice-Hall, Inc. 15 – 17
Harry’s Auto Tire Example
◼ Graphical representation of the cumulative
probability
distribution
1.00
for radial 1.00 – – 00
tires 0.85
– 86
0.80 – 85
Represents 4
Cumulative Probability
Tires Demanded
0.65
– 66
0.60 – 65
Numbers
Random
0.40 – 0.35
– 36
35
0.20 – 0.15
– 16
15 Represents 1
0.05 06 Tire Demanded
–
– 05
Figure 15.2 0.00 – – 01
0 1 2 3 4 5
Daily Demand for Radials
© 2009 Prentice-Hall, Inc. 15 – 18
Harry’s Auto Tire Example
◼ Assignment of random number intervals for
Harry’s Auto Tire
DAILY DEMAND PROBABILITY CUMULATIVE INTERVAL OF
PROBABILITY RANDOM NUMBERS
0 0.05 0.05 01 to 05
1 0.10 0.15 06 to 15
2 0.20 0.35 16 to 35
3 0.30 0.65 36 to 65
4 0.20 0.85 66 to 85
5 0.15 1.00 86 to 00
Table 15.4
© 2009 Prentice-Hall, Inc. 15 – 19
Harry’s Auto Tire Example
Step 4: Generating random numbers
◼ Random numbers can be generated in several
ways
◼ Large problems will use computer program to
generate the needed random numbers
◼ For small problems, random processes like
roulette wheels or pulling chips from a hat may
be used
◼ The most common manual method is to use a
random number table
◼ Because everything is random in a random
number table, we can select numbers from
anywhere in the table to use in the simulation
© 2009 Prentice-Hall, Inc. 15 – 20
Harry’s Auto Tire Example
◼ Table of random numbers (partial)
52 06 50 88 53 30 10 47 99 37
37 63 28 02 74 35 24 03 29 60
82 57 68 28 05 94 03 11 27 79
69 02 36 49 71 99 32 10 75 21
98 94 90 36 06 78 23 67 89 85
96 52 62 87 49 56 59 23 78 71
33 69 27 21 11 60 95 89 68 48
50 33 50 95 13 44 34 62 64 39
88 32 18 50 62 57 34 56 62 31
90 30 36 24 69 82 51 74 30 35
Table 15.5
© 2009 Prentice-Hall, Inc. 15 – 21
Harry’s Auto Tire Example
Step 5: Simulating the experiment
◼ We select random numbers from Table 15.5
◼ The number we select will have a corresponding
range in Table 15.4
◼ We use the daily demand that corresponds to the
probability range aligned with the random
number
© 2009 Prentice-Hall, Inc. 15 – 22
Harry’s Auto Tire Example
◼ Ten-day simulation of demand for radial tires
DAY RANDOM NUMBER SIMULATED DAILY DEMAND
1 52 3
2 37 3
3 82 4
4 69 4
5 98 5
6 96 5
7 33 2
8 50 3
9 88 5
10 90 5
39 = total 10-day demand
3.9 = average daily demand for tires
Table 15.6
© 2009 Prentice-Hall, Inc. 15 – 23
Harry’s Auto Tire Example
◼ Note that the average demand from this
simulation (3.9 tires) is different from the
expected daily demand
Expected 5
daily = (Probability of i tires)(Demand of i tires)
demand i =0
= (0.05)(0) + (0.10)(1) + (0.20)(2) + (0.30)(3)
+ (0.20)(4) + (0.15)(5)
= 2.95 tires
◼ If this simulation were repeated hundreds or
thousands of times it is much more likely the
average simulated demand would be nearly the
same as the expected demand
© 2009 Prentice-Hall, Inc. 15 – 24
Using QM for Windows for Simulation
◼ Monte Carlo simulation of Harry’s Auto Tire using
QM for Windows
Program 15.1
© 2009 Prentice-Hall, Inc. 15 – 25
Simulation with Excel Spreadsheets
◼ Using Excel to simulate tire demand for Harry’s
Auto Tire Shop
Program 15.2A
© 2009 Prentice-Hall, Inc. 15 – 26
Simulation with Excel Spreadsheets
◼ Excel simulation results for Harry’s Auto Tire
Shop showing a simulated average demand of 2.8
tires per day
Program 15.2B
© 2009 Prentice-Hall, Inc. 15 – 27
Simulation with Excel Spreadsheets
◼ Generating normal random numbers in Excel
Program 15.3A
© 2009 Prentice-Hall, Inc. 15 – 28
Simulation with Excel Spreadsheets
◼ Excel output with normal random numbers
Program 15.3B
© 2009 Prentice-Hall, Inc. 15 – 29
Simulation and Inventory Analysis
◼ We have seen deterministic inventory models
◼ In many real-world inventory situations, demand
and lead time are variables
◼ Accurate analysis is difficult without simulation
◼ We will look at an inventory problem with two
decision variables and two probabilistic
components
◼ The owner of a hardware store wants to establish
order quantity and reorder point decisions for a
product that has probabilistic daily demand and
reorder lead time
© 2009 Prentice-Hall, Inc. 15 – 30
Simkin’s Hardware Store
◼ The owner of a hardware store wants to find a
good, low cost inventory policy for an electric
drill
◼ Simkin identifies two types of variables,
controllable and uncontrollable inputs
◼ The controllable inputs are the order quantity and
reorder points
◼ The uncontrollable inputs are daily demand and
variable lead time
◼ The demand data for the drill is shown in Table
15.7
© 2009 Prentice-Hall, Inc. 15 – 31
Simkin’s Hardware Store
◼ Probabilities and random number intervals for
daily Ace drill demand
(1) (2) (4) (5)
DEMAND FOR FREQUENCY (3) CUMULATIVE INTERVAL OF
ACE DRILL (DAYS) PROBABILITY PROBABILITY RANDOM NUMBERS
0 15 0.05 0.05 01 to 05
1 30 0.10 0.15 06 to 15
2 60 0.20 0.35 16 to 35
3 120 0.40 0.75 36 to 75
4 45 0.15 0.90 76 to 90
5 30 0.10 1.00 91 to 00
300 1.00
Table 15.7
© 2009 Prentice-Hall, Inc. 15 – 32
Simkin’s Hardware Store
◼ Probabilities and random number intervals for
reorder lead time
(1) (2) (4) (5)
LEAD TIME FREQUENCY (3) CUMULATIVE RANDOM NUMBER
(DAYS) (ORDERS) PROBABILITY PROBABILITY INTERVAL
1 10 0.20 0.20 01 to 20
2 25 0.50 0.70 21 to 70
3 15 0.30 1.00 71 to 00
50 1.00
Table 15.8
© 2009 Prentice-Hall, Inc. 15 – 33
Simkin’s Hardware Store
◼ The third step is to develop a simulation model
◼ A flow diagram, or flowchart, is helpful in this
process
◼ The fourth step in the process is to specify the
values of the variables that we wish to test
◼ The first policy Simkin wants to test is an order
quantity of 10 with a reorder point of 5
◼ The fifth step is to actually conduct the
simulation
◼ The process is simulated for a 10 day period
© 2009 Prentice-Hall, Inc. 15 – 34
Simkin’s Hardware Store
◼ Flow diagram
for Simkin’s Start
inventory Begin day
example of simulation
Has Increase current
order Yes
inventory by
arrived? quantity ordered
No
Select random number
to generate today’s
demand
Is
demand greater Yes Record
than beginning number of
inventory lost sales
?
No
Figure 15.3 (a)
© 2009 Prentice-Hall, Inc. 15 – 35
Simkin’s Hardware Store
◼ Flow diagram
for Simkin’s No
inventory Compute ending inventory
= Beginning inventory
Record ending
inventory = 0
example – Demand
Is Has
Yes No
ending inventory order been Place
less than reorder placed that hasn’t order
point? arrived yet
?
No Have Yes
Select random
No enough days number to
of this order policy generate lead
been simulated time
?
Yes
Compute average ending inventory,
average lost sales, average number
of orders placed, and corresponding End
costs
Figure 15.3 (b)
© 2009 Prentice-Hall, Inc. 15 – 36
Simkin’s Hardware Store
◼ Using the table of random numbers, the
simulation is conducted using a four-step process
1. Begin each day by checking whether an ordered
inventory has arrived. If it has, increase the current
inventory by the quantity ordered.
2. Generate a daily demand from the demand probability
by selecting a random number
3. Compute the ending inventory every day. If on-hand
inventory is insufficient to meet the day’s demand,
satisfy as much as possible and note the number of
lost sales.
4. Determine whether the day’s ending inventory has
reached the reorder point. If necessary place an order.
© 2009 Prentice-Hall, Inc. 15 – 37
Simkin’s Hardware Store
◼ Simkin Hardware’s first inventory simulation
ORDER QUANTITY = 10 UNITS REORDER POINT = 5 UNITS
(2) (3) (4) (6) (7) (9) (10)
(1) UNITS BEGINNING RANDOM (5) ENDING LOST (8) RANDOM LEAD
DAY RECEIVED INVENTORY NUMBER DEMAND INVENTORY SALES ORDER NUMBER TIME
1 … 10 06 1 9 0 No
2 0 9 63 3 6 0 No
3 0 6 57 3 3 0 Yes 02 1
4 0 3 94 5 0 2 No
5 10 10 52 3 7 0 No
6 0 7 69 3 4 0 Yes 33 2
7 0 4 32 2 2 0 No
8 0 2 30 2 0 0 No
9 10 10 48 3 7 0 No
10 0 7 88 4 3 0 Yes 14 1
Total 41 2
Table 15.9
© 2009 Prentice-Hall, Inc. 15 – 38
Analyzing Simkin’s Inventory Cost
◼ The objective is to find a low-cost solution so
Simkin must determine what the costs are
◼ Equations for average daily ending inventory,
average lost sales, and average number of orders
placed
Average 41 total units
ending = = 4.1 units per day
inventory 10 days
Average 2 sales lost
= = 0.2 units per day
lost sales 10 days
Average 3 orders
number of = = 0.3 orders per day
orders placed 10 days
© 2009 Prentice-Hall, Inc. 15 – 39
Analyzing Simkin’s Inventory Cost
◼ Simkin’s store is open 200 days a year
◼ Estimated ordering cost is $10 per order
◼ Holding cost is $6 per drill per year
◼ Lost sales cost $8
Daily order cost = (Cost of placing one order)
x (Number of orders placed per day)
= $10 per order x 0.3 orders per day = $3
Daily holding cost = (Cost of holding one unit for one day) x
(Average ending inventory)
= $0.03 per unit per day x 4.1 units per day
= $0.12
© 2009 Prentice-Hall, Inc. 15 – 40
Analyzing Simkin’s Inventory Cost
◼ Simkin’s store is open 200 days a year
◼ Estimated ordering cost is $10 per order
◼ Holding cost is $6 per drill per year
◼ Lost sales cost $8
Daily stockout cost = (Cost per lost sale)
x (Average number of lost sales per
day)
= $8 per lost sale x 0.2 lost sales per day
= $1.60
Total daily
inventory cost = Daily order cost + Daily holding cost
+ Daily stockout cost
= $4.72
© 2009 Prentice-Hall, Inc. 15 – 41
Analyzing Simkin’s Inventory Cost
◼ For the year, this policy would cost approximately
$944
◼ This simulation should really be extended for
many more days, perhaps 100 or 1,000 days
◼ Even after a larger simulation, the model must be
verified and validated to make sure it truly
represents the situation on which it is based
◼ If we are satisfied with the model, additional
simulations can be conducted using other values
for the variables
◼ After simulating all reasonable combinations,
Simkin would select the policy that results in the
lowest total cost
© 2009 Prentice-Hall, Inc. 15 – 42
Simulation of a Queuing Problem
◼ Modeling waiting lines is an important application
of simulation
◼ The assumptions of queuing models are quite
restrictive
◼ Sometimes simulation is the only approach that
fits
◼ In this example, arrivals do not follow a Poisson
distribution and unloading rates are not
exponential or constant
© 2009 Prentice-Hall, Inc. 15 – 43
Port of New Orleans
◼ Fully loaded barges arrive at night for unloading
◼ The number of barges each night varies from 0 - 5
◼ The number of barges vary from day to day
◼ The supervisor has information which can be
used to create a probability distribution for the
daily unloading rate
◼ Barges are unloaded first-in, first-out
◼ Barges must wait for unloading which is
expensive
◼ The dock superintendent wants to do a simulation
study to enable him to make better staffing
decisions
© 2009 Prentice-Hall, Inc. 15 – 44
Port of New Orleans
◼ Overnight barge arrival rates and random number
intervals
NUMBER OF CUMULATIVE RANDOM
ARRIVALS PROBABILITY PROBABILITY NUMBER INTERVAL
0 0.13 0.13 01 to 13
1 0.17 0.30 14 to 30
2 0.15 0.45 31 to 45
3 0.25 0.70 46 to 70
4 0.20 0.90 71 to 90
5 0.10 1.00 91 to 00
Table 15.10
© 2009 Prentice-Hall, Inc. 15 – 45
Port of New Orleans
◼ Unloading rates and random number intervals
DAILY UNLOADING CUMULATIVE RANDOM
RATE PROBABILITY PROBABILITY NUMBER INTERVAL
1 0.05 0.05 01 to 05
2 0.15 0.20 06 to 20
3 0.50 0.70 21 to 70
4 0.20 0.90 71 to 90
5 0.10 1.00 91 to 00
1.00
Table 15.11
© 2009 Prentice-Hall, Inc. 15 – 46
Port of New Orleans
◼ Queuing simulation of barge unloadings
(1) (2) (3) (4) (5) (6) (7)
DAY NUMBER DELAYED RANDOM NUMBER OF TOTAL TO BE RANDOM NUMBER
FROM PREVIOUS DAY NUMBER NIGHTLY ARRIVALS UNLOADED NUMBER UNLOADED
1 — 52 3 3 37 3
2 0 06 0 0 63 0
3 0 50 3 3 28 3
4 0 88 4 4 02 1
5 3 53 3 6 74 4
6 2 30 1 3 35 3
7 0 10 0 0 24 0
8 0 47 3 3 03 1
9 2 99 5 7 29 3
10 4 37 2 6 60 3
11 3 66 3 6 74 4
12 2 91 5 7 85 4
13 3 35 2 5 90 4
14 1 32 2 3 73 3
15 0 00 5 5 59 3
20 41 39
Total delays Total arrivals Total unloadings
Table 15.12 © 2009 Prentice-Hall, Inc. 15 – 47
Port of New Orleans
◼ Three important pieces of information
Average number of barges 20 delays
=
delayed to the next day 15 days
= 1.33 barges delayed per day
Average number of 41 arrivals
= = 2.73 arrivals
nightly arrivals 15 days
Average number of barges 39 unloadings
unloaded each day = = 2.60 unloadings
15 days
© 2009 Prentice-Hall, Inc. 15 – 48
Using Excel to Simulate the Port of
New Orleans Queuing Problem
◼ An Excel model for the Port of New Orleans
queuing simulation
Program 15.4A
© 2009 Prentice-Hall, Inc. 15 – 49
Using Excel to Simulate the Port of
New Orleans Queuing Problem
◼ Output from the Excel formulas in Program 15.4A
Program 15.4B
© 2009 Prentice-Hall, Inc. 15 – 50
Fixed Time Increment and Next Event
Increment Simulation Models
◼ Simulation models are often classified into fixed
time increment models and next event increment
models
◼ The terms refer to the frequency in which the
system status is updated
◼ Fixed time increments update the status of the
system at fixed time intervals
◼ Next event increment models update only when
the system status changes
◼ Fixed event models randomly generate the
number of events that occur during a time period
◼ Next event models randomly generate the time
that elapses until the next event occurs
© 2009 Prentice-Hall, Inc. 15 – 51
Simulation Model for a
Maintenance Policy
◼ Simulation can be used to analyze different
maintenance policies before actually
implementing them
◼ Many options regarding staffing levels, parts
replacement schedules, downtime, and labor
costs can be compared
◼ This can including completely shutting down
factories for maintenance
© 2009 Prentice-Hall, Inc. 15 – 52
Three Hills Power Company
◼ Three Hills provides power to a large city through
a series of almost 200 electric generators
◼ The company is concerned about generator
failures because a breakdown costs about $75
per generator per hour
◼ Their four repair people earn $30 per hour and
work rotating 8 hour shifts
◼ Management wants to evaluate the
1. Service maintenance cost
2. Simulated machine breakdown cost
3. Total cost
© 2009 Prentice-Hall, Inc. 15 – 53
Three Hills Power Company
◼ There are two important maintenance system
components
◼ Time between successive generator breakdowns
which varies from 30 minutes to three hours
◼ The time it takes to repair the generators which
ranges from one to three hours in one hour
blocks
◼ A next event simulation is constructed to study
this problem
© 2009 Prentice-Hall, Inc. 15 – 54
Three Hills Power Company
◼ Three Hills Start
flow diagram
Generate random number
for “Time Between
Breakdowns”
Record actual clock time
of breakdown
Examine time previous
repair ends
Is the
repairperson No Wait until previous
free to begin repair is completed
repair?
Figure 15.4 (a)
Yes
© 2009 Prentice-Hall, Inc. 15 – 55
Three Hills Power Company
◼ Three Hills Generate random number
for repair time required
flow diagram
Compute time repair
completed
Compute hours of machine
downtime = Time repair
completed – Clock time
of breakdown
No Enough
breakdowns
simulated?
Yes
Figure 15.4 (b) Compute downtime and
comparative cost data End
© 2009 Prentice-Hall, Inc. 15 – 56
Three Hills Power Company
◼ Time between generator breakdowns at Three
Hills Power
TIME BETWEEN
RECORDED NUMBER RANDOM
MACHINE OF TIMES CUMULATIVE NUMBER
FAILURES (HRS) OBSERVED PROBABILITY PROBABILITY INTERVAL
0.5 5 0.05 0.05 01 to 05
1.0 6 0.06 0.11 06 to 11
1.5 16 0.16 0.27 12 to 27
2.0 33 0.33 0.60 28 to 60
2.5 21 0.21 0.81 61 to 81
3.0 19 0.19 1.00 82 to 00
Total 100 1.00
Table 15.13
© 2009 Prentice-Hall, Inc. 15 – 57
Three Hills Power Company
◼ Generator repair times required
NUMBER RANDOM
REPAIR TIME OF TIMES CUMULATIVE NUMBER
REQUIRED (HRS) OBSERVED PROBABILITY PROBABILITY INTERVAL
1 28 0.28 0.28 01 to 28
2 52 0.52 0.80 29 to 80
3 20 0.20 1.00 81 to 00
Total 100 1.00
Table 15.14
© 2009 Prentice-Hall, Inc. 15 – 58
Three Hills Power Company
◼ Simulation of generator breakdowns and repairs
(5) (6) (9)
TIME REPAIR- RANDOM NUMBER
(2) (3) PERSON IS NUMBER (7) (8) OF
(1) RANDOM TIME (4) FREE TO FOR REPAIR TIME HOURS
BREAKDOWN NUMBER FOR BETWEEN TIME OF BEGIN THIS REPAIR TIME REPAIR MACHINE
NUMBER BREAKDOWNS BREAKDOWNS BREAKDOWN REPAIR TIME REQUIRED ENDS DOWN
1 57 2 02:00 02:00 07 1 03:00 1
2 17 1.5 03:30 03:30 60 2 05:30 2
3 36 2 05:30 05:30 77 2 07:30 2
4 72 2.5 08:00 08:00 49 2 10:00 2
5 85 3 11:00 11:00 76 2 13:00 2
6 31 2 13:00 13:00 95 3 16:00 3
7 44 2 15:00 16:00 51 2 18:00 3
8 30 2 17:00 18:00 16 1 19:00 2
9 26 1.5 18:30 19:00 14 1 20:00 1.5
10 09 1 19:30 20:00 85 3 23:00 3.5
11 49 2 21:30 23:00 59 2 01:00 3.5
12 13 1.5 23:00 01:00 85 3 04:00 5
13 33 2 01:00 04:00 40 2 06:00 5
14 89 3 04:00 06:00 42 2 08:00 4
15 13 1.5 05:30 08:00 52 2 10:00 4.5
Total 44
Table 15.15
© 2009 Prentice-Hall, Inc. 15 – 59
Cost Analysis of Simulation
◼ The simulation of 15 generator breakdowns
covers 34 hours of operation
◼ The analysis of this simulation is
Service
maintenance = 34 hours of worker service time
cost x $30 per hour
= $1,020
Simulated machine
breakdown cost = 44 total hours of breakdown
x $75 lost per hour of downtime
= $3,300
Total simulated
maintenance cost of = Service cost + Breakdown cost
the current system = $1,020 + $3,300
= $4,320
© 2009 Prentice-Hall, Inc. 15 – 60
Cost Analysis of Simulation
◼ The cost of $4,320 should be compared with
other alternative plans to see if this is a “good”
value
◼ The company might explore options like adding
another repairperson
◼ Strategies such as preventive maintenance might
also be simulated for comparison
© 2009 Prentice-Hall, Inc. 15 – 61
Building an Excel Simulation Model
for Three Hills Power Company
◼ An Excel spreadsheet model for simulating the
Three Hills Power Company maintenance problem
Program 15.5A
© 2009 Prentice-Hall, Inc. 15 – 62
Building an Excel Simulation Model
for Three Hills Power Company
◼ Output from Excel spreadsheet in Program 15.5A
Program 15.5B
© 2009 Prentice-Hall, Inc. 15 – 63
Two Other Types of
Simulation Models
◼ Simulation models are often broken into
three categories
◼ The Monte Carlo method
◼ Operational gaming
◼ Systems simulation
◼ Though theoretically different,
computerized simulation has tended to
blur the differences
© 2009 Prentice-Hall, Inc. 15 – 64
Operational Gaming
◼ Operational gaming refers to simulation involving
two or more competing players
◼ The best examples of this are military games and
business games
◼ These types of simulation allow the testing of
skills and decision-making in a competitive
environment
© 2009 Prentice-Hall, Inc. 15 – 65
Systems Simulation
◼ Systems simulation is similar in that allows users
to test various managerial policies and decisions
to evaluate their effect on the operating
environment
◼ This models the dynamics of large systems
◼ A corporate operating system might model sales,
production levels, marketing policies,
investments, union contracts, utility rates,
financing, and other factors
◼ Economic simulations, often called econometric
models, are used by governments, bankers, and
large organizations to predict inflation rates,
domestic and foreign money supplies, and
unemployment levels
© 2009 Prentice-Hall, Inc. 15 – 66
Systems Simulation
◼ Inputs and outputs of a typical economic system
simulation
Income Tax Gross National
Levels Product
Corporate Tax Inflation Rates
Rates Econometric Model Unemployment
Interest Rates (in Series of Rates
Mathematical
Government Equations) Monetary
Spending Supplies
Foreign Trade Population
Policy Growth Rates
Figure 15.5
© 2009 Prentice-Hall, Inc. 15 – 67
Verification and Validation
◼ It is important that a simulation model be checked
to see that it is working properly and providing
good representation of the real world situation
◼ The verification process involves determining
that the computer model is internally consistent
and following the logic of the conceptual model
◼ Verification answers the question “Did we build
the model right?”
◼ Validation is the process of comparing a
simulation model to the real system it represents
to make sure it is accurate
◼ Validation answers the question “Did we build the
right model?”
© 2009 Prentice-Hall, Inc. 15 – 68
Role of Computers in Simulation
◼ Computers are critical in simulating complex
tasks
◼ Three types of computer programming languages
are available to help the simulation process
◼ General-purpose languages
◼ Special-purpose simulation languages
1. These require less programming
2. Are more efficient and easier to check for errors
3. Have random number generators built in
◼ Pre-written simulation programs built to handle a
wide range of common problems
◼ Excel and add-ins can also be used for
simulation problems
© 2009 Prentice-Hall, Inc. 15 – 69
Chapter 13
Project Management
To accompany
Quantitative Analysis for Management, Tenth Edition,
by Render, Stair, and Hanna © 2008 Prentice-Hall, Inc.
Power Point slides created by Jeff Heyl © 2009 Prentice-Hall, Inc.
Learning Objectives
After completing this chapter, students will be able to:
1. Understand how to plan, monitor, and control
projects with the use of PERT and CPM
2. Determine earliest start, earliest finish, latest
start, latest finish, and slack times for each
activity, along with the total project
completion time
3. Reduce total project time at the least total
cost by crashing the network using manual
or linear programming techniques
4. Understand the important role of software in
project management
© 2009 Prentice-Hall, Inc. 13 – 2
Chapter Outline
13.1 Introduction
13.2 PERT/CPM
13.3 PERT/Cost
13.4 Project Crashing
13.5 Other Topics in Project
Management
© 2009 Prentice-Hall, Inc. 13 – 3
Introduction
Most realistic projects are large and complex
Tens of thousands of steps and millions of dollars
may be involved
Managing large-scale, complicated projects
effectively is a difficult problem and the stakes are
high
The first step in planning and scheduling a project
is to develop the work breakdown structure
Time, cost, resource requirements, predecessors,
and people required are identified for each activity
Then a schedule for the project can be developed
© 2009 Prentice-Hall, Inc. 13 – 4
Project Management
It is nothing more (or less) than knowing
what the status of a project is:
when it should be done
how much (and if) it has slipped from the
original schedule
what the bottlenecks are
what you might drop to save some time
© 2009 Prentice-Hall, Inc. 13 – 5
Project Management Models
History
One of the earliest techniques was the Gantt chart
(used by US Navy).
This type of chart shows the start and finish times
of one or more activities, as shown below:
© 2009 Prentice-Hall, Inc. 13 – 6
Project Planning, Controlling
and Scheduling
Project Planning:
1. Setting goals.
2. Defining the project.
3. Tying needs into timed project activities.
4. Organizing the team.
Project Scheduling:
1. Tying resources to specific activities.
2. Relating activities to each other.
3. Updating and revising on regular basis.
Before Project Project Controlling:
1. Monitoring resources, costs, quality
During Project and budgets.
2. Revising and changing plans.
3. Shifting resources to meet demands.
© 2009 Prentice-Hall, Inc. 13 – 7
Project Management Models
PERT
PERT/Cost
Critical Path Method (CPM)
© 2009 Prentice-Hall, Inc. 13 – 8
Introduction
The program evaluation and review technique
(PERT)
PERT and the critical path method (CPM)
CPM are two
popular quantitative analysis techniques to help
plan, schedule, monitor, and control projects
They were developed because there was a critical
need for a better way to manage.
Originally the approaches differed in how they
estimated activity times
PERT used three time estimates to develop a
probabilistic estimate of completion time
CPM was a more deterministic technique
They have become so similar they are commonly
considered one technique, PERT/CPM
© 2009 Prentice-Hall, Inc. 13 – 9
Six Steps of PERT/CPM
1. Define the project and all of its significant
activities or tasks
2. Develop the relationships among the activities
and decide which activities must precede others
3. Draw the network connecting all of the activities
4. Assign time and/or cost estimates to each
activity
5. Compute the longest time path through the
network; this is called the critical path
6. Use the network to help plan, schedule, monitor,
and control the project
The critical path is important since any delay in
these activities can delay the completion of the
project
© 2009 Prentice-Hall, Inc. 13 – 10
PERT/CPM
Given the large number of tasks in a project,
it is easy to see why the following questions
are important
1. When will the entire project be completed?
2. What are the critical activities or tasks in the
project, that is, the ones that will delay the
entire project if they are late?
3. Which are the non-critical activities, that is,
the ones that can run late without delaying
the entire project’s completion?
4. If there are three time estimates, what is the
probability that the project will be completed
by a specific date?
© 2009 Prentice-Hall, Inc. 13 – 11
PERT/CPM
5. At any particular date, is the project on
schedule, behind schedule, or ahead of
schedule?
6. On any given date, is the money spent equal
to, less than, or greater than the budgeted
amount?
7. Are there enough resources available to
finish the project on time?
8. If the project is to be finished in a shorter
amount of time, what is the best way to
accomplish this at the least cost?
© 2009 Prentice-Hall, Inc. 13 – 12
General Foundry Example of
PERT/CPM
General Foundry, Inc. has long been trying to
avoid the expense of installing air pollution
control equipment
The local environmental protection group has
recently given the foundry 16 weeks to install a
complex air filter system on its main smokestack
General Foundry was warned that it will be forced
to close unless the device is installed in the
allotted period
They want to make sure that installation of the
filtering system progresses smoothly and on time
© 2009 Prentice-Hall, Inc. 13 – 13
General Foundry Example of
PERT/CPM
Activities and immediate predecessors for
General Foundry
IMMEDIATE
ACTIVITY DESCRIPTION PREDECESSORS
A Build internal components —
B Modify roof and floor —
C Construct collection stack A
D Pour concrete and install frame B
E Build high-temperature burner C
F Install control system C
G Install air pollution device D, E
H Inspect and test F, G
Table 13.1
© 2009 Prentice-Hall, Inc. 13 – 14
Drawing the PERT/CPM Network
There are two common techniques for drawing
PERT networks
Activity-on-node (AON)
AON where the nodes
represent activities
Activity-on-arc (AOA)
AOA where the arcs are used to
represent the activities
The AON approach is easier and more commonly
found in software packages
One node represents the start of the project, one
node for the end of the project, and nodes for
each of the activities
The arcs are used to show the predecessors for
each activity
© 2009 Prentice-Hall, Inc. 13 – 15
General Foundry Example of
PERT/CPM
Network for General Foundry
A C F
Build Internal Construct Install Control
Components Collection Stack System
E H
Start Inspect Finish
Build Burner and Test
B D G
Modify Roof Pour Concrete Install Pollution
and Floor and Install Frame Device
Figure 13.1
© 2009 Prentice-Hall, Inc. 13 – 16
Activity Times
In some situations, activity times are known with
certainty
CPM assigns just one time estimate to each
activity and this is used to find the critical path
In many projects there is uncertainty about activity
times
PERT employs a probability distribution based on
three time estimates for each activity
A weighted average of these estimates is used for
the time estimate and this is used to determine the
critical path
© 2009 Prentice-Hall, Inc. 13 – 17
Activity Times
The time estimates in PERT are
Optimistic time (a) = time an
activity will take if everything goes
as well as possible. There should be
only a small probability (say, 1/100) of
this occurring.
Pessimistic time (b) = time an
activity would take assuming very
unfavorable conditions. There
should also be only a small
probability that the activity will really
take this long.
Most likely time (m) = most
realistic time estimate to complete
the activity
© 2009 Prentice-Hall, Inc. 13 – 18
Activity Times
PERT often assumes time estimates follow a beta probability
distribution
The beta probability distribution is often used when there is no solid
historical data upon which to activity time base estimates
Found to be appropriate in many cases for determining an expected value
and variance for activity completion times
Probability of 1 in 100
of a Occurring
Probability
Probability of 1 in 100
of b Occurring
Activity Time
Most Most Most
Optimistic Likely Pessimistic
Time Time Time
(a) (m) (b)
Figure 13.2
© 2009 Prentice-Hall, Inc. 13 – 19
Activity Times
To find the expected activity time (t), the beta
distribution weights the estimates as follows
a + 4m + b
t=
6
To compute the dispersion or variance of activity
completion time,
time we use the formula
2
b−a
Variance =
6
© 2009 Prentice-Hall, Inc. 13 – 20
Activity Times
Time estimates (weeks) for General Foundry
MOST EXPECTED VARIANCE,
OPTIMISTIC, PROBABLE, PESSIMISTIC, TIME,
ACTIVITY a m b t = [(a + 4m + b)/6] [(b – a)/6]2
A 1 2 3 2 4/36
B 2 3 4 3 4/36
C 1 2 3 2 4/36
D 2 4 6 4 16/36
E 1 4 7 4 36/36
F 1 2 9 3 64/36
G 3 4 11 5 64/36
H 1 2 3 2 4/36
25
Table 13.2
© 2009 Prentice-Hall, Inc. 13 – 21
How to Find the Critical Path
We accept the expected completion time for
each task as the actual time for now
The total of 25 weeks in Table 13.2 does not take
into account the obvious fact that some of the
tasks could be taking place at the same time
To find out how long the project will take we
perform the critical path analysis for the network
The critical path is the longest path through the
network
© 2009 Prentice-Hall, Inc. 13 – 22
How to Find the Critical Path
General Foundry’s network with expected activity
times
A 2 C 2 F 3
E 4 H 2
Start Finish
B 3 D 4 G 5
Figure 13.3
© 2009 Prentice-Hall, Inc. 13 – 23
How to Find the Critical Path
To find the critical path, need to determine the
following quantities for each activity in the
network
1. Earliest start time (ES):
ES the earliest time an
activity can begin without violation of immediate
predecessor requirements
2. Earliest finish time (EF):
EF the earliest time at
which an activity can end
3. Latest start time (LS):
LS the latest time an activity
can begin without delaying the entire project
4. Latest finish time (LF):
LF the latest time an activity
can end without delaying the entire project
© 2009 Prentice-Hall, Inc. 13 – 24
How to Find the Critical Path
In the nodes, the activity time and the early and
late start and finish times are represented in the
following manner
ACTIVITY t
ES EF
LS LF
Earliest times are computed as
Earliest finish time =
Earliest start time
+ Expected activity time
EF = ES + t
Earliest start = Largest of the earliest finish times of
immediate predecessors
ES = Largest EF of immediate predecessors
© 2009 Prentice-Hall, Inc. 13 – 25
How to Find the Critical Path
At the start of the project we set the time to zero
Thus ES = 0 for both A and B
A t=2
ES = 0 EF = 0 + 2 = 2
Start
B t=3
ES = 0 EF = 0 + 3 = 3
© 2009 Prentice-Hall, Inc. 13 – 26
How to Find the Critical Path
General Foundry’s ES and EF times
A 2 C 2 F 3
0 2 2 4 4 7
E 4 H 2
Start 4 8 13 15 Finish
B 3 D 4 G 5
0 3 3 7 8 13
Figure 13.4
© 2009 Prentice-Hall, Inc. 13 – 27
How to Find the Critical Path
Latest times are computed as
Latest start time =
Latest finish time
– Expected activity time
LS = LF – t
Latest finish time = Smallest of latest start times
for following activities
LF = Smallest LS of following activities
For activity H
LS = LF – t = 15 – 2 = 13 weeks
© 2009 Prentice-Hall, Inc. 13 – 28
How to Find the Critical Path
General Foundry’s LS and LF times
A 2 C 2 F 3
0 2 2 4 4 7
0 2 2 4 10 13
E 4 H 2
Start 4 8 13 15 Finish
4 8 13 15
B 3 D 4 G 5
0 3 3 7 8 13
1 4 4 8 8 13
Figure 13.5
© 2009 Prentice-Hall, Inc. 13 – 29
How to Find the Critical Path
Once ES, LS, EF, and LF have been determined, it
is a simple matter to find the amount of slack time
that each activity has
Slack = LS – ES, or Slack = LF – EF
From Table 13.3 we see activities A, C, E, G, and H
have no slack time
These are called critical activities and they are said
to be on the critical path
The total project completion time is 15 weeks
Industrial managers call this a boundary timetable
© 2009 Prentice-Hall, Inc. 13 – 30
How to Find the Critical Path
General Foundry’s schedule and slack times
EARLIEST EARLIEST LATEST LATEST ON
START, FINISH, START, FINISH, SLACK, CRITICAL
ACTIVITY ES EF LS LF LS – ES PATH?
A 0 2 0 2 0 Yes
B 0 3 1 4 1 No
C 2 4 2 4 0 Yes
D 3 7 4 8 1 No
E 4 8 4 8 0 Yes
F 4 7 10 13 6 No
G 8 13 8 13 0 Yes
H 13 15 13 15 0 Yes
Table 13.3
© 2009 Prentice-Hall, Inc. 13 – 31
How to Find the Critical Path
General Foundry’s critical path
A 2 C 2 F 3
0 2 2 4 4 7
0 2 2 4 10 13
E 4 H 2
Start 4 8 13 15 Finish
4 8 13 15
B 3 D 4 G 5
0 3 3 7 8 13
1 4 4 8 8 13
Figure 13.6
© 2009 Prentice-Hall, Inc. 13 – 32
Probability of Project Completion
The critical path analysis helped determine the
expected project completion time of 15 weeks
But variation in activities on the critical path can
affect overall project completion, and this is a
major concern
If the project is not complete in 16 weeks, the
foundry will have to close
PERT uses the variance of critical path activities
to help determine the variance of the overall
project
Project variance = ∑ variances of activities
on the critical path
© 2009 Prentice-Hall, Inc. 13 – 33
Probability of Project Completion
From Table 13.2 we know that
ACTIVITY VARIANCE
A 4/36
B 4/36
C 4/36
D 16/36
E 36/36
F 64/36
G 64/36
H 4/36
Hence, the project variance is
Project variance = 4/36 + 4/36 + 36/36 + 64/36 + 4/36 = 112/36 = 3.111
© 2009 Prentice-Hall, Inc. 13 – 34
Probability of Project Completion
We know the standard deviation is just the
square root of the variance, so
Project standard deviation = σ T = Project variance
= 3.11 = 1.76 weeks
We assume activity times are independent and
total project completion time is normally
distributed
© 2009 Prentice-Hall, Inc. 13 – 35
Probability of Project Completion
The project’s expected completion date is 15 weeks.
Assume that the total project completion time follows a
normal probability distribution
Chart tells us that there is a 50% chance of completing the
entire project in less than 15 weeks and a 50% chance it will
exceed 15 weeks
Standard Deviation = 1.76 Weeks
15 Weeks
(Expected Completion Time)
Figure 13.7
© 2009 Prentice-Hall, Inc. 13 – 36
Probability of Project Completion
The standard normal equation can be applied as
follows
Due date − Expected date of completion
Z=
σT
16 weeks − 15 weeks
= = 0.57
1.76 weeks
From Appendix A we find the probability of
0.71566 associated with this Z value
That means there is a 71.6% probability this
project can be completed in 16 weeks or less
© 2009 Prentice-Hall, Inc. 13 – 37
Probability of Project Completion
Probability of General Foundry meeting the 16-
week deadline
Expected Time is 15 Weeks 0.57 Standard Deviations
Probability
(T ≤ 16 Weeks)
is 71.6%
15 16
Time
Weeks Weeks
Figure 13.8
© 2009 Prentice-Hall, Inc. 13 – 38
What PERT Was Able to Provide
PERT has been able to provide the project manager
with several valuable pieces of information
The project’s expected completion date is 15 weeks
There is a 71.6% chance that the equipment will be in
place within the 16-week deadline
Five activities (A, C, E, G, H) are on the critical path
If any one of the critical activities is delayed for any
reason, the entire project will be delayed.
Three activities (B, D, F) are not critical but have some
slack time built in
They can borrow from their resources, if needed, possibly
to speed up the entire project.
A detailed schedule of activity starting and ending
dates has been made available
© 2009 Prentice-Hall, Inc. 13 – 39
Sensitivity Analysis and
Project Management
The time required to complete an activity can vary
from the projected or expected time
If the activity is on the critical path, the
completion time of the project will change
This will also have an impact on ES, EF, LS, and
LF times for other activities
The exact impact depends on the relationship
between the various activities
A predecessor activity is one that must be
accomplished before the given activity can be
started
A successor activity is one that can be started
only after the given activity is finished
© 2009 Prentice-Hall, Inc. 13 – 40
Sensitivity Analysis and
Project Management
Impact of an increase (decrease) in an activity
time for a critical path activity
SUCCESSOR PARALLEL PREDECESSOR
ACTIVITY TIME ACTIVITY ACTIVITY ACTIVITY
Earliest start Increase (decrease) No change No change
Earliest finish Increase (decrease) No change No change
Latest start Increase (decrease) Increase (decrease) No change
Latest finish Increase (decrease) Increase (decrease) No change
Slack No change Increase (decrease) No change
Table 13.4
© 2009 Prentice-Hall, Inc. 13 – 41
PERT/COST
Although PERT is an excellent method of
monitoring and controlling project length, it does
not consider the very important factor of project
cost
PERT/Cost is a modification of PERT that allows a
manager to plan, schedule, monitor, and control
cost as well as time
Using PERT/Cost to plan, schedule, monitor, and
control project cost helps accomplish the sixth
and final step of PERT
© 2009 Prentice-Hall, Inc. 13 – 42
Planning and Scheduling Project Costs:
Budgeting Process
The overall approach in the budgeting
process of a project is to determine how
much is to be spent every week or month
This can be accomplished in four basic
budgeting steps
© 2009 Prentice-Hall, Inc. 13 – 43
Four Steps of the Budgeting Process
1. Identify all costs associated with each of the
activities then add these costs together to get
one estimated cost or budget for each activity
2. In large projects, activities can be combined into
larger work packages. A work package is simply
a logical collection of activities.
3. Convert the budgeted cost per activity into a
cost per time period by assuming that the cost of
completing any activity is spent at a uniform rate
over time
4. Using the ES and LS times, find out how much
money should be spent during each week or
month to finish the project by the date desired
© 2009 Prentice-Hall, Inc. 13 – 44
Budgeting for General Foundry
The Gantt chart in Figure 13.9 illustrates this
project
The horizontal bars shown when each activity will
be performed based on its ES-EF times
We determine how much will be spent on each
activity during each week and fill these amounts
into a chart in place of the bars
The following two tables show the activity costs
and budgeted cost for the General Foundry
project
© 2009 Prentice-Hall, Inc. 13 – 45
Budgeting for General Foundry
Gantt chart General Foundry project
A
B
C
Activity
D
E
F
G
H
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Week
Figure 13.9
© 2009 Prentice-Hall, Inc. 13 – 46
Budgeting for General Foundry
Activity costs for General Foundry
EARLIEST LATEST TOTAL BUDGETED
START, START, EXPECTED BUDGETED COST PER
ACTIVITY ES LS TIME, t COST ($) WEEK ($)
A 0 0 2 22,000 11,000
B 0 1 3 30,000 10,000
C 2 2 2 26,000 13,000
D 3 4 4 48,000 12,000
E 4 4 4 56,000 14,000
F 4 10 3 30,000 10,000
G 8 8 5 80,000 16,000
H 13 13 2 16,000 8,000
Total 308,000
Table 13.5
© 2009 Prentice-Hall, Inc. 13 – 47
Budgeting for General Foundry
Budgeted cost for General Foundry
WEEK
ACTIVITY 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 TOTAL
A 11 11 22
B 10 10 10 30
C 13 13 26
D 12 12 12 12 48
E 14 14 14 14 56
F 10 10 10 30
G 16 16 16 16 16 80
H 8 8 16
308
Total per week 21 21 23 25 36 36 36 14 16 16 16 16 16 8 8
Total to date 21 42 65 90 126 162 198 212 228 244 260 276 292 300 308
Table 13.6
© 2009 Prentice-Hall, Inc. 13 – 48
Budgeting for General Foundry
It is also possible to prepare a budget based on
the latest starting time
This budget will delay the expenditure of funds
until the last possible moment
The following table shows the latest start budget
for the General Foundry project
The two tables form a budget range
Any budget can be chosen between these two
values depending on when the company wants to
actually spend the money
The budget ranges are plotted in Figure 13.10
© 2009 Prentice-Hall, Inc. 13 – 49
Budgeting for General Foundry
Late start budgeted cost for General Foundry
WEEK
ACTIVITY 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 TOTAL
A 11 11 22
B 10 10 10 30
C 13 13 26
D 12 12 12 12 48
E 14 14 14 14 56
F 10 10 10 30
G 16 16 16 16 16 80
H 8 8 16
308
Total per week 11 21 23 23 26 26 26 26 16 16 26 26 26 8 8
Total to date 11 32 55 78 104 130 156 182 198 214 240 266 292 300 308
Table 13.7
© 2009 Prentice-Hall, Inc. 13 – 50
Budgeting for General Foundry
Total
Budgeted A manager can
Cost
$300,000 – choose any budget
Budget Using
that falls between
250,000 – Earliest Start the budgets
Times, ES presented in the
200,000 –
two tables
The two tables
150,000 –
form feasible
Budget Using budget ranges
Latest Start
100,000 – Times, LS
50,000 –
0–
Figure 13.10
| | | | | | | | | | | | | | |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Weeks
© 2009 Prentice-Hall, Inc. 13 – 51
Monitoring and Controlling
Project Costs
Costs are monitored and controlled to ensure the
project is progressing on schedule and that cost
overruns are kept to a minimum
The status of the entire project should be checked
periodically
The project is now in it’s 6th week of 15 weeks
Activities A,B, and C have completed at costs of
$20,000, $36,000 and $26,000 respectively
Activity D is only 10% complete at a cost of $6,000
Activity E is 20% complete at a cost of $20,000
Activity F is 20% complete with a cost of $4,000
What is the value of the work completed?
Are there any cost overruns?
© 2009 Prentice-Hall, Inc. 13 – 52
Monitoring and Controlling
Project Costs
Monitoring and controlling budgeted cost
VALUE OF
TOTAL WORK ACTIVITY
BUDGETED PERCENT OF COMPLETED ACTUAL DIFFERENCE
ACTIVITY COST ($) COMPLETION ($) COST ($) ($)
A 22,000 100 22,000 20,000 –2,000
B 30,000 100 30,000 36,000 6,000
C 26,000 100 26,000 26,000 0
D 48,000 10 4,800 6,000 1,200
E 56,000 20 11,200 20,000 8,800
F 30,000 20 6,000 4,000 –2,000
G 80,000 0 0 0 0
H 16,000 0 0 0 0
Total 100,000 112,000 12,000
Table 13.8 Overrun
© 2009 Prentice-Hall, Inc. 13 – 53
Monitoring and Controlling
Project Costs
The value of work completed, or the cost to date
for any activity, can be computed as follows
Value of work (Percentage of work complete)
=
completed x (Total activity budget)
The activity difference is also of interest
Activity difference = Actual cost
– Value of work completed
A negative activity difference is a cost underrun
and a positive activity difference is a cost overrun
© 2009 Prentice-Hall, Inc. 13 – 54
Monitoring and Controlling
Project Costs
Value completed is $100,000 while actual cost
is $112,000; cost overrun of $12,000
Using the earliest start times budget, by the
end of the 6th week we should have completed
75% of D (vs 10%), 50% of E (vs 20%) and 66.7% of
F (vs 20%) and spent $162,000 so the project is
behind schedule
Using the latest start times budget, by the end
of the 6th week we should have completed
50% of D (vs 10%), 50% of E (vs 20%) and 0% of F
(vs 20%) and spent $130,000 so the project is also
behind schedule
© 2009 Prentice-Hall, Inc. 13 – 55
Project Crashing
Projects will sometimes have deadlines
that are impossible to meet using normal
procedures
By using exceptional methods it may be
possible to finish the project in less time
than normally required
However, this usually increases the cost
of the project
Reducing a project’s completion time is
called crashing
© 2009 Prentice-Hall, Inc. 13 – 56
Project Crashing
Crashing a project starts with using the
normal time to create the critical path
The normal cost is the cost for completing
the activity using normal procedures
If the project will not meet the required
deadline, extraordinary measures must be
taken
The crash time is the shortest possible
activity time and will require additional
resources
The crash cost is the price of completing the
activity in the earlier-than-normal time
© 2009 Prentice-Hall, Inc. 13 – 57
Four Steps to Project Crashing
1. Find the normal critical path and identify
the critical activities
2. Compute the crash cost per week (or
other time period) for all activities in the
network using the formula
Crash cost – Normal cost
Crash cost/Time period =
Normal time – Crash time
© 2009 Prentice-Hall, Inc. 13 – 58
Four Steps to Project Crashing
3. Select the activity on the critical path
with the smallest crash cost per week
and crash this activity to the maximum
extent possible or to the point at which
your desired deadline has been reached
4. Check to be sure that the critical path
you were crashing is still critical. If the
critical path is still the longest path
through the network, return to step 3. If
not, find the new critical path and return
to step 2.
© 2009 Prentice-Hall, Inc. 13 – 59
General Foundry Example
General Foundry has been given 14 weeks instead of
16 weeks to install the new equipment
The critical path for the project is 15 weeks
What options do they have?
The normal and crash times and costs are shown in
Table 13.9
Crash costs are assumed to be linear and Figure 13.11
shows the crash cost for activity B
Crashing activity A will shorten the completion time to
14 but it creates a second critical path B,D,G,H
because when you recalculate the LF and LS times for
B and D they now match the EF and ES
Any further crashing must be done to both critical
paths
© 2009 Prentice-Hall, Inc. 13 – 60
General Foundry Example
Normal and crash data for General Foundry
TIME (WEEKS) COST ($)
CRASH
COST PER CRITICAL
ACTIVITY NORMAL CRASH NORMAL CRASH WEEK ($) PATH?
A 2 1 22,000 23,000 1,000 Yes
B 3 1 30,000 34,000 2,000 No
C 2 1 26,000 27,000 1,000 Yes
D 4 3 48,000 49,000 1,000 No
E 4 2 56,000 58,000 1,000 Yes
F 3 2 30,000 30,500 500 No
G 5 2 80,000 86,000 2,000 Yes
H 2 1 16,000 19,000 3,000 Yes
Table 13.9
© 2009 Prentice-Hall, Inc. 13 – 61
General Foundry - QM
© 2009 Prentice-Hall, Inc. 13 – 62
General Foundry - QM
© 2009 Prentice-Hall, Inc. 13 – 63
Revised Path After Crashing
After crashing the project by 1 week, this is the new
network
Two critIcal paths
A-C-E-G-H
B-D-G-H NODE Time ES EF LS LF
A 1 0 1 0 1
B 3 0 3 0 3
C 2 1 3 1 3
D 4 3 7 3 7
E 4 3 7 3 7
F 3 3 6 9 12
G 5 7 12 7 12
H 2 12 14 12 14
© 2009 Prentice-Hall, Inc. 13 – 64
Chapter 5
Forecasting
To accompany
Quantitative Analysis for Management, Eleventh Edition,
by Render, Stair, and Hanna
Power Point slides created by Brian Peterson
Forecasting Models
Forecasting
Techniques
Qualitative Time-Series Causal
Models Methods Methods
Delphi Moving Regression
Methods Average Analysis
Jury of Executive Exponential Multiple
Opinion Smoothing Regression
Sales Force Trend
Composite Projections
Figure 5.1
Consumer
Decomposition
Market Survey
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-2
Qualitative Models
Qualitative models incorporate judgmental
or subjective factors.
These are useful when subjective factors
are thought to be important or when
accurate quantitative data is difficult to
obtain.
Common qualitative techniques are:
Delphi method.
Jury of executive opinion.
Sales force composite.
Consumer market surveys.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-3
Qualitative Models
Delphi Method – This is an iterative group
process where (possibly geographically
dispersed) respondents provide input to decision
makers.
Jury of Executive Opinion – This method collects
opinions of a small group of high-level managers,
possibly using statistical models for analysis.
Sales Force Composite – This allows individual
salespersons estimate the sales in their region
and the data is compiled at a district or national
level.
Consumer Market Survey – Input is solicited from
customers or potential customers regarding their
purchasing plans.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-4
Time-Series Models
Time-series models attempt to predict
the future based on the past.
Common time-series models are:
Moving average.
Exponential smoothing.
Trend projections.
Decomposition.
Regression analysis is used in trend
projections and one type of
decomposition model.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-5
Scatter Diagrams
Wacker Distributors wants to forecast sales for
three different products (annual sales in the table,
in units):
TELEVISION COMPACT DISC
YEAR RADIOS
SETS PLAYERS
1 250 300 110
2 250 310 100
3 250 320 120
4 250 330 140
5 250 340 170
6 250 350 150
7 250 360 160
8 250 370 190
Table 5.1 9 250 380 200
10 250 390 190
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-6
Scatter Diagram for TVs
(a)
Sales appear to be
Annual Sales of Televisions
330 –
constant over time
250 –
Sales = 250
200 –
A good estimate of
150 – sales in year 11 is
100 – 250 televisions
50 –
| | | | | | | | | |
0 1 2 3 4 5 6 7 8 9 10
Time (Years)
Figure 5.2a
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-7
Scatter Diagram for Radios
(b)
420 – Sales appear to be
400 – increasing at a
Annual Sales of Radios
380 – constant rate of 10
360 –
radios per year
340 –
Sales = 290 + 10(Year)
320 –
300 –
A reasonable
280 –
estimate of sales in
year 11 is 400
| | | | | | | | | |
0 1 2 3 4 5 6 7 8 9 10 radios.
Time (Years)
Figure 5.2b
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-8
Scatter Diagram for CD
Players
This trend line may
(c) not be perfectly
Annual Sales of CD Players
200 –
accurate because
180 –
of variation from
160 – year to year
140 – Sales appear to be
increasing
120 –
A forecast would
100 – probably be a
| | | | | | | | | |
larger figure each
0 1 2 3 4 5 6 7 8 9 10 year
Time (Years)
Figure 5.2c
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-9
Measures of Forecast Accuracy
We compare forecasted values with actual values
to see how well one model works or to compare
models.
Forecast error = Actual value – Forecast value
One measure of accuracy is the mean absolute
deviation (MAD):
MAD
forecast error
n
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-10
Measures of Forecast Accuracy
Using a naïve forecasting model we can compute the MAD:
ACTUAL
SALES OF ABSOLUTE VALUE OF
CD FORECAST ERRORS (DEVIATION),
YEAR PLAYERS SALES (ACTUAL – FORECAST)
1 110 — —
2 100 110 |100 – 110| = 10
3 120 100 |120 – 110| = 20
4 140 120 |140 – 120| = 20
5 170 140 |170 – 140| = 30
6 150 170 |150 – 170| = 20
7 160 150 |160 – 150| = 10
8 190 160 |190 – 160| = 30
9 200 190 |200 – 190| = 10
10 190 200 |190 – 200| = 10
Table 5.2 11 — 190 —
Sum of |errors| = 160
MAD = 160/9 = 17.8
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-11
Measures of Forecast Accuracy
Using a naïve forecasting model we can compute the MAD:
ACTUAL ABSOLUTE VALUE OF
SALES OF CD FORECAST ERRORS (DEVIATION),
YEAR PLAYERS SALES (ACTUAL – FORECAST)
1 110 — —
2 100 110 |100 – 110| = 10
3 120 100 |120 – 110| = 20
MAD
4
5
forecast error 140
170
120
140
160
17.8
|140 – 120| = 20
|170 – 140| = 30
6 150 n 170 9 |150 – 170| = 20
7 160 150 |160 – 150| = 10
8 190 160 |190 – 160| = 30
9 200 190 |200 – 190| = 10
10 190 200 |190 – 200| = 10
11 — 190 —
Sum of |errors| = 160
MAD = 160/9 = 17.8
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-12
Table 5.2
Measures of Forecast Accuracy
There are other popular measures of forecast
accuracy.
The mean squared error:
MSE
( error ) 2
n
The mean absolute percent error:
error
actual
MAPE 100 %
n
And bias is the average error.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-13
Time-Series Forecasting Models
A time series is a sequence of evenly
spaced events.
Time-series forecasts predict the future
based solely on the past values of the
variable, and other variables are
ignored.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-14
Components of a Time-Series
A time series typically has four components:
1. Trend (T) is the gradual upward or downward
movement of the data over time.
2. Seasonality (S) is a pattern of demand
fluctuations above or below the trend line that
repeats at regular intervals.
3. Cycles (C) are patterns in annual data that
occur every several years.
4. Random variations (R) are “blips” in the data
caused by chance or unusual situations, and
follow no discernible pattern.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-15
Decomposition of a Time-Series
Product Demand Charted over 4 Years, with Trend
and Seasonality Indicated
Demand for Product or Service
Trend
Component
Seasonal Peaks
Actual
Demand
Line
Average Demand
over 4 Years
| | | |
Figure 5.3 Year Year Year Year
1 2 3 4
Time
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-16
Decomposition of a Time-Series
There are two general forms of time-series
models:
The multiplicative model:
Demand = T x S x C x R
The additive model:
Demand = T + S + C + R
Models may be combinations of these two
forms.
Forecasters often assume errors are
normally distributed with a mean of zero.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-17
Moving Averages
Moving averages can be used when
demand is relatively steady over time.
The next forecast is the average of the
most recent n data values from the time
series.
This methods tends to smooth out short-
term irregularities in the data series.
Sum of demands in previous n periods
Moving average forecast
n
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-18
Moving Averages
Mathematically:
Yt Yt 1 ... Yt n1
Ft 1
n
Where:
Ft 1 = forecast for time period t + 1
Yt = actual value in time period t
n = number of periods to average
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-19
Wallace Garden Supply
Wallace Garden Supply wants to
forecast demand for its Storage Shed.
They have collected data for the past
year.
They are using a three-month moving
average to forecast demand (n = 3).
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-20
Wallace Garden Supply
MONTH ACTUAL SHED SALES THREE-MONTH MOVING AVERAGE
January 10
February 12
March 13
April 16 (10 + 12 + 13)/3 = 11.67
May 19 (12 + 13 + 16)/3 = 13.67
June 23 (13 + 16 + 19)/3 = 16.00
July 26 (16 + 19 + 23)/3 = 19.33
August 30 (19 + 23 + 26)/3 = 22.67
September 28 (23 + 26 + 30)/3 = 26.33
October 18 (26 + 30 + 28)/3 = 28.00
November 16 (30 + 28 + 18)/3 = 25.33
December 14 (28 + 18 + 16)/3 = 20.67
January — (18 + 16 + 14)/3 = 16.00
Table 5.3
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-21
Weighted Moving Averages
Weighted moving averages use weights to put
more emphasis on previous periods.
This is often used when a trend or other pattern is
emerging.
Ft 1
( Weight in period i )( Actual value in period)
( Weights)
Mathematically:
w1Yt w2Yt 1 ... wnYt n1
Ft 1
w1 w2 ... wn
where
wi = weight for the ith observation
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-22
Wallace Garden Supply
Wallace Garden Supply decides to try a
weighted moving average model to forecast
demand for its Storage Shed.
They decide on the following weighting
scheme:
WEIGHTS APPLIED PERIOD
3 Last month
2 Two months ago
1 Three months ago
3 x Sales last month + 2 x Sales two months ago + 1 X Sales three months ago
6
Sum of the weights
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-23
Wallace Garden Supply
THREE-MONTH WEIGHTED
MONTH ACTUAL SHED SALES MOVING AVERAGE
January 10
February 12
March 13
April 16 [(3 X 13) + (2 X 12) + (10)]/6 = 12.17
May 19 [(3 X 16) + (2 X 13) + (12)]/6 = 14.33
June 23 [(3 X 19) + (2 X 16) + (13)]/6 = 17.00
July 26 [(3 X 23) + (2 X 19) + (16)]/6 = 20.50
August 30 [(3 X 26) + (2 X 23) + (19)]/6 = 23.83
September 28 [(3 X 30) + (2 X 26) + (23)]/6 = 27.50
October 18 [(3 X 28) + (2 X 30) + (26)]/6 = 28.33
November 16 [(3 X 18) + (2 X 28) + (30)]/6 = 23.33
December 14 [(3 X 16) + (2 X 18) + (28)]/6 = 18.67
January — [(3 X 14) + (2 X 16) + (18)]/6 = 15.33
Table 5.4
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-24
Exponential Smoothing
Exponential smoothing is a type of moving
average that is easy to use and requires little
record keeping of data.
New forecast = Last period’s forecast
+ (Last period’s actual demand
– Last period’s forecast)
Here is a weight (or smoothing constant) in
which 0≤≤1.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-25
Exponential Smoothing
Mathematically:
Ft 1 Ft (Yt Ft )
Where:
Ft+1 = new forecast (for time period t + 1)
Ft = pervious forecast (for time period t)
= smoothing constant (0 ≤ ≤ 1)
Yt = pervious period’s actual demand
The idea is simple – the new estimate is the old
estimate plus some fraction of the error in the
last period.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-26
Exponential Smoothing Example
In January, February’s demand for a certain
car model was predicted to be 142.
Actual February demand was 153 autos
Using a smoothing constant of = 0.20, what
is the forecast for March?
New forecast (for March demand) = 142 + 0.2(153 – 142)
= 144.2 or 144 autos
If actual demand in March was 136 autos, the
April forecast would be:
New forecast (for April demand) = 144.2 + 0.2(136 – 144.2)
= 142.6 or 143 autos
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-27
Selecting the Smoothing Constant
Selecting the appropriate value for is key
to obtaining a good forecast.
The objective is always to generate an
accurate forecast.
The general approach is to develop trial
forecasts with different values of and
select the that results in the lowest MAD.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-28
Exponential Smoothing
Port of Baltimore Exponential Smoothing Forecast
for =0.1 and =0.5.
ACTUAL
TONNAGE FORECAST FORECAST
QUARTER UNLOADED USING =0.10 USING =0.50
1 180 175 175
2 168 175.5 = 175.00 + 0.10(180 – 175) 177.5
3 159 174.75 = 175.50 + 0.10(168 – 175.50) 172.75
4 175 173.18 = 174.75 + 0.10(159 – 174.75) 165.88
5 190 173.36 = 173.18 + 0.10(175 – 173.18) 170.44
6 205 175.02 = 173.36 + 0.10(190 – 173.36) 180.22
7 180 178.02 = 175.02 + 0.10(205 – 175.02) 192.61
8 182 178.22 = 178.02 + 0.10(180 – 178.02) 186.30
9 ? 178.60 = 178.22 + 0.10(182 – 178.22) 184.15
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-29
Table 5.5
Exponential Smoothing
Absolute Deviations and MADs for the Port of
Baltimore Example
ACTUAL ABSOLUTE ABSOLUTE
TONNAGE FORECAST DEVIATIONS FORECAST DEVIATIONS
QUARTER UNLOADED WITH = 0.10 FOR = 0.10 WITH = 0.50 FOR = 0.50
1 180 175 5….. 175 5….
2 168 175.5 7.5.. 177.5 9.5..
3 159 174.75 15.75 172.75 13.75
4 175 173.18 1.82 165.88 9.12
5 190 173.36 16.64 170.44 19.56
6 205 175.02 29.98 180.22 24.78
7 180 178.02 1.98 192.61 12.61
8 182 178.22 3.78 186.30 4.3..
Sum of absolute deviations 82.45 98.63
Σ|deviations|
MAD = = 10.31 MAD = 12.33
n
Table 5.6
Best choice
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-30
Trend Projections
Trend projection fits a trend line to a
series of historical data points.
The line is projected into the future for
medium- to long-range forecasts.
Several trend equations can be
developed based on exponential or
quadratic models.
The simplest is a linear model developed
using regression analysis.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-31
Trend Projection
The mathematical form is
Yˆ b0 b1 X
Where
Ŷ = predicted value
b0 = intercept
b1 = slope of the line
X = time period (i.e., X = 1, 2, 3, …, n)
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-32
Midwestern Manufacturing
Midwest Manufacturing has a demand for
electrical generators from 2004 – 2010 as given
in the table below.
YEAR ELECTRICAL GENERATORS SOLD
2004 74
2005 79
2006 80
2007 90
2008 105
2009 142
Table 5.7
2010 122
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-33
Midwestern Manufacturing
Company Example
The forecast equation is
Yˆ 56.71 10.54 X
To project demand for 2011, we use the coding
system to define X = 8
(sales in 2011) = 56.71 + 10.54(8)
= 141.03, or 141 generators
Likewise for X=9
(sales in 2012) = 56.71 + 10.54(9)
= 151.57, or 152 generators
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-34
Midwestern Manufacturing
Electrical Generators and the Computed Trend Line
Figure 5.4
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-35
Trend projections
Observing a time series for the GDP of
Jordan from 2001 to 2011, annual data
20
18
16
14
12
JD bl
10
8
6
4
2
0
Souce:
0
2000 IMF
2 World
4 Development
6
2005 8 Database
10
2010 12
year
36
Trend projections
Observing a time series for the GDP of
Jordan from 2001 to 2011, annual data
20
18 y = 1.4431x + 3.2106
16 R² = 0.9355
14
12
JD bl
10
8
6
4
2
0
Souce:
0
2000 IMF
2 World
4 Development
6
2005 8 Database
10
2010 12
year
37
Trend projections
Different trend models can be
compared
E.g.
20 linear versus quadratic
18
16
y = 4.9615x0.4884
14
R² = 0.8158
12
y = 1.4431x + 3.2106
JD bl
10
R² = 0.9355
8
0
0 2 4 6 8 10 12
2000 2005
year 2010
38
Trend projections
Different trend models can be
compared
20
18 E.g. linear versus
y =exponential
5.092e0.1294x
16 R² = 0.977
14
12 y = 1.4431x + 3.2106
JD bl
10 R² = 0.9355
8
0
0 2 4 6 8 10 12
2000 2005
year 2010
39
Seasonal Variations
Recurring variations over time may
indicate the need for seasonal
adjustments in the trend line.
A seasonal index indicates how a
particular season compares with an
average season.
When no trend is present, the seasonal
index can be found by dividing the
average value for a particular season by
the average of all the data.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-40
Eichler Supplies
Eichler Supplies sells telephone
answering machines.
Sales data for the past two years has
been collected for one particular model.
The firm wants to create a forecast that
includes seasonality.
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-41
Eichler Supplies Answering Machine
Sales and Seasonal Indices
AVERAGE
SALES DEMAND
AVERAGE TWO- MONTHLY SEASONAL
MONTH YEAR 1 YEAR 2 YEAR DEMAND DEMAND INDEX
January 80 100 90 94 0.957
February 85 75 80 94 0.851
March 80 90 85 94 0.904
April 110 90 100 94 1.064
May 115 131 123 94 1.309
June 120 110 115 94 1.223
July 100 110 105 94 1.117
August 110 90 100 94 1.064
September 85 95 90 94 0.957
October 75 85 80 94 0.851
November 85 75 80 94 0.851
December 80 80 80 94 0.851
Total average demand = 1,128
1,128 Average two-year demand
Average monthly demand = = 94 Seasonal index = Average monthly demand
12 months
Table 5.9
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-42
Seasonal Variations
The calculations for the seasonal indices are
1,200 1,200
Jan. 0.957 96 July 1.117 112
12 12
1,200 1,200
Feb. 0.851 85 Aug. 1.064 106
12 12
1,200 1,200
Mar. 0.904 90 Sept. 0.957 96
12 12
1,200 1,200
Apr. 1.064 106 Oct. 0.851 85
12 12
1,200 1,200
May 1.309 131 Nov. 0.851 85
12 12
1,200 1,200
June 1.223 122 Dec. 0.851 85
12 12
Copyright ©2012 Pearson Education, Inc. publishing as Prentice Hall 5-43