Module OR Final
Module OR Final
CHAPTER ONE
Objectives:
1.1. Introduction
Many people still remain in the bondage of self-incurred tutelage. Tutelage is a person's
inability to make his/her own decisions. Self-incurred is this tutelage when its cause lies
not in lack of reason but in lack of resolution and courage to use it without wishing to
have been told what to do by something or somebody else.
The difficulty in life is the choice. Good decision-making brings about a better life. A bad
decision may force you to make another one. A good decision is never an accident; it is
always the result of high intention, sincere effort, intelligent direction and skillful
execution; it represents the wise choice of many alternatives. One must appreciate the
difference between a decision and an objective. A good decision is the process of
optimally achieving a given objective.
When decision making is too complex or the interests at stake are too important, quite
often we do not know or are not sure what to decide. In many instances, we resort to
informal decision support techniques such as tossing a coin, asking an oracle, visiting an
astrologer, etc. However formal decision support from an expert has many advantages.
Business Science focuses on the formal model-driven decision support techniques such as
mathematical programs for optimization, and decision tree analysis for risky decisions.
Such techniques are now part of our everyday life. For example, when a bank must
decide whether a given client will obtain credit or not, a technique, called credit scoring,
is often used.
Rational decisions are often made unwillingly, perhaps unconsciously. We may start the
process of consideration. It is best to learn the decision-making process for complex,
important and critical decisions. Critical decisions are those that cannot and must not be
wrong.
The aim of this course is to make you a better decision maker by learning the decision-
making process. Decision-making is a complicated process. This complication arises
from the fact that your present goal (including wants, resources, and abilities) dictates
your choices; however, your choices will change your goals. This influential-cycle keeps
the decision-maker busy all the time. Selecting your goals and your criteria for success is
a dynamic process and changes over time.
On a daily basis a manager has to make many decisions. Some of these decisions are
routine and inconsequential, while others have drastic impacts on the operations of the
firm for which he/she works. Some of these decisions could involve large sums of money
being gained or lost, or could involve whether or not the firm accomplishes its mission
and its goals.
In our increasingly complex world, the tasks of decision-makers are becoming more
challenging with each passing day. The decision-maker (i.e., the responsible manager)
must respond quickly to events that seem to take place at an ever-increasing pace. In
addition, a decision-maker must incorporate a sometimes-bewildering array of choices
and consequences into his or her decision. Routine decisions are often made quickly,
perhaps unconsciously without the need for a detailed process of consideration. However,
for complex, critical or important managerial decisions it is necessary to take time to
decide systematically.
Making strategic decisions requires that one takes a structured approach following a
formal decision making process. Otherwise, it will be difficult to be sure that one has
considered all the key aspects of the decision.
It is generally agreed that OR come into existence as a discipline during world war II
when there was a critical need to manage scarce resources. However, a particular model
and technique of OR can be traced back to much earlier times. The scientific approach to
management can be traced back to the era of Industrial Revolution and even periods
before that. But operations research, as it exists today, was born during the Second World
War when the British military management called upon a group of scientists to examine
the strategies and tactics of various military operations with the intention of efficient
allocation of scarce resources to the war effort.
The war involved strategic and tactical problems which were greatly complicated to
expect adequate solution from individuals or specialists in a single discipline. Therefore
groups of individuals who collectively were considered specialists in mathematics,
economics, statistics and probability theory, engineering, behavioral and physical science
were formed as a special units within the armed forces to deal with strategic and tactical
problems of various military operations.
Problems of this latter kind, often called systems problems, were studied rigorously for
the first time by scientists, not in industry but in the management of operations in the
1939 - 45 war.
One such team (the Radar Operational Research team - which gave OR its name) was
responsible for implementing the installation and operation of radar after the technical
development work was complete.
Later on another team examined the relative ineffectiveness of the Allied Forces at
destroying the German U-boats, which were sinking the food convoy ships supplying
Britain.
The OR team played an important part in reducing shipping losses and in the ultimate
defeat of the U-boats. With no precedent in the application of Science to the management
function, the scientists were still able to use their scientific approach-collecting
information and developing hypothesis in order to come up with practical plans for the
improvement of these wartime operations.
After the end of the war, the same approach was used with great success in the industrial
and commercial field, developing most rapidly in Britain and the USA, followed by
Europe and the rest of the world. This has provided the basis for the evolution of
Operational Research (OR) as a separate and independent subject of Science. It is not
based on any single academic engineering, social science, economics, statistics and
computing but is none of these.
There are three different factors behind the rapid development in the use and
development of operations research approach.
The economic and industrial boom after World War II resulted in continuous
mechanization, automation, decentralization of operations and decisions of
management functions. This industrialization also resulted in complex managerial
problems and therefore application of operations research to managerial decision
making become popular.
Many operations researchers continued their research after war. A key person in
the post war development of OR was George B. Dantzing. In 1947 he developed
the concept of linear programming and its solution by a method known as
Simplex method.
Analytic power was made available by high-speed computers. The development
of electronic digital computers, with their ability to perform arithmetic
calculations thousands or even millions of times faster than a human being can,
was a tremendous boom to OR.
1.3. Meaning of operations research
Definition
In the USA, where it is called Operations Research, the OR Society of America says
more briefly;
"OR is concerned with scientifically deciding how to best design and operate man-
machine systems, usually under conditions requiring the allocation of scarce
resources"
From the concept and definition given above, Operations Research is:
As the term implies, OR involves research on military operations. The operation research
approach is particularly useful in balancing conflicting objectives where there are many
alternative courses of action available to the decision maker. In a theoretical sense, the
optimum decision (global optimum) must be one that is best for the organization as a
whole.
Operations research attempts to resolve the conflicts of interest among various sections of
the organization and seeks the optimal solution which may not be acceptable to one
department but is in the interest of the organization as a whole.
Existence of a problem
Intension to solve a problem
Application of system concept and system analysis to the problem
Scientific approach to the problem where research methods are used
Formation of group consists of different specialists
Assist management to make scientific decision
Solution must meet the objective within the constraints
OR is concerned with providing the decision maker with decision aids (or rules) derived
from:
Interdisciplinary Approach
Methodological approach
Wholistic Approach
While arriving at a decision, an operation research team examines the relative importance
of all conflicting and multiple objectives, and the validity of claims of various
departments of the organization from the perspective of the whole organization.
Objectivistic Approach
An operation research approach seeks to obtain an optimal solution to the problem under
analysis. For this a measure of desirability (effectiveness) is defined, based on the
objective(s) of the organization. A measure of desirability so defined is then used to
compare alternative courses of action with respect to their outcomes.
Observation
Problem definition
Feed
Model construction
Back Management science technique
Solution
Information
Implementation
Dear learners, here we will discuss models and modeling in operations research.
The approximation or abstraction, maintaining only the essential elements of the system,
which may be constructed with various forms by establishing relationships among
specified variables and parameters of the system, is called a model.
the innumerable and changing characteristics of the real life problems to be represented.
Instead, they are limited approximation of reality.
A model is constructed to analyze and understand the given system for the purpose of
improving its performance. The reliability of the solution obtained from a model depends
on the validity of the model in representing the system under study.
The key to model building lies in abstracting only the relevant variables that affect the
criteria of the measures of performance of the given system and expressing the
relationship in a situation form. A model should be as simple as possible so as to give the
desired result.
There are many ways to classify models. Classification schemes can also provide a useful
frame of reference for modelers. There are five classification schemes for models.
a) Physical models
The models provide a physical appearance of the real object under study either reduced in
size or scaled up. Physical models are useful only in design problems because they are
easy to observe, build and describe.
The chief merit of an iconic model is that it is concrete and specific. It resembles visually
the thing it represents and therefore, there are likely to be fewer problems in translating
any “findings” from the model into the real-life situation. However, the disadvantage of
such models is that they often do not lead themselves to manipulation for experimental
purpose.
Use symbols (letters, numbers) and functions to represent variables and their
relationships to describe the properties of the system. These models are also used to
represent relationships which can be represented in a physical form.
ii. Mathematical models: - involve the use of mathematical symbols, letters, numbers
and mathematical operations (+, –, ÷, ×) to represent relationships among
variables of the system to describe its properties or behaviors.
Example: P = 2L + 2W2
Subject ¿:
x 1+ x 2<60
x 1+2 x 2< 45
x 1<15
x 1 , x 2> 0
b) Predictive models: - these models indicate “if these occurs then what will follow”.
They relate dependent and independent variables and permit trying out “what if
questions’. Used to predict the outcomes due to a given set of alternatives for the
problem.
c) Normative (optimization) models: - provide the best or optimal solution to problems
subject to certain limitations on the use of resources. Are also called prescriptive
models because they provide what the decision maker out to do.
a. Static models: - static models represent a system at some specified time and do
not account for changes over time.
b. Dynamic models; - time is considered as one of the variables and allows the
impact of changes due to change in time. Thus the sequence of interrelated
decisions over a period of time is made to select the optimal course of action to
optimize the given objective.
a) Deterministic models
If all the parameters, constants and functional relationships are assumed to be known with
certainty when the decision is made, the model is said to be deterministic. Thus, in such
case, the outcome associated with a particular course of action is known. That is for a
specific set of input values, there is a uniquely determined output which represents the
solution of the model under considerations of certainty. The results of the models assume
single value.
Models in which at least one parameter or decision variable is a random variable are
called probabilistic (or stochastic) models. Since at least one decision variable is random,
therefore, an independent variable which is the function of dependent variable(s) will also
be random. This means consequence or payoff due to certain changes in the independent
variable cannot be predicted with certainty. However, it is possible to predict a pattern of
values of both the variables by their probability distribution.
Example: insurance against risk of fire, accidents, sickness, where the pattern of events is
studied in the form of probability distribution.
a. Heuristic models: These models employ some set of rules which, though
perhaps not optimal, do facilitate solution of problems when applied in a
consistent manner.
b. Analytical models: it has specific analytical structure and thus can be solved
by known analytical or mathematical techniques. Any optimization model
Chapter Summary
In our increasing complex world, the tasks of decision-makers are becoming more
challenging with each passing day. The decision-maker is expected to respond quickly to
events that seem to take place at an ever increasing pace.
Operations research attempts to resolve the conflicts of interest among various sections of
the organization and seeks the optimal solution which may not be acceptable to one
department but is in the interest of the organization as a whole.
It is generally agreed that OR come into existence as a discipline during the second
world war when there was a critical need to manage scarce resources.
The operation research problem solving approach involves five major steps: observation,
definition of the problem, model construction and implementation of solution results.
The approximation or abstraction, maintaining only the essential elements of the system,
which may be constructed with various forms by establishing relationships among
specified variables and parameters of the system, is called a model.
A model is constructed to analyze and understand the given system for the purpose of
improving its performance. The reliability of the solution obtained from a model depends
on the validity of the model in representing the system under study. There are different
classifications of models which are constructed to represent and resolve real world
business problem.
Exercises
a) Controllable c) Parameters
b) Uncontrollable d) None of the above
3. A model is
a) Multi-disciplinary
b) Scientific
c) Intuitive
d) All of the above
CHAPTER TWO
Linear Programming
Learning objectives:
On successful completion of this chapter, you should be able to:
all the courses of options available to an organization are known & the objective of the
firm along with its constraints are quantified. That course of action is chosen out of all
possible alternatives which yield the optimal results. Linear Programming can also be
used as a verification and checking mechanism to ascertain the accuracy and
the reliability of the decisions which are taken solely on the basis of manager's
experience- without the aid of a mathematical model.
Some of the definitions of Linear Programming are as follows:
"Linear Programming is a method of planning and operation involved in the construction
of a model of a real-life situation having the following elements:
Even though Linear Programming has wide & diverse' applications, yet all LP problems
have the following properties in common:
The objective is always the same (i.e.; profit Maximization or cost
minimization).
Presence of constraints which limit the extent to which the objective can be
achieved.
Availability of alternatives i.e.; different courses of action to choose from,
The objectives and constraints can be expressed in the form of linear relation.
Regardless of the size or complexity, all LP problems take the same form i.e.
allocating scarce resources among various competing alternatives. Irrespective of the
manner in which one defines Linear Programming, a problem must have certain
basic characteristics before this technique can be utilized to find the optimal values.
6. Linearity Criterion. The relationship among the various decision variables must be
directly proportional. Both the objective and the constraint, must be expressed in
terms of linear equations or inequalities. For example. If one of the factor inputs
(resources like material, labor, plant capacity etc.) increases, then it should result
in a proportionate manner in the final output. These linear equations can graphically
be presented as a straight line.
7. Additively. It is assumed that the total profitability and the total amount of each
resource utilized would be exactly equal to the sum of the respective individual amounts.
Thus the function or the activities must be additive - and the interaction among
the activities of the resources does not exist.
8. Mutually Exclusive Criterion. All decision parameters and the variables are assumed to
be mutually exclusive. In other words, the occurrence of anyone variable rules out the
simultaneous occurrence of other such variables.
9. Divisibility. Variables may be assigned fractional values. i.e.; they need not
necessarily always be in whole number. If a fraction of a product cannot be produced, an
integer programming problem exists. Thus, the continuous values of the decision
variables and resources must be permissible in obtaining an optimal solution.
10. Certainty. It is assumed that conditions of certainty exist i.e.; all the relevant
parameters or coefficients in the Linear Programming model are ful1y and
completely known and that they don't change during the period. However, such an
assumption may not hold good at all times.
11. Finiteness. LP assumes the presence of a finite number of activities and constraints
without which it is not possible to obtain the best or the optimal solution.
Diagrammatically,
Resource
constraints
Objective Constraints
s
Non-negativity
Constraints
Optimization
Maximize Minimization
2. Evaluation of All Possible Alternatives. Most of the problems faced by the present
organizations are highly complicated - which cannot be solved by the traditional
approach to decision making. The technique of Linear Programming ensures that all
possible solutions are generated - out of which the optimal solution can be selected.
those situations where a given problem can clearly be represented in the form of
linear relationship between different decision variables. Hence it is based on the
implicit assumption that the objective as well as all the constraints or the limiting factors
can be stated in term of linear expressions - which may not always hold good in
real life situations. In practical business problems, many objective function & constraints
cannot be expressed linearly. Most of the business problems can be expressed quite easily
in the form of a quadratic equation (having a power 2) rather than in the terms of
linear equation.
3. No Scope for Fractional Value Solutions. There is absolutely no certainty that the
solution to a LP problem can always be quantified as an integer quite often, Linear
Programming may give fractional-varied answers, which are rounded off to the next
integer. Hence, the solution wouldn’t be the optimal one. For example, in finding out the
number of men and machines required to perform a particular job, a fractional Larson-
integer solution would be meaningless.
In the last few decades since 1960s, no other mathematical tool or technique has had such
a profound impact on the management's decision making criterion as Linear
Programming well and truly it is one of the most important decision making tools of the
last century which has transformed the way in which decisions are made and businesses
are conducted.
Some of the important application areas of Linear Programming are the following:
assembly process in the optimal (best possible) sequence so that the total elapsed time
could be minimized.
7. Mixing or Blending. Such problems arise when the same product can be
produced with the help of a different variety of available raw-materials each
having a fixed composition and cost. Here the objective is to determine the
minimum cost blend or mix (i.e. the cost minimizations) and the various
constraints that operate are the availability of raw materials and restrictions on
some of the product constituents.
10. Profit Planning & Contract. Linear Programming is also quite useful in profit
planning and control. The objective is to maximize the profit margin from
investment in the plant facilities and machinery.
11. Traveling Salesmen Problem. Traveling salesmen problem refers to the problem
of a salesman to find the shortest route originating from a particular city, visiting
each of the specified cities and then returning back to the originating point of
departure. The restriction being that no city must be visited more than once during
a particular tour. Such types of problems can quite easily be solved with the help
of Linear Programming.
12. Media Selection/Evaluation. Media selection means the selection of the optimal
media-mix so as to maximize the effective exposure. The various constraints in
this case are: Budget limitation, different rates for different media (i.e.; print
media, electronic media like radio and T.V. etc.) and the minimum number of
repeated advertisements in the various media. The use of Linear Programming
facilities like the decision making process.
13. Staffing. Staffing or the man-power costs are substantial for a typical
organization which make its products or services very costly. Linear
Programming techniques help in allocating the optimum employees (man-power
or man-hours) to the job at hand. The overall objective is to minimize the total
man-power or overtime costs.
14. Job Analysis. Linear Programming is frequently used for evaluation of jobs in an
organization and also for matching the right job with the right worker.
Linear programming is one of the most useful techniques for effective decision making. It
is an optimization approach with an emphasis on providing the optimal solution for
resource allocation. How best to allocate the scarce organizational or national resources
among different competing and conflicting needs (or uses) forms the core of its working.
The scope for application of linear programming is very wide and it occupies a central
place in many diversified decisional problems. The effective use and application of linear
programming requires the formulation of a realistic model which represents accurately
the objectives of the decision making subject to the constraints in which it is required to
be made.
The basic steps in formulating a linear programming model are presented as follows:
Step I Identification of the decision variables. The decision variables (parameters)
having a bearing on the decision at hand shall first be identified, and then expressed or
determined in the form of linear algebraic functions or in equations.
a. Express each constraint in words. For this first see whether the constraint is of the
form ≥ (at least, as large as), or of the form ≤ (no larger than) or = (exactly equal)
b. Then express the objective function verbally.
c. Step (a) and (b) should then allow you to verbally identify the decision variables.
Step II Identification of the constraints. All the constraints in the given problem which
restrict the operation of a firm at a given point of time must be identified in this stage.
Further these constraints should be broken down as linear functions in terms of the pre-
defined decision variables.
Step III Formulate the objective function. In the last stage, the objective which is
required to be optimized (i.e., maximized or memorized) must be dearly identified and
expressed in terms of the pre-defined decision variables.
Example 1
3F furniture Ltd. manufactures two products, tables & chair. Both the products have to be
processed through two machines Ml & M2 the total machine-hours available are: 200
hours ofM1 and 400 hours of M2 respectively. Time in hours required for producing a
chair and a table on both the machines is as follows:
Time in Hours:
M1 7 4
M2 5 5
Profit from the Sale of table is Birr 40 and that of a chair is Birr 30.
Required: formulate the LP Problem (LPP)?
Solution:
Step I Identification of the decision variables. Let x1 = Number of tables produced and
∴ 7 x1 + 4 x 2 ≤ 200
(Since it takes 7 hours to produce a table & 4 hours to produce a chair on machine M1)
∴ 5 x 1 +5 x 2 ≤ 400
(Since it takes 5 hours to produce both a table & a chair on machine M2)
Step III. The objective function for maximizing the profit is given by
Maximize Z=40 x 1 +30 x 2 (objecti ve function)
(Since profit per unit from a table and a chair is Birr 40 & Birr. 30 respectively).
Presenting the problem as LPP, the given problem can now be formulated as a linear
programming model as follows:
Maximize Z ¿ 40 x 1 +30 x 2
(Since if x 1∧x 2 < 0 it means that negative quantities of products are being manufactured
which has no meaning).
Example 2
Alpha Limited produces & sells 2 different products under the brand name black &
white. The profit per unit on these products is Birr 50 & Birr 40 respectively. Both black
& white employ the same manufacturing process which has a fixed total capacity of
50,000 man-hours. As per the estimates of the marketing research department of Alpha
Limited, there is a market demand for maximum 8,000 units of Black & 10,000 units of
white. Subject to the overall demand, the products can be sold in any possible
combination. If it takes 3 hours to produce one unit of black & 2 hours to produce one
unit of white. Formulate the about as a linear programming model?
Solution:
Step I Identification of the decision variables. Let x1& x2 denote the number of units
(Since it takes 3 hours to produce 1 unit of x 1 & 2 hours to produce 1 unit of x 2 & the
total available man – hours are 50,000).
b) Marketing constraints:
x 1 ≤ 8,000
(Since maximum 8,000 units of x1 can be sold)
x 2 ≤ 10,000
(Since maximum 10,000 units of x2 can be sold).
Step III. The objective function for maximizing the profit is given by
Maximize Z=50 x1 + 40 x 2 ( objective function)
Presenting the problem as LPP Now, the given problem can be written as a linear
programming model in the follows:
x 1 ≤ 8,000
x 2 ≤ 10,000
Further; x 1∧x 2 ≥ 0
(Since if x1, x2< 0, it means that negative Quantities of products are being manufactured
I1 I2 O1 O2
X 2 6 3 7
Y 4 8 5 9
Maximum 500 units of Input I1 and 300 units of I2 are available to ABC Chemicals in
the local market. The forecasted demand for outputs O I & O2 are at least 5,000 units &
7,000 units respectively. The respective profits from process X & Y are Birr 1,000 & Birr
2,000 per production run.
Required: formulate the above problem as a linear programming model.
There are two types of finding a solution for Linear programming problems.
Dear learners, during the previous section, you have seen how to formulate a given
problem as a Linear Programming model. And you are also familiar with the method of
solving LP models by using graphics method (see module of mathematics for
management mgt221). So in this module we try to remind you some points of graphic
solution method.
The next step, after the formulation of LP model is to devise effective methods to solve
the model and ascertain the optimal solution.
Dear learners, we start with the graphical method of solving LPP and now let we move
on to simplex algorithm for solving the Linear Programming model.
The graphic solution procedure is one of the methods of solving two variable linear
programming problems.
Step I Defining the problem. Formulate the problem mathematically. Express it in terms
of several mathematical constraints & an objective function. The objective function
relates to the optimization aspect, i.e. Maximization or minimization Criterion.
Step II Plot the constraints Graphically. Each inequality in the constraint equation has to
be treated as an equation. An arbitrary value is assigned to one variable & the value of the
other variable is obtained by solving the equation. In the similar manner, a different
arbitrary value is again assigned to the variable & the corresponding value of other
variable is easily obtained. These 2 sets of values are now plotted on a graph and
connected by a straight line. The same procedure has to be repeated for all the
constraints. Hence, the total straight lines would be equal to the total no of equations,
each straight line representing one constraint equation.
Step III Locate the solution space. Solution space or the feasible region is the graphical
area which satisfies all the constraints at the same time. Such a solution point (x, y)
always occurs at the corner Points of the feasible Region the feasible region is determined
as follows:
a) For "greater than" & "greater than or equal to" constraints, the feasible region or the
solution space is the area that lays above the constraint lines.
b) For" Less Then" &" Less than or equal to" constraint, the feasible region or the
solution space is the area that lays below the constraint lines.
Step IV Selecting the graphic solution technique. Select the appropriate graphic technique
to be used for generating the solution. There are two graphic techniques to find solution;
I. Since the solution point (x, y) always occurs at the corner point of the feasible or
solution space, identify each of the extreme points or corner points of the feasible
region by the method of simultaneous equations.
II. By putting the value of the corner point's co-ordinates [e.g. (2, 3)] into the objective
function, calculate the profit (or the cost) at each of the corner points.
III. In a Maximize problem, the optimal solution occurs at that corner point which gives
the highest profit. In a minimization problem, the optimal solution occurs at that
corner point which gives the lowest profit.
Dear learners, let us now turn our attention to the important theorems which are used in
solving a linear programming problem.
Note. Convex set is a polygon "Convex" implies that if any two points of the polygon are
selected arbitrarily then straight line segment joining these two points lays completely
within the polygon. The extreme points of the convex set are the basic solution to the
linear programming problem.
Important Terms
Some of the important terms commonly used in linear programming are disclosed as
follows:
Example 1
XYZ Ltd. Co. Wishes to purchase a maximum of 3600 units of two types of product, A &
B are available in the market. Product A occupies a space of 3 cubic feet & cost Birr 9
whereas B occupies a space of 1 cubic feet & cost Birr 13 per unit. The budgetary
constraints of the company do not allow spending more than Birr 39,000. The total
availability of space in the company's god own is 6000 cubic feet. Profit margin of both
the product A & B is Birr 3 & Birr 4 respectively. Formulate the above problem as a
linear programming model and solve it by using graphical method. You are required to
ascertain the best possible combination of purchase of A & B so that the total profits are
maximized.
Solution
The 1st constraint is: x 1+ x2 ≤3600 c h ang e this in to x 1+ x2 =3600 and find the X and Y
intercept as follows.
Put x 1=0 ⇒ x 2=3600 ∴ t h e point is(0 , 3600)
Put x 2=0 ⇒ x 1=3600 ∴ t h e point is(3600 , 0)
Draw the graph with x1 on x-axis & x2 on y-axis as shown in the following figure.
nd
The 2 constraint is: 3 x 1+ x 2 ≤6000 change this in to 3 x 1+ x 2=6000 and find X and Y intercept.
Draw the graph with x1 on x-axis & x2 on y-axis as shown in the following figure.
rd
The 3 constraint is: 9 x 1+ 13 x 2 ≤39,000 c h ang e this in to 9 x 1+ 13 x 2=39,000 and find
the X and Y intercept as follows.
Put x 1=0 ⇒ x 2=3000 ∴ t h e point is(0 , 3000)
Put x 2=0 ⇒ x 1=4333.3 ∴ t h e point is(4333.3 , 0)
Draw the graph with x1 on x-axis & x2 on y-axis as shown in the following figure.
Step II Determine the set of the points which satisfy the constraint:
x 1+ x2 ≤3600
This can easily be done by verifying whether the origin (0, 0) satisfies the constraint.
Here, 0+ 0<3600 hence the solution is true; all the points below the line will satisfy the
constraint.
Like it’s in the above iteration, determine the set of points which satisfy the constraint
3 x 1+ x 2 ≤6000 . At origin;
0+ 0<6000. Hence, all the points below the line will satisfy the constraint.
Like it’s in the above iteration, determine the set of points which satisfy the constraint.
9 x 1+ 13 x 2 ≤39,000. At origin;
0+ 0<39,000. Hence, all the points below the line will satisfy the constraint.
The intersection of the above graphic denotes the feasible region for the given problem.
a. For ≥ constraint the feasible region will be the area, which lays above the constraint lines,
and for ≤ constraints, it will lays below the constraint lines. This would be useful in
identifying the feasible region.
b. According to a theorem on linear programming, an optimal solution to a problem (if it
exists) is found at a corner point of the solution space.
X1X1
In order to get the value of point B apply simultaneous equation by taking the two
intersection lines. Solve the point.
3 x 1+ x 2 ≤6000 (Storage area constraint)
9 x 1+ 13 x 2 ≤39000 (Budgetary constraint)
(∴ A+B these two lines are intersecting)
3x1 x 2 6000 …1
9 x 1+ 13 x 2=39000 …2
Multiply equation (1) by -3 on both sides:
-9x1 - 3x 2 -18000 …3
9 x 1+ 13 x 2=39000 …4
__ __ __
10X2 = 21,000 ∴ x 2=2100
Put the Value of x2 in first equation:
⇒X1=1300
At point (1300, 2100) Maximize Z=3 x 1+ 4 x 2
Z=3 ×1300+ 4 ×2100
= 12,300 which is the maximum value
Result
The optimal solution is:
No of units of product A=1300
No of units of product B=2100
Total profit, = 12, 300 which is the maximum
Example 2:
Suppose that a machine shop has two different types of machines; machine 1 and
machine 2, which can be used to make a single product .These machine vary in the
amount of product produced per hour, in the amount of labor used and in the cost of
operation.
Assume that at least a certain amount of product must be produced and that we would
like to utilize at least the regular labor force. How much should we utilize each machine
in order to utilize total costs and still meets the requirement? The resource used, the cost
and the required hour is given in the following table.
Resource used
Labor/hr 2 3 15
st
The 1 constraint is: 20 x 1+15 x 2 ≥1 00 change this in to 20 x 1+15 x 2=100 and find X and
Y intercept.
20
Put x 1=0 ⇒ x 2= ∴ t h e point is(0 , 20/3)
3
Put x 2=0 ⇒ x 1=5 ∴ t h e point is(5 , 0)
Draw the graph with x1 on x-axis & x2 on y-axis as shown in the following figure.
The 2nd constraint is 2X1+3X2≥15 change this in to 2X1+3X2=15 and find X and Y
intercept.
i. For ≥ constraint the feasible region will be the area, which lays above the constraint lines,
and for ≤ constraints, it will lays below the constraint lines. This would be useful in
identifying the feasible region.
ii. According to a theorem on linear programming, an optimal solution to a problem (if it
exists) is found at a corner point of the solution space.
X2
X1 =0
A (0, 20/3) 2X1 + 3X2 = 15
Feasible Region
AMU, CBE,5Department of Management December 2010 Page 42
B (2.5, 3.33)
X2 =0
Operations Research
Feasible solution
Corner Coordinates Objective function Value
point
Minimize . Z=25 X 1 +30 X 2
Use the same stapes to find the value of B by using the simultaneous equation as
presented previously.
Result
The optimal solution is:
X1 =2.5, X2=3.33 and Minimum cost Z= 162.5
• Alternative Optima
• Infeasible Solution
• Unbounded solution
Example1
Maximize z=2 x 1+ 4 x 2
Subject ¿:
x 1+2 x 2≤ 5
x 1+ x 2≤ 4∧x 1 , x 2≥ 0
The above Figure demonstrates how alternative optima can arise in LP model when the
objective function is parallel to a binding constraint. Any point on the line segment BC
represents an alternative optimum with the same objective value z = 10. Mathematically,
we can determine all the points (x 1, x2) on the line segment BC as a nonnegative
C: x1= 3, x2= 1
Observe that when α=0, (x1, x2) = (3, 1), which is point C. When α=1,
(x1, x2) = (0, 5/2), which is point B. For values of α between 0 and 1 (x 1, x2) lays
between B and C
Example 2:
_____________________________________________________________________
_____________________________________________________________________
Cutting 3 6 900
Assembly 1 1 200
_____________________________________________________________________
Assume that the company has a marketing constraint on selling products B and therefore
it can sale a maximum of 125 units of this product.
Required:
X2 X2=0
(0, 200)
FR X1=0
A (0, 0) X1
(200, 0) (300, 0)
X 1 +X 2 ≤200
Corners Coordinates Maximize Z=8 X1 + 16X2
A (0, 0) 0
Interpretation:
Both C and D are optimal solutions. Any point on the line segment CD will also lead to
the same optimal solution.
==>Multiple optimal solutions provide more choices for management to reach their
objectives.
2. Infeasible Solution
A solution is called feasible if it satisfies all the constraints and the constraints and non-
negativity condition. However, it is sometimes possible that the constraints may be
inconsistent so that there is no feasible solution to the problem. Such a situation is called
infeasibility.
Example:
Maximize Z=20X1+30X2
Subject to:
2X1+X2< 40
4X1+X2< 60
X1 > 30 and X1, X2 > 0
Solution:
X2 X1=0
(0, 60) X1=30
4X1+X2= 60
(0, 40)
2X1+X2= 40
X2=0
X1
(15, 0) (20, 0) (30, 0)
Note: -In the above graph, there is no common point in the shaded area.
-All constraints cannot be satisfied simultaneously and there is no feasible solution
to the problem.
3. Unbounded Solution
1. Maximize Z=3X1+4X2
Subject to:
X1-X2<-1==> -X1+X2>1 since the quantity solution is positive
-X1+X2<0
X2 X1, X2 > 0
X1-X2 =-1
X1+X2 =0
X1
2. Maximize Z=3X1+2X2
Subject to:
X1-X2<1
X1+X2≥3
X1, X2 > 0
X2
A (0, 3) Unbounded
Feasible Region
X1-X2=1
Note: here that the two corners of the region are A(0,3) and .B(2,1).The value of
Maximize Z(A)=6 and Maximize Z(B)=8. But there exist number of points in the shaded
region for which the value of the objective function is more than 8.For example, the point
(10, 12) lies in the region and the function value at this point is 70 which is more than 8.
Remark: An unbounded solution does not mean that there is no solution to the given
LPP, but implies that there exits an infinite number of solutions.
Exercise 1:
X2> 10
X1, X2 > 0
Exercise2.
-X1+X2< 5 X1 + 4 X2 > 12
II.A manufacturer produces two different models; X and Y, of the same product .The raw
materials r1 and r2 are required for production. At least 18 Kg of r1 and 12 Kg of r2 must
be used daily. Almost at most 34 hours of labor are to be utilized .2Kg of r1 are needed
for each model X and 1Kg of r1 for each model Y. For each model of X and Y, 1Kg of r2
is required. It takes 3 hours to manufacture a model X and 2 hours to manufacture a
model Y. The profit realized is $50 per unit from model X and $30 per unit from model
Y. How many units of each model should be produced to maximize the profit?
Answer: 10 units of model X, 2 units of model Y and the maximum profit is $ 560.
III.A manufacturing firm produces two machine parts P1 and P2 using milling and
grinding machines .The different machining times required for each part, the machining
times available on different machines and the profit on each machine part are as given
below:
____________________________________________________________________
Machine P1 P2
________________________________________________________________________
Lathe 10 5 25,000
_____________________________________________________________________
Dear learners, all of you have by now understand the graphical method of solving a
linear programming model and the various special cases of graphic solution of linear
programming model.
Let us see what the limitations of graphic solution are, and how can these be deal with?
Limitations of the Graphical Method
The applicability of the graphical method is very limited in scope. This is due to the fact that
it is quite simple to identity all the corner points & then tests them for optimality-in the case
of a two-variable problem. As a result, the graphical method cannot be always employed to
solve the real-life practical Linear programming models which involve more than two
decision-variables.
The above limitation of the graphical method is tackled by what is known as the simplex
method. Developed in 1940 by George B-Dantizg, it remains a widely applicable method for
solving complex LP problems. It can be applied to any LP problem which can be expressed
in terms of a Linear Objective function subject to a set of Linear Constraints. As such, no
theoretical restrictions are placed on the number of decision variables or constraints
contained in a LPP.
The Simplex method is an iterative or “step by step” method or repetitive algebraic
approach that moves automatically from one basic feasible solution to another basic
feasible solution improving the situation each time until the optimal solution is reached
at.
Note: The Simplex method starts with a corner that is in the solution space or feasible region and
moves to another corner, the solution space improving the value of the objective function each
time until optimal solution is reached at the optimal corner.
The use of the Simplex method to solve an LP problem requires that the problem be
converted into its standard form. The standard form of the LP problem should have the
following characteristics:
[ ]
a11 a12 … a1 n
a21 a22 … a2 n
And
… … …
am 1 am 2 amn
Are added in the given LP problem to convert it in to the standard form for the following
reasons:
Types of Extra variable needed Coefficient of extra variables In Presence of extra variable in
constraint the objective function initial solution mix
Max Z Min Z
variable is added
The summery of the extra variable to be added in the given LP problem to convert it in to
standard form is given in the following table.
Remark: a slack variable represents unused resource, either in the form of time on a
machine, labor hour, money, warehouse space or any number of such resources in various
business problems. Since variable yields no profit, therefore such variable are added to
the original objective function with zero coefficients.
A surplus variable represents amount by which solution values exceed a resource. These
variables are also called negative slack variables. Surplus variables carry a zero
coefficient in the objective function.
This method utilizes the property of a LP problem of having optimal solution only at the
corner point of the feasible solution space. It systematically generates corner point
solutions & evaluates them for optimality. The method stops when an optimal solution is
found. Hence, it is an iterative (repetitive) technique.
If we get more variables & less equation, we can set extra variables equal to zero, to
obtain a system of equal variables & equal equations. Such solution is called basic
solution.
The variables having positive values in a basic feasible solution are called basic variable
while the variables which are set equal to zero, so as to define a corner point are called
non-basic variables.
Slack variables are the fictitious variables which indicate how much of a particular
resource remains unused in any solution. These variables can’t be assigned negative
values. A zero value indicates that all the resources are fully used up in the production
process.
Z j Row denotes the contribution margin lost if one unit is brought into the solution.
Hence, it represents the opportunity cost. (Opportunity cost is the cost of sacrifice i.e., the
opportunity foregone by selecting a particular course of action out of a number of
different available alternatives).
C j−Z j Row denotes the Net Potential contribution or the Net unit margin potential, per unit.
The rules used under simplex method, for solving a linear programming problem are as
follows:-
1. Convert the LP to the following form:
Convert the given problem into Standard maximization Problem i.e. minimization
problem into a maximization one (by multiplying the objective function by -1). All
variables must be non-negative. All right hand side values must be non-negative
(multiply both sides by -1, if needed). All constraints must be in ≤ form (except the
non-negativity conditions). No strictly equality or ≥ constraints are allowed.
2. Convert all ≤ constraints to equalities by adding a different slack variable for each
one of them.
3. Construct the initial simplex tableau with all slack variables in the Basic Variables.
The last row in the table contains the coefficient of the objective function (rowC j).
4. Determine whether the current tableau is optimal. That is: If all the right hand side
values are non-negative (called the feasibility condition). If all elements of the last
row, that is C j rows are non-positive (called, the optimality condition)
If the answers to both of these two questions are yes, then stop. The current tableau
contains an optimal solution. Otherwise, go to the next step.
5. If the current BASIC VRIABLES is not optimal, determine, which non basic variable
should become a basic variable and, which basic variable should become a non basic
variable. To find the new BASIC VRIABLES with the better objective function value,
perform the following tasks:
Identify the entering variable: The entering variable is the one with the
largest positive C j value (In case of a tie, select the variable that corresponds
to the left most of the columns).
Identify the outgoing variable: The outgoing variable is the one with
smallest non-negative column ratio (to find the column ratios, divide the RHS
column by the entering variable column, wherever possible). In case of a tie
select the variable that corresponds to the upper most of the tied rows.
Generate the new tableau
a. Select the largest value of C j−Z j row. The column, under which this
value falls, is the pivot-column.
b. Pivot-row selection rule. Find the ratio of quantity to the corresponding
pivot-column coefficient. The pivot-row selected is the variable having
the least ratio.
c. The coefficient, which is in both, the pivot-row & the pivot column, is
called the pivot-element or pivot-number.
d. Up-dating Pivot-row. Pivot-row, also called replaced rows, are updated
as under all elements of old-row divided by Pivot-element. Now, in the
basic activities column, write the pivot-column variable in place of the
pivot-row variable. i.e.; the pivot-row variable is to be replaced by the
pivot-column variable.
e. Up-Dating all other rows. Update all other rows by updating the
formulae.
Up-dating Z j∧C j −Z j rows. Each Z j , is obtained as the sum of the products of the C j
column coefficients multiplied by the corresponding coefficient in the j th column. (i.e.)
the Quantity column). It is then subtracted from C j−Z j row values to get C j−Z j values.
This pivoting is to be repeated till no positive coefficients exist in the C j−Z j row, the
optimal solution is known.
Dear learners let us consider some examples to test our understanding of the solution
algorithm that has been discussed so far.
Example 1
Smart Limited manufactures two types of cement which are sold under the brand name
quick and Tuff. Each product consumes the same raw materials but in varying
proportions. The following table depicts the amount of raw materials along with thei r
respective cost.
A 400 50 100
P 400 50 100
1000kg 1000kg
Quick can be blended @ 1000 kg /hour where as the blending rate for Tuff is 1250
kg/hour. Their respective selling prices are Birr. 1010 & Birr 845 You may assume the
variable costs to be Birr 500 per hour of plant production time. The maximum availability
of raw materials is:
N 1000
A 300
P 250
I 1800
Formulate as a linear programming model and find out the optimal units of quick & tuff
to be produced so as to maximize the profits.
Solution
Step I List the objective & constraint equations.
Step II Introduce the slack variables
Step III Arrange in the form of initial table.
Step IV Find out the profit-margins from given sales price.
Step V Generate solutions.
The detailed solutions are as under:
Step I List the objective & constraint equations.
Xj ≥0 for all j
Remarks
2. Contribution margin represents the profit which remains if after deducting variable costs
from sales i.e.; It covers fixed costs & net profit i.e. it fixed cost = Birr 200 X five – 100
kg of quick are sold then the net profit = 5 x 150 –200 = Rs. 550
Now, we have to find those values of x 1 & x2 for which the contribution is maximum.
Here, the constraint is the availability of raw material. Hence, the problem is formulated
as;
Maximize Z= 150 x 1+125 x 2
Subject to: 350 x 1+200 x 2 ≤1000
x1 , x2 ≥ 0
==>n=5 variables (x1 , x2, s1, s2, and s3) and m=3 constraints (Labor, machine and
marketing constraints), excluding non-negativity.
Therefore, n-m=5-3=2 variables(x1 and x2) are set equal to zero in the 1 st simplex tableau.
These are non-basic variables. 3 Variables (s1, s2, and s3) are basic variables (in the 1st
simplex tableau) because they have non-zero solution values.
Converting inequalities into equalities by using slack-variable
Maximize Cont-margin = 150 x 1+125 x 2
Subject to: 350 x 1+200 x 2 + s1=1000
50 x 1+100 x 2 +s 2=250
550 x 1+600 x 2 + s3=1800
si X j≥ 0 for all i∧ j
Note. We have dropped: 50 x 1+100 x 2 ≤ 300 as it is already contained in
50 x 1+100 x 2 ≤ 250
Slack can’t be larger than the constant. On right hand side (RHS) (i.e. S l=1000 , S 2=250
etc) It is to be so, and then some other variable must be negative for equality to exist, which
is not possible.
(2) Rewriting as:
Maximize Cont-margin = 150 x 1+125 x 2 +0 s 1 +0 s 2 +0 s 3
Subject to: 350 x 1+200 x 2 + s1 +0 s 2 +0 s 3=1000
50 x 1+100 x 2 +0 s 1 +s 2 +0 s 3=250
550 x 1+600 x 2 +0 s 1 +0 s 2 +s 3=1800
si∧ X j ≥ 0 for all i∧ j
Zj(Rs .) 0 0 0 0 0 0
0 S2
0 S3
Old row element−Corresponding Pivot column element × Updated corresponding Pivot row element=New
50 - (50 x 1) = 0
1- (50 x 0) =1
0- (50 x 0) =0
550 - (550 x 1) = 0
0- (550 x 0) =0
1- (550 x 0) =1
Zj is calculated as under :
Zj value for QTY Column=150 x 2.857+0 x 107.145+ 0 x 228.595=428.55
For x 1=150 x 1+0 x 0+0 x 0=150
x 2=150 x 0.5714 +0+0=85.71
S 1=150 x 0.0029+ 0+0=0.435
S 2=0
S 3=0
Entering
variable
Completed 3rd Tableau
Birr x1 x2 S1 S2 S3
0 S2 50 0 0 0.225 1 -0.25
Note
(1)x 1row is up−dated as :
2.857−¿
1−(0.5714 x 0)=1
Since no positive co-efficient exists in the C1 – z1 row, this is the optimal solution.
X1 = 2.4 X2 = 0.8
Z1= 460.
Example 2
X Ltd. Produces two products p1, p2 having profit of Rs. 4 $ Rs. 3 each p 1, p2 require 4
hrs. & 2 hours of machining respectively, the total available machining time is 10 hours.
P1, p2 consume 2 units & 8/3 units of raw material respectively subject to a total of
maximum 8 units. Any no. of p 2 j can be produced & sold but the no. of p 1 must not be
Solution
Maximise=4 x 1+3 x 2+ os 2+os 3
Subject ¿:
8
2 x 1+ x 2 +0 s 1 +s 2+ 0 s 3=8 x 1 ; x 1 , x 2 ≥ 0
3
0 S1 10 4 2 1 0 0 10/4= 2.5
0 S2 8 2 8/3 0 1 0 8/2= 4
0 S3 6 1 0 0 0 1 6/1= 6
Zj(Birr) 0 0 0 0 0 0
Cj−Zj(Birr) 4 3 0 0 0
Pivot Column = x1 because in the Cj−Zj raw x1 have the highest value.
( 52 )=3
8− 2 × ( 52 )= 72
8− 2 ×
5
Qty.4× =10
2
8
3 ( )
1 5
− 2× =
2 3 ( )
0− 2 ×
1
2
=−1
1
X2=4× =2
2
( 41 )=−12
0− 2 × ( 14 )=−14
0− 1 ×
1
S1=4× =1
4
Second tableau:
s1 Qty./ x1
Birr x1 x2 s2 s3
Zj(Birr) 10 4 2 -1/4 1 0
Cj−Zj(Birr) 0 1 1 1 0
Pivot row = s2 the smallest value in the min ratio column (9/5). Don’t take negative and
zero value.
(
5 1 9 8
− × =
2 2 5 5 ) 2
− (
7 −1 9 22
× =
2 5 5 ) 8 3× 9
Qty.4× +
5 5
=11.
6
1
( )
1− ×0 =1
2
1
(
0− 0×− =0
2 )
1 1
( )
− ×1 =0
2 2
−1 −1
2
−
2
×1 ( )
(
1 1 1 2
− × =
4 2 4 5 ) 4
− (
−1 −1 −3 −2
2
×
4
=
5 )
0− ( 12 × 35 )= −310 0− ( −12 × 35 )= 103
0− ( 12 × 0)=0 1− ( −12 ×0)=1
3rd Tableau
Cj (Birr) Basic activities Qty. 4 3 0 0 0
Birr x1 x2 s1 s2 s3
3 x2 9/5 0 1 ¼ - 3/5
3/10
Cj−Zj(Birr) 0 0 1 - -3/5
7/10
Example 3:
Subject ¿:
2 x 1+ x 2< 40(Labor )
x 1<12( Marketing)
x 1 , x 2> 0
Solution
i.e. Convert constraint inequality into equality form by introducing a variable called
Sack variable.
Slack Variables:
A sack variable(s) is added to the left hand side of a < constraint to covert the constraint
inequality in to equality. The value of the slack variable shows unused resource.
A slake variable emerges when the LPP (linear programming problem) is a maximization
problem.
Slack variables represent unused resource or idle capacity. Thus, they don’t produce any product
and their contribution to profit is zero.
Slack variables are added to the objective function with zero coefficients.
Let that, s2 and s3 be unused labor, machine and marketing hours, respectively.
Let that s1, s2, and s3 are unused labor, machine and marketing hours respectively.
Subject to:
2 x1 + x 2 +s 1 +0 s 2 +0 s 3=40
x 1+ 3 x 2 +0 s1 + s2 +0 s 3=45
x 1+ 0 s 1+ 0 s 2+ s 3=12
x 1 , x 2 , s1 , s 2 , s3 >0
To represent the data, the simplex method uses a table called the simplex table or the
simplex matrix.
==> In constructing the initial simplex tableau, the search for of the optimal solution
begins at the origin. Indicating that nothing can be produced;
==> x 1+ 0 s 1+ 0 s 2+ s 3=12
0+ 0 s 1+0 s2 + s3 =12
Note:
In general, whenever there are n variables and m constraints (excluding the non-
negativity), where m is less than n (m<n), n-m variables must be set equal to zero before
the solution can be solved algebraically.
Or: basic variables are variables that are in the basic solution. Basic variables have 0
values in the C j−Z j row.
Or: non-basic variables are variables that are out of the solution.
==>n=5 variables (x1 , x2, s1, s2, and s3) and m=3 constraints (Labor, machine and
marketing constraints), excluding non-negativity.
Therefore, n-m=5-3=2 variables(x1 and x2) are set equal to zero in the 1 st simplex tableau.
These are non-basic variables. 3 Variables (s1, s2, and s3) are basic variables (in the 1st
simplex tableau) because they have non-zero solution values.
Step 3 Construct the initial simplex tableau. Dear learners alternatively you may solve LP
simplex method by using the following method (you may solve the problem by using the
following table arrangement).
Cj 300 250 0 0 0
Profit per unit
Row
SV X1 X2 S1 S2 S3 Q
0 S1 2 1 1 0 0 40 R1
0 S3 1 0 0 0 1 12 R3
Gross Profit row
Zj 0 0 0 0 0 0
Net Profit row /Indicato
Cj - Zj 300 250 0 0 0
Note:
The entering variable is the variable that has the most positive value in the C j−Z j row
also called as indicator row. Or the entering variable is the variable that has the highest
contribution to profit per unit.
==> In this step, we determine the variable that will leave the solution for X1 (or entering
variable)
Note: The row with the minimum or lowest positive (non-negative) replacement ratio
shows the variable to leave the solution.
Note: RR>0
The variable leaving the solution is called leaving variable or outgoing variable.
The row associated with the leaving variable is called key or pivot row (s 3 column
in our case)
The element that lies at the intersection of the pivot column and pivot row is
called pivot element(No 1 in our case)
Step 6: Repeat step 3-5 till optimum basic feasible solution is obtained.
Or: repeat step 3-5 till no positive value occurs in the C j−Z j row.
Note:
Divide each element of the pivot row by the pivot element to find new values in the
key or pivot row.
Perform row operations to make all other entries for the pivot column equal to
zero.
Cj 300 250 0 0 0
SV X1 X2 S1 S2 S3 Q
0 S1 0 1 1 0 -2 16
R’1=R1-2R2
0 S2 0 3 0 1 -1 33
R’2=R2-R3
300 X1 1 0 0 0 1 12
R’3=R3
Cj - Zj 0 250 0 0 -300
Cj 300 250 0 0 0
SV X1 X2 S1 S2 S3 Q
0 S1 0 0 1 -1/3 -5/3 5
R’’1=R’1-R’2
250
R’’2=R2/3
X2 0 1 0 1/3 -1/3 11
R’’3=R’3
300 X1 1 0 0 0 1 12
Cj - Zj 0 0 0 -250/3 - 650/3
2. Conversion method
Example:
2. 5x1+3x2 < 45
3. 5x1+2x2 >20
4. 2x1+x2 >40
Thus, in order to avoid the mathematical contradiction, we have to add artificial variable (A)
Artificial variable is a variable that has no meaning in a physical sense but acts as a tool to create
an initial feasible LP solution.
Note:
The Big-M Method is a method which is used in removing artificial variables from the
basis .In this method; we assign coefficients to artificial variables, undesirable from the
objective function point of view. If objective function Z is to be minimized, then a very
large positive price (called penalty) is assigned to each artificial variable. Similarly, if Z
is to be maximized, then a very large negative price (also called penalty) is assigned to
each of these variables.
Example:
Subject to:
x1, x2 >0
Solution
Subject to:
The initial basic feasible solution is obtained by setting x1= x2= s1= s2=0
Cj 25 30 0 0 M M
SV X1 X2 S1 S2 A1 A2 Q
RR (min ratio)
M A1 20 15 -1 0 1 0 100
M A2 2 3 0 -1 0 1 15 100/20=5
X1 is entering variable because in the C j – Z j row X1 have the highest negative value. And
A1 is leaving or outgoing variable because the RR is lowest positive in the A 1 row.
Note:
Once an artificial variable has left the basis, it has served its purpose and can therefore
be removed from the simplex tableau. An artificial variable is never considered for re-
entry into the basis.
Cj 25 30 0 0 M
SV X1 X2 S1 S2 A2 Q
M A2 0 3/2 1/10 -1 1 5
R’2=R2-2
Zj 25 75/4+3/2M -5/4+1/10M -M M 125+5 M
Cj - Zj 0 45/4-3/2M 5/4-1/10 M M 0
Cj 25 30 0 0
SV X1 X2 S1 S2 Q
Cj - Zj 0 0 1/2 15/2
X1=5/2
Note:
As long as an “A” variable is available in the solution variable column, the solution is
infeasible.
Example 2. Use the penalty (Big-M) method to solve the following LPP
Subject to:
2x1+4x2 < 12
2x1+ 2x2 = 10
x1, x2 >0
Solution
Cj 5 3 0 0 M M
SV X1 X2 S1 S2 A1 A2 Q
RR
0 S1 2 4 1 0 0 0 12 6
5
M A1 2 2 0 0 1 0 10
M A2 5 2 0 -1 0 1 10 2
Zj 7M 4M 0 M M M 20 M
Cj - Zj 5 -7M 3- 4M 0 -M 0 0
Cj 5 3 0 0 M
SV X1 X2 S1 S2 A1 Q
0 S1 0 16/5 1 2/5 0 8
nd
2 simplex tableau
M A1 0 6/5 0 2/5 1 6
5 X1 1 2/5 0 -1/5 0 2
Cj - Zj 0 -6/5M +1 0 -2/5M+1 0
Operations Research
Cj 5 3 0 0 M
SV X1 X2 S1 S2 A1 Q
RR
3 X2 0 1 5/16 1/8 0 2.5
20
M A1 0 0 -3/8 1/4 1 3
12
5 X1 0 0 -1/8 -1/4 0 1
-
Zj 5 3 -3/8M +5/6 M/4-7/8 M 12.5+3 M
Cj 5 3 0 0
SV X1 X2 S1 S2 Q
3 X2 0 1 1/2 0 1
0 S2 0 0 -3/2 1 12
5 X1 0 0 -1/2 0 4
Zj 5 3 -1 0 23
Cj - Zj 0 0 1 0
Exercise
Subject ¿:
2 x 1+5 x 2 ≥ 150
3 x 1+ x 2≥ 120
x1, x2≥0
Subject ¿:
x 1+2 x 2> 80
3 x 1+ x 2>75
x 1 , x 2> 0
1. Mixed constraints
Example
Subject to:
x2 < 4
x1+ x2 = 9
x1, x2 >0
Standard form
St:
x2 + s1 =4
x1+ x2 + A2 =9
Standard form
6x1+2x2 - s3 + A3 =24
Cj 6 8 0 0 -M -M
SV X1 X2 S1 S3 A2 A3 Q
0 S1 0 1 1 0 0 0 4
-M A2 1 1 0 0 1 0 9
-M A3 6 2 0 -1 0 1 4
Zj -7M -3M 0 +M -M -M 24
Cj - Zj 7M +6 3M+8 0 -M 0 0
Note:
For the initial basis, use artificial variables for constraints that have them. Otherwise, use
a constraint slack variable. Hence, surplus variables will not appear in an initial solution.
In order to break this tie, the selection for the key column (entering variable) can be made
arbitrary. However; the number of solution can be minimized by adopting the following
rules:
1. If there is a tie between two decision variables, then the selection can be made
arbitrary.
2. If there is a tie between a decision variable and a slack (or surplus) variable, then
select the decision variable to enter into basis first.
3. If there is a tie between slack or surplus variable, then selection can be made
arbitrary.
If a non- basic variable corresponding to which C j−Z j =0 the problem have multiple
optimal solution.
Example1
As before, we add slacks, and we solve by the Simplex method, using tableau representation.
Max Z =6 x 1+ 4 x 2
Subject to: 2 x1 +3 x 2 ≤ 30
3 x 1+2 x 2 ≤ 24
x 1+ x2 ≥ 3
x 1∧x 2 ≥ 0
As before, we add slacks s1, s2 and surplus s 3∧artificial A 1, and we solve by the simplex
method, using tableau representation.
3 x 1+2 x 2 + s 2=24
x 1+ x2 + A1 =3
x 1 , x 2 , s1 , s 2 , s3 + A1 ≥ 0
The optimal solution for this LP problem is presented in the following table
Cj 6 4 0 0 0 RHS
BV x1 x2 s1 s2 s3
0 s1 0 5/3 1 - 0 14
2/3
0 s3 0 -1/3 0 1/3 1 5
6 x1 1 2/3 0 1/3 0 8
Zj 6 4 0 2 0 48
C j−Z j 0 0 0 -2 0
The optimal solution is x1=8, x2=0, s1=14, and s3= 5 and max Z=48
But the C j−Z j row shows a value zero corresponding to a non basic variable X 2. Thus
an alternative optimal solution can also be obtained by entering variable X 2 in to the basis
and removing S1 from the basis.
Alternative solution
Cj 6 4 0 0 0 RHS
BV x1 x2 s1 s2 s3
Zj 6 4 0 2 0 48
C j−Z j 0 0 0 -2 0
The optimal solution is x1=12/5, x2=42/5, s1=0, and s3= 39/5 and max Z=48
3. Unbounded solution
In maximization LP problem, if C j−Z j >0 (C j−Z j <0 for a minimization case) for a
column not in the basis and all entries in this column are negative, then for determining
key row, we have to calculate minimum ration corresponding to each basic variable have
negative or zero value in the denominator.
Maximize Z= 3 x 1+5 x 2
Subject to x 1−2 x2 ≤ 6
x 1 ≤ 10
x2≥ 1
x 1∧x 2 ≥ 0
add slack s1, s2 and surplus s 3∧artific ial A 1, and we solve by the simplex method, using
tableau representation.
Subject to x 1−2 x2 + s1 =6
x 1 ≤ s 2=10
x 2−s 3+ A 1=1
x 1 , x 2 , s1 , s 2 , s3∧ A 1 ≥ 0
Cj 3 5 0 0 0 -M RHS Min
ratio
BV x1 x2 s1 s2 s3 A1
0 s1 1 -2 1 0 0 0 6 ___
0 s2 1 0 0 1 0 0 10 ___
-M A1 0 1 0 0 -1 1 1 1
Zj 0 -M 0 0 M -M 48
C j−Z j 3 5+M 0 0 -M 0
Cj 3 5 0 0 0 -M RHS Min
ratio
BV x1 x2 s1 s2 s3 A1
0 s1 1 -2 1 0 0 0 6 ___
0 s2 1 0 0 1 0 0 10 ___
-M A1 0 1 0 0 -1 1 1 1
Zj 0 -M 0 0 M -M -M
C j−Z j 3 5+M 0 0 -M 0
Second table
Cj 3 5 0 0 0 -M RHS Min
BV x1 x2 s1 s2 s3 A1 ratio
0 s1 1 0 1 0 -2 2 6 ___
0 s2 1 0 0 1 0 0 10 ___
5 x2 0 1 0 0 -1 1 1 1
Zj 0 5 0 0 -5 5 5
C j−Z j 3 0 0 0 5 -M-5
Variable s3 should inter to the basis but it may be noted that coefficient in the S 3 column
ar all negative or zero. This indicates that S 3 cannot be entered in to the basis. However,
the value of S3 can be increased infinitely without removing any one of the basic variable.
Infeasible solution
In final simplex table if at least one artificial variable appear with a positive value, no
feasible solution exists, because it is not possible to remove such an artificial variable
from the basis using the simplex algorithm.
Max Z =6 x 1+ 4 x 2
Subject to: x 1+ x2 ≤5
x2 ≥ 8
x 1∧x 2 ≥ 0
x 2 + s 2+ A 1=8
x 1 , x 2 , s1 , s 2∧ A1 ≥ 0
Cj 6 4 0 0 -M RHS
BV x1 x2 s1 s2 A1
4 x2 1 1 1 0 0 5
-M A1 -1 0 -1 -1 1 3
Since the C j−Z j ≤ 0 the solution is optimal. But the solution not feasible for the given
problem since it has X1=0, X2=5 (recall that in the second constraint X 2 ≥ 8). The fact that
artificial variable A1 =3 is in the solution also indicates that the final solution violates the
second constraint (X2 ≥ 8) by 3 unites.
Introduction
For every LP formulation there exists another unique linear programming formulation
called the 'Dual' (the original formulation is called the 'Primal'). The - Dual formulation
can be derived from the same data from which the primal was formulated. The Dual
formulated can be solved in the same manner in which the Primal is solved since the Dual
is also a LP formulation.
The Dual can be considered as the 'inverse' of the Primal in every respect.
The columns coefficients in the Primal constrain become the row coefficients in
the Dual constraints.
The coefficients in the Primal objective function become the right hand side
constraints in the Dual constraints.
The column of constants on the right hand side of the Primal constraints becomes
the row of coefficients of the dual objective function.
The 'directions of the inequalities are reversed.
If the primal objective function is a 'Maximization' function then the dual
objective function is a 'Minimization' function and vice-versa.
The concept of duality is very much useful to obtain additional information about the
variation in the optimal solution when certain changes are effected in the constraint
coefficient, resource availabilities and objective function coefficients. This is termed as
post optimality or sensitivity analysis.
Note:
The dual contains economic information useful to management, and it may also be
easier to solve, in terms of less computation, than the primal problem.
Corresponding to every LPP, there is another LPP.
Table
Primal-Dual Relationship
Primal Dual
Objective is minimization Objective is maximization and vice versal
No of columns No of rows
No of rows No of columns
c x
Duality Advantage
2. The dual reduces the computational difficulties associated with some formulation
Dual Formation
Following are the steps adopted to convert primal problem into its dual.
Step1. For each constraint in primal problem there is an associated variable in dual
problem.
Step2. The elements of right hand side of the constraints will be taken as the co-
efficient of the objective function in the dual problem.
Step3. If the primal problem is maximization, then its dual problem will be
minimization and vice versa.
Step4. The inequalities of constraints should be interchanged from ≥ to ≤ and vice versa
and the variables in both the problems and non-negative.
Step5. The rows of primal problem are changed to columns in the dual problem. In other
words the matrix A of the primal problem will be changed to its transpose (A) for the
dual problem.
Step6. The co-efficient of the objective function will be taken the right hand side of the
constraints of the dual problem.
An example will clarify the concept basis. Consider the following 'Primal' LP formulation
Example 1
Mathematical formulation the primal and the dual linear programming problems of the
given problem are
1. The column coefficient in the Primal constraint namely (2,2) and (3,1) have become the
row co-efficient in the Dual constraints.
2. The co-efficient of the Primal objective function namely, 12 and 10 have become the
constants in the right hand side of the Dual constraints.
3. The constants of the Primal constraints, namely 18 and 14, have become the coefficient in
the Dual objective function.
4. The directions of the inequalities have been reserved. The Primal constraints have the in
equalities of ≤; while the Dual constraints have the inequalities of ≥.
5. While the Primal is a 'Maximization' problem the Dual is a 'Minimization' problem
Example: 2
a. Maximize Z=5x1+6x2
Subject to:
x1, x2 > 0
Solution
u1 2 3 < 3000
u2 5 7 < 1000
u3 1 1 < 500
Maximize Z 5 6
By referring the above table, dual for this can be stated as:
Subject ¿:
3 u 1+ 7 u 2+u 3> 6
u 1 ,u 2 , u 3>0
Note:
Example3.
_____________________________________________________________________
A B
I 30 20 300
II 5 10 110
Profit/unit ($) 6 8
i. Formulate the LPP model for the primal and determine the optimal solution using simplex
method.
ii. Formulate the LPP model for the dual and determine the optimal solution from the final
primal simplex tableau
Solution
Subject ¿:
30 X 1 +20 X 2 ≤300
5 X 1 +10 X 2 ≤110
X 1 , X 2 ≥0
Cj 6 8 0 0
SV X1 X2 S1 S2 Q
6 X1 1 0 1/20 - 4
1/10
8 X2 0 1 -1/40 9
3/20
6 8 -1/10
Zj 96
6/10
0 0 -1/10 -
C j - Zj
6/10
Subject ¿:
30 u1 +5 u2 ≥ 6
20 u1 +5 u2 ≥ 8
u1 ,u 2 ≥ 0
NOTE:
The value in the Cj - Zj row under columns of the slack /surplus variables with change in
sign give directly the optimal values of the dual/primal basic variables.
Therefore,
u j=−(C j−Z j)
Note: The column s1 corresponds to u1 an the column s2 corresponds to u2
Note: The dual variables uj are also called as the shadow prices
Exercise1.
Subject to:
2 X1+ X2 <40
X1 <12
X1 , X2 > 0
Cj 300 250 0 0 0
SV X1 X2 S1 S2 S3 Q
0 S1 0 0 1 -1/3 -5/3 5
300 X1 1 0 0 0 1 12
Cj - Z j 0 0 0 -250/3 -650/3
Required:
I. What are the solution of the dual variables (i.e. Shadow prices) u1, u2 and u3?
II. What is the optimal dual cost?
Solution:
u j=−(C j−Z j)
Sensitivity Analysis is concerned with the study of ‘Sensitivity ‘of the optimal solution of
an LPP with discretion variables (changes) in parameters .The degree of sensitivity of the
solution due to those variations can range from no change at all to a substantial change in
the optimal solution of the given LPP. Thus, insensitivity analysis, we determine the
range over which the LP model parameters can change without affecting the current
optimal solution. The process of studying the sensitivity of the optimal solution of an
LPP is called post-optimal analysis.
Five types of discrete changes in the original LP model may be investigated during the
sensitivity analysis:
Note:
Instead of resolving the entire problem as a new problem with new parameters, we may
take the original optimal solution table as an initial solution table for the purpose of
knowing ranges both lower and upper within which a parameter may assume value.
The range of optimality is the range over which a basic decision variable coefficient in
the objective function can change without changing the optimal solution mix. However,
this change will change only the optimal value of the objective function.
Example:
Subject to:
5x1+6.4x2+15x3 < 77
Cj 5 4.5 1 0 0 0
SV X1 X2 X3 S1 S2 S3 Q
5 X1 1 1.053 0 0.067 0 0 10
Determine the range of optimality for the coefficients of the basic-decision variables.
Solution:
The analysis will be conducted on products on X1 and X3 which are in the basic solution.
Divide each Cj - Zj row entry for variables not in the solution (for instance, by X2, S1 and
S2 values) by the associated variable aij from X1or X3 row.
I. Analysis of X1
Steps:
a. Take the Cj - Zj row of the optimal solution of the non-basic variables
b. Take the X1 row of the non-basic variables
c. Cj - Zj row
X1 row
X2 S1 S2 Non-basic variables
X1 row
Upper Limit
The smallest positive number in the Cj - Zj row tells how much the profit X1 of X1 can
be increased before the solution is changed.
Cj - Zj row
-1.26 14.13 -1
X1 row
=5+∞=∞ X1 row
Lower Limit
Therefore, the range of optimality for the coefficient of X1 is 4.2< Cj (for X1) < ∞ (The
coefficient of X1 in the objective function can change between 4.2 and ∞ without changing
the optimal solution mix X1=10, X3=1.8 and S3=15.12)
II. Analysis of X3
Therefore, the range of optimality for the coefficient of X3 is 0 < Cj (for X3) < 15.13 (The
coefficient of X3 in the objective function can change between 0 and 15.13 without
changing the optimal solution mix X 1=10, X3=1.8 and S3=15.12).However, this change
will change only the optimal value of the objective function(i.e.MaxZ will change)
Exercise:
Max.Z=50x1 +120x2
Subject to:
2 x1+4x2 < 80
3x1+x2< 60
x1, x2 > 0
Determine the range of optimality for the coefficient of the basic variables.
Optimal Solution
Cj
SV X1 X2 S1 S2 Q
5 X1 1/2 1 1/4 0 20
0 S2 5/2 0 -1/4 1 40
Zj 60 120 30 0 $2,400
$-10 $0 $-30
Cj - Z j
$0
Ans: The range of optimality for X2’s profit coefficient is: $100 < Cj (for X2) < ∞
If there is a variable Cj, not participating in the optimal basis, then, in order for the
variable to be included in the optimal solution, its coefficient in the objective function
will have to change from the existing Cj to a new level Cj(new).
Cj(new)> Zj
The range of insignificance is the range over which Cj rates for non-basic variables can
vary without causing a change in the optimal solution mix (variable) is called the range of
insignificance.
Example:
Subject to:
5x1+6.4x2+15x3 < 77
2.8x2+11.8x3 < 36
x1, x2 , x3 > 0
Cj 5 4.5 1 0 0 0
SV X1 X2 X3 S1 S2 S3 Q
5 X1 1 1.053 0 0.067 0 0 10
Solution
Cj(new for X2)>5.342==> Cj(new for X2).If the profit contribution of X2 is greater than
5.342,then X2 will be included in the solution.
Thus, ∞< Cj(new for X2)< 5.342 is the range of insignificance for X2.
Cj (new for X2) can vary with in the given range without causing a change in the optimal
solution mix(X1=10, X1=0, X3=1.8, S1= S2=0 and S3=15.12).
Or
Shadow prices:
Shadow prices signify the change in the optimal value of the objective function for 1 unit
increases in the value of the RHS of the constraint that represent the availability of scarce
resources.
The negative of the number of Cj - Zj row in its slack variable columns provide as with
shadow prices. Or: shadow prices are found in the Zj row of the final simplex tableau in
the slack variable columns.
RHS ranging is the range over which shadow prices remain valid.
Example:
Max.Z=3x1+4x2
Subject to:
3 x1+5x2 < 15
2x1 + x2 < 8
x2 < 2
x1, x2 > 0
Cj 3 4 0 0 0
SV X1 X2 S1 S2 S3 Q
Cj - Zj 0 0 -0.714 -0.428 0
Required:
Q S1 Q/ S1
=15-3=12
=15-(-3.99) =18.99
Therefore, 12< b1< 18.99 (The range of resource 1 over which the shadow price $0.714
per unit is valid).
Q S2 Q/ S2
3.57 0.714 5
1.143 0.428 2.67
0.857 -0.428 -2
=8-(2.67)=5.33
=8-(-2) =10
Therefore, 5.33< b1< 10 (The range of resource 2 over which the shadow price $0.428
per unit is valid).
Q S3 Q/ S3 Lower
Limit=b 3-the
=2-
0.857 0 -
(1.143)= 0.857
=2-(-∞) =∞
Therefore, 0.857< b3< ∞ (The range of resource 3 over which the shadow price $0 per
unit is valid).
Exercise
1. Solve the following LPP using simplex method. A firm that manufactures both lawn
mowers and snow blowers:
Max.Z=30x1+80x2
Subject to:
x1, x2 > 0
Answer:
b. What are the shadow prices? When the optimal solution has been reached, which resource
has the highest marginal value?
Answer:
The shadow price for 1 additional snow blowers engine made available =$20
Thus, snow blower engine have the highest marginal value at the optimal solution.
c. Over what range in each of the RHS values are these shadows valid?
Answer:
The shadow price for labor hours is valid from 800 hours to 1,066.66 hours
The shadow price for pounds of steel is valid from 1,000pounds up to an infinite number
of pounds. The shadow price for snow blower engines ranges from 180 engines up to 250
engines
d. What are the ranges over which the objective function coefficients can vary for each of the
two decision variables?
Answer:
Without changing the current solution mix, the profit coefficient for the mowers can
range from $0 to 40, while the coefficient for the blowers can range from $60 to infinity.
Chapter Summary
The basic steps in formulating a linear programming model are presented as follows:
Step II Identification of the constraints. All the constraints in the given problem which
restrict the operation of a firm at a given point of time must be identified in this stage.
Step III formulate the objective function. In the last stage, the objective which is
required to be optimized (i.e., maximized or memorized) must be dearly identified and
expressed in terms of the pre-defined decision variables.
There are two types of finding a solution for Linear programming problems this are
Graphic solution and Simplex solution method.
The graphic solution procedure is one of the methods of solving two variable linear
programming problems. Steps in graphic solution method:-
Step I Defining the problem. Formulate the problem mathematically. Express it in terms
of several mathematical constraints & an objective function.
Step II Plot the constraints graphically. Each inequality in the constraint equation has to
be treated as an equation.
Step III Locate the solution space. Solution space or the feasible region is the graphical
area which satisfies all the constraints at the same time. Such a solution point (x, y)
always occurs at the corner Points of the feasible Region the feasible region is determined
as follows:
a. For "greater than" & "greater than or equal to" constraints, the feasible region or the
solution space is the area that lays above the constraint lines.
b. For" Less Then" &" Less than or equal to" constraint, the feasible region or the
solution space is the area that lays below the constraint lines.
Step IV Selecting the graphic solution technique. Select the appropriate graphic
technique to be used for generating the solution.
Optimal feasible solution: Any basic feasible solution which optimizes (i.e. maximize
or minimizes) the objective function of a LP models is known as the optimal feasible
solution to that linear programming model.
Degenerate Solution: A basic solution to the system of equations is termed as
degenerate if one or more of the basic variables become equal to zero.
The second method to find a solution for LPP is using the Simplex method. It is an
iterative or “step by step” method or repetitive algebraic approach that moves
automatically from one basic feasible solution to another basic feasible solution
improving the situation each time until the optimal solution is reached at.
The summery of the extra variable to be added in the given LP problem to convert it in to
standard form is given in the following table.
7. LP is a?
A. Constrained optimization technique
B. Technique for economic allocation of limited resources
C. Mathematical technique
D. All
8. Which of the following is limitation of LP
A. Linear Relationship
B. Constant Value of objective & Constraint Equations.
C. No Scope for Fractional Value Solutions.
D. Degree of Complexity.
E. Multiplicity of Goals.
F. All
9. The role of artificial variable in the simplex table is?
A. To aid in finding an initial solution
B. To obtain optimal dual prices in the final simplex table
C. To start phases of simplex method
D. All
10. For a minimization problem, the objective function coefficient for an artificial
variable is?
A. –M C. Zero
B. +M D. None
11. A variable which does not appear in the basic variable column of simplex table is ?
12. We must add artificial variable to write the standard form of LPP is?
15. the RHS constant of a constraint in a primal problem appears in the corresponding
dual as
A. a coefficient in the objective function
B. a RHS constant of constraint
C. an input- output coefficient
D. none
16. shadow price indicates how much one unit change in the resource value will change
the
A. optimality range of an objective function
B. optimal value of the objective function
C. value of basic variables in the optimal solution
D. none
20. the entering variable in the sensitivity analysis of objective function coefficient is
always
A. decision variable
B. non basic variable
C. basic variable
D. slack variable
E. none
21. while performing sensitivity analysis, the upper bound infinity on the value of the
RHS of the constraint means that
A. there is slack in constraint
B. the constraint is redundant
C. the shadow price for that constraint is zero
D. None
-X1+X2< 5 X1 + 4 X2 > 12
3. A manufacturer produces two different models; X and Y, of the same product .The raw
materials r1 and r2 are required for production. At least 18 Kg of r1 and 12 Kg of r2 must
be used daily. Almost at most 34 hours of labor are to be utilized .2Kg of r1 are needed
for each model X and 1Kg of r1 for each model Y. For each model of X and Y, 1Kg of r2
is required. It takes 3 hours to manufacture a model X and 2 hours to manufacture a
model Y. The profit realized is $50 per unit from model X and $30 per unit from model
Y. How many units of each model should be produced to maximize the profit?
4. A manufacturing firm produces two machine parts P1 and P2 using milling and grinding
machines .The different machining times required for each part, the machining times
available on different machines and the profit on each machine part are as given below:
____________________________________________________________________
Machine P1 P2
________________________________________________________________________
Lathe 10 5 25,000
_____________________________________________________________________
CHAPTER THREE
Introduction
All linear programming can be solved by the Simplex method. However due to their
specialized structure, certain classes of linear programming problems lend themselves to
solutions by other techniques that are computationally more efficient than the simplex
methods. These types of linear problems can be represented by the so called
transportation and assignment model. Let as first describe the transformation model and
the assignment model will be discussed later in this chapter.
Basic transportation problem was stated by F.L. Hitchcock in 1941, and later expanded
by T.C. Koopmans and W.W. Cooper. Several extensions of transportation model and
methods have been subsequently developed. The transportation problem is a particular
case of linear programming problem.
The structure of transportation problem involves a large number of shipping routes from
several supply origins to several demand destinations. The objective is to determine the
number of units which should be shipped from an origin a destination in order to satisfy
the required quantity of goods or services at each demand destination, within the limited
quantity of goods or services available at each supply origin, at the minimum
transportation cost and/or time.
The transportation algorithm discussed in this chapter is applied to minimize the total
cost of transporting a homogeneous commodity (product) from supply origin to demand
destination. However, it can also be applied to the maximization of some total value or
utility, for example, financial resources are distributed in such away that the profitable
Example1.
A company has three production facilities S1, S2, and S3 with production capacity of 7,
9, and 18 units (in 100s) per week of product respectively. These units are to be shipped
to four ware houses D1, D2, D3, and D4 with requirement of 5, 8, 7 and 14 units (in100s)
per week respectively. The transportation costs in birr per unit between factories to
warehouses are given in the table below:
(Production facilities) D1 D2 D3 D4
Solution
The First step is identifying the information needed to formulate the transportation table.
Generally, three information are needed to formulate the transportation table: demand
(also called destination or requirement), supply (also called source or capacity), and unit
transportation cost. These in formations are arranged as follows:
S1 7
S2 9
S3 18
warehouses demand(requirement
D1 5
D2 8
D3 7
D4 14
(Production facilities) D1 D2 D3 D4
To formulate the table, list the demand vertically (in column) and supply horizontally (in
rows). Thus, the production facilities (sources) are listed down the left side of the table,
and their respective supply’s quantities are listed down the right side of the table. The
destinations (warehouses) are listed across the top of the table and their respective
demands are listed across the bottom of the table. The unit shipping costs are shown in
the upper right hand corner of each sell, which represent the shipping route.
Transportation table
To D1 D2 D3 D4 Supply
From
S1 19 30 50 10 7
S2 70 30 40 60 9
S3 40 8 70 20 18
Demand 5 8 7 14 34
34
*Remark: when the total supply equals total demand, the problem is called balanced
transportation problem, otherwise an unbalanced transportation problem. The unbalanced
Remark: the quantity shipped cannot exceed the available supply in rows, and it should
not exceed the amount of demand (column total). Such solution is a feasible solution. The
starting point of the transportation method is a feasible solution. A feasible solution is
one in which assignments are made in such a way that all supply and demand
requirements are satisfied.
In this section, we shall discuss the following three different methods to obtain an initial
solution:
The initial solution obtained by any of the three methods must satisfy the following
conditions:
The solution must be feasible i.e. it must satisfy all the supply and demand
constraints
The number of occupied cell should equal one less than the sum of the number of
rows(m) and the number of columns(n)in a transportation table i.e. (m + n) – 1. For
example, in the case of a table with 3 rows and 3 columns, the number of occupied
cell should be: (3 +3) – 1 = 5 in order to be able to use the transportation algorithm.
Any solution that satisfies the above conditions is called non degenerate basic feasible
solution, otherwise, degenerate solution.
Dear student, now let as see the above three methods of finding initial feasible solution
one by one beginning with the NWCM.
The north –west corner method is a systematic approach for developing an initial basic
feasible solution. Its chief advantages are that it is simple to use and easy to understand.
Its chief drawback is that it does not take transportation cost in to account. Consequently,
such solution may require much additional effort to obtain the optimal solution.
The north-west corner method gets its name because the starting point for the allocation
process is the upper left hand (North- West) corner of transportation table.
Step1. Start with the sell at the upper left (north-west) corner of the transportation table
and allocate as many units as possible to that cell. This will be the smaller of the row
supply and the column demand. Adjust the row and column quantities to reflect the
allocation.
Step2. - If allocation made in step 1 is equal to the supply available in the first row, then
move vertically down to the cell in the second row and first column.
If allocation made in step1 is equal to the demand of the first destination (in first
column), then move horizontally to the cell in the first brow and second column and
apply step 1 again for next allocation.
Step3. Continue the procedure step by step till an allocation is made in the south east
corner cell of the transportation table.
Example2. Let as return to Example1 transportation table and see how to find the initial
feasible solution using NWCM.
D1 D2 D3 D4 Supply
S1 19 30 50 10 7
S2 70 30 40 60 9
S3 40 8 70 20 18
Demand 5 8 7 14 34
34
Required: consider the above table and then find the initial basic feasible solution (total
transportation cost) using NWCM.
Solution
To find the initial basic feasible solution using NWCM, track the following steps.
1. Beginning in cell (S1, D1) which is the north-west corner cell in the given
transportation table. The value for row S1 (i.e.7) and column D1 (i.e.5) are
compared. The smaller of the two, i.e. 5 is assigned as the first allocation;
otherwise it will violate the feasibility condition. This means that 5 unit of
commodities are to be transported from production facility 1 (S1) to warehouse 1
(D1). However, this allocation leaves a supply of 2 (7-5 =2) units of commodity
at S1.
2. Staying in row S1, move horizontally to cell (S1, D2) where supply is now 2 and
demand is 8. Allocate 2 units to this cell, exhausting the supply of row S1 and
leave 6 units of demand inD2 column.
3. Staying in column D2, move down to cell S2 D2, where supply is now 9 and
demand is 6. Allocate 6 units to this cell, exhausting the demand in column D2
and leaving 3 units of supply in row S2.
4. Staying in row S2, move to cell (S2, D3) and allocate 3 units to this cell,
exhausting that row’s supply and leaving 4 units in columnD3.
5. Staying in column D3, move down to cell (S3, D3) and allocate 4 units to this
cell, exhausting D3 columns and leaving 14 units of supply.
6. Staying in column S3 move to cell (S3, D4) and allocate 14 units to this cell,
exhausting both the rows and columns quantities.
D1 D2 D3 D4 Supply
S1 19 30 50 10 7
5 2
S2 70 30 40 60 9
6 3
S3 40 8 70 20 18
4 14
Demand 5 8 7 14 34
34
The total transportation cost of the initial solution derived by the NWCM is obtained by
multiplying the quantity in the occupied cells with the corresponding unit costs and
adding. Thus,
30 X 2 = 60
30 X 6 = 180
40 X 3 = 120
70 X 4 = 280
20 X 14 = 280
Birr 1015
In terms of minimizing total transportation cost, this solution may or may not be optimal.
We shall make that determination very shortly under ‘test for optimality’ section. At this
point, our target is simply computing basic initial feasible solution (i.e. total cost =1015).
As noted earlier, the main drawback of the north –west corner method is that it does not
consider cell (route) costs in making the allocation. Consequently, if this allocation is
optimal, that can be attributed to chance rather than the method used. To overcome this
drawback, we can use the LCM (which consider the unit transportation cost) to find the
initial feasible solution.
This approach, also called the minimum –cost method, uses lowest cell as the basis for
selecting route. Generally, this method takes in to account the minimum unit cost of
transportation for obtaining initial solution and can be summarized as follows:
Step1. Select the cell with the lowest unit cost in the entire transportation table and
allocate as much unit as possible to this cell and eliminate (line out) that row or column
in which either supply or demand is exhausted.
If both a row and column are satisfied simultaneously, only one may be crossed
out arbitrarily
In case the smallest unit cost cell is not unique, then select the cell where
maximum allocation can be made.
Step2. after adjusting the supply and demand for all uncrossed –out rows and columns
repeat the procedure with the next lowest unit cost among the remaining rows and
columns of the transportation table and allocate a quantity to this cell that is equal to the
lower of the available supply of the row and the demand for the column and eliminate
(cross out) that row or column in which either demand or supply is exhausted.
Step3. Repeat the procedure until the entire available supply at various sources and
demand at various destinations is satisfied. The solution obtained need to be non-
degenerate.
Example3
Taking data of the above example, the initial solution using LCM and total transportation
cost is calculated as follows:
The cell with lowest unit cost is 8 at cell (S3, D2). The maximum units which we
can allocate to this cell are selected from (18 and 8) whichever is smallest.
Therefore, at this cell we must allocate 8. This meets the complete demand of D2
and leaves 10 units with S3. Thus, crossed out column D2.
In the reduced table without column D2 (note: D2 cannot be used because it has
been crossed out), the next smallest unit transportation cost is 10 in cell (S1, D4).
The maximum unit which can be allocated to this cell can be 7(i.e. compare 14
and 7 and then assign the smallest unit). This exhausts the capacity of S1 and
leaves 7 units withD4 as unsatisfied demand. Thus, eliminate S1 row.
The next smallest cost, without considering S1, is 20 units at cell (S3, D4). The
maximum allocation in this cell could be 7 units. This satisfies the entire demand
of D4 and leaves 3 units with S3 as remaining supply. Thus line out D4.
The next smallest unit cost cell is not unique. That is, there are two cells (S2, D3)
and (S3, D1), having the same unit transportation cost of 40 birr. Allocate 7 units
in cell (S2, D3) first because it can accommodate more units as compared to cell
(S3, D1) and allocate 3 units (only supply left with S3) in cell (S3, D1). The
remaining demand of 2 units of D1 is fulfilled from S2. Since supply and demand
at each origin in destination is exhausted, the initial solution is arrived at. The
following table summarizes the discussion.
D1 D2 D3 D4 Supply
S1 19 30 50 10 7
S2 70 30 40 60 9
2 7
S3 40 8 70 20 18
3 8 7
Demand 5 8 7 14 34
34
70 X 2 = 140
40 X 7 = 280
40 X 3 = 120
8X8 = 64
20 X 7 = 140
Birr 814
allocation in certain cells with minimum unit transportation costs were missed. In this
method allocations are made so that the penalty cost is minimized. The advantage of this
method is that it gives an initial solution which is nearer to an optimal solution or is the
optimal solution itself. In VAM method, we look for cells having lowest cost and next
lowest cost in each row and column.
Step1. Calculate penalties for each row /column by taking the difference between the
smallest and next smallest unit transportation cost in the same row/column and note down
their difference. This difference indicates the penalty or extra cost which has to be paid if
one fails to allocate to the cell with the minimum unit transportation cost.
Step2. Select the row or column with the largest penalty (with maximum differences) and
allocate as much unit as possible in the cell having the least cost in the selected row or
column satisfying the rim condition. If there is a tie in the value of penalties, it can be
broken by selecting the cell where maximum allocation can be made.
Step3. Adjust the supply and demand and cross out the satisfied row or column and re-
draw the reduced matrix. If a row and column are satisfied simultaneously, only one of
them is crossed out and the remaining row or column is assigned a zero supply or
demand. Any row /column with zero supply or demand should not be used in computing
future penalties and hence should be eliminated.
Step4. Repeat step 1 to 3 until the entire available supply at various sources and demand
at various destinations are satisfied.
Example4
Taking data of the above example, the initial solution using VAM and total transportation
cost is calculated as follows:
As shown in the following table, differences (penalty costs) for each row and
column have been calculated. In the first round, the maximum penalty, 22 occur
in column D2. Thus, the cell (S3, D2) having the least transportation cost 8 is
chosen for allocation. The maximum possible allocation in this cell is 8 (compare
demand i.e. 8 and supply i.e.18 and assign the smallest unit) and it satisfies
demand in column D2. Adjust the supply of S3 from 18 to 10 (18-8=10) and
cross out column D2.
Calculate the new row and column penalties except for column D2 because its
demand has been satisfied. The penalties are:
- D1=21 (i.e. the smallest 19 subtracted from the next smallest 40 units)
- D3=10 (i.e. the smallest 40 subtracted from the next smallest 50 units)
- D4=10 (i.e. the smallest 10 subtracted from the next smallest 20 units)
- S1 = 9 (i.e. the smallest 10 subtracted from the next smallest 19 units)
- S2 =20 (i.e. the smallest 40 subtracted from the next smallest 60 units)
- S3 = 20(the smallest 20 subtracted from the next smallest 40 units)
The second round allocation is made in column D1 with target penalty 21 in the
same way as in the first round as shown in cell (S1, D1) of the table.
In the third round, the maximum penalty 50 occurs at S3. The maximum possible
allocation of 10 units is made in cell (S3, D4) having least transportation cost 20 as
shown in the table.
The process is continued with new allocations till a complete solution is obtained.
The initial solution using VAM is in the table.
S1 19 30 50 10 7 9 9 40 40
5 2
S2 70 30 40 60 9 10 20 20 20
7 2
S3 40 8 70 20 18 12 20 50 --
8 10
Demand 5 8 7 14 34
34
column 21 22 10 10
difference
21 -- 10 10
-- -- 10 10
-- -- 10 50
10 X 2 = 20
40 X 7 = 280
60 X 2 = 120
8X8 = 64
20 X 10 = 200
Birr 779
Once the initial solution is obtained, the next step is to check its optimality. An optimal
solution is one where there is no other set of transportation routes (allocations) that will
further reduce the total transportation cost.
The test of optimality for initial feasible solution involves a cost evaluations of
unoccupied or empty cell (i.e. routes to which no units have been allocated) to see if an
improved solution is possible. Thus, we have to evaluate each unoccupied cell in the
transportation table in terms of an opportunity of reducing total transportation cost.
In this chapter, we shall consider two methods for cell evaluation: the stepping stone
method and the MODI method.
Mr. A. Charnes developed stepping stone methods for solving transportation problems.
Subsequent improvements led to the developments of MODI methods in 1955. Stepping
stone methods, even though now superseded by MODI method, is taken up for discussion
in this section, for better understanding the logic behind its solutions. Now let as see first
the stepping stone method and later we will see the MODI in detail.
The stepping stone method involves tracing a series of closed paths in the transportation
table, using one such path for each empty cell. The name stepping stone relates to an
analogy of crossing a pond or scream by moving from stone to stone; in the case of
transportation problem the stones are the completed cell.
The steps for evaluating the optimality of initial feasible solution using stepping stone
method include the following:
Trace a ‘closed path’ that begins at the unoccupied cell, and moves alternately in
horizontal and vertical directions pivoting only on occupied cells and terminates
on the unoccupied cells.
A ‘+’ sign is assigned to the unoccupied cell, and succeeding corner points on the
path are alternately assigned a minus sign ‘–’ and ‘+’ sign.
The plus and minus sign indicates the necessary adjustments for satisfying the
row and column requirements. See the following example.
Example 5
Now let as take the initial feasible solutions found using LCM in Example 5 and test for
optimality using the stepping stone method.
A B C Supply
S1 5 7 (-)9 500
(+) 500
S2 6 2 4 600
500 100
600 200
1900
The index for each unoccupied cell can be determined by tracing a closed path.
Unoccupied cell Closed path with alternate ‘+’ and ‘–’ sign
S1, A +5 – 9 + 5 – 3 = -2
S1, B +7 – 9 +4 – 2 = 0
S2, A +6 – 4 +5 – 3 = 4
S3, B +9 – 5 +4– 2 = 6
Step2. If a better solution exists, determine which cell should be occupied (entering
cell). For minimization problems, a better solution exists if there are any negative
improvement indices. An optimal solution has been found when all improvement indices
are nonnegative (zero or positive). Once the optimal solution has been found, the
existence of one or more improvement indices at zero level indicates alternative optimal
solutions.
NB: If more than one negative index exists simultaneously, occupied (assign) to the cell
with largest negative number (i.e. largest number with minus sign). Based on the above
example, cell (S1, A) has negative indices (i.e.-2) indicating that total transportation cost
can be further reduced by 2Q (where Q is the number of units assigned to this cell).
Step3. Determine the departing cell and the number of units to assign the entering cell
The departing cell is identified by the smallest cell in a minus position on the closed path
for the entering cell. The entering cell is assigned the number of units previously given
the departing variable.
Based on the above example, the cells with minus signs are :( S1,C) & (S3 , A) with a
quantity of 500 unit and 600 units respectively.
The departing cell is the cell with the smallest quantity in a negative position.
Therefore, cell (S1, C) is the departing cell (since 500 is less than 600). Thus, the number
of unit assigned to the new entering cell (i.e.S1, A) is 500 units.
The new solution can be determined by adding 500 units to the cell with + sign (i.e. cell
S1, A and S3, C) in the closed path and by deducting 500 from the cell with minus sign
(i.e. S1, C, and S3, A).
(S1, A), (S2, B), (S2, C), (S3,A), (S3, C). See the following table.
A B C Supply
500
S2 6 2 4 600
500 100
100 700
1900
500 X 2 = 100
100 X 4 = 400
100 X 3 = 300
700 X5 = 3500
Birr 7700
To prove whether this figure is an optimal solution or not see the following. We stated
that total transportation cost can be reduced by 2Q where Q is the number of units
assigned to this cell. Now we identified that Q=500 units and the total transportation cost
reduction = 500X2= birr 1000
Unlike stepping stone method, the MODI algorithm eliminates the need to trace closed
paths for each unoccupied cell. The MODI method of evaluating a transportation
solution for optimality involves the use of Index numbers that are established for the row
and columns. These are based on the unit cost of the occupied cells. The index number
can be used to obtain the cell evaluation for empty cells without the use of stepping
stone-paths. There is one index number for each column and one for each row.
The steps to evaluate unoccupied cells using MODI method are as follows:
Step1. For an initial basic feasible solution with m + n -1 occupied cells, calculate Ui for
rows(i= 1,2,3……m) and Vj for columns(j=1,2,3……n).
-Start the process by assigning a value of zero for the index number of row1.
-then complete the calculation of Ui’s and Vj’s for other rows and columns using the
relations: cij = Ui + Vj where cij is the transportation cost of the cell.
Step2. For unoccupied cell calculate an opportunity cost by using the relation:
-if dij = 0, then current basic feasible solution will remain un affected but an alternative
solution exist
Step4. Construct a closed path for the unoccupied cell with largest negative opportunity
cost. start the closed path with the selected unoccupied cell and mark a plus sign (+) in
this cell, trace a path along the rows (or columns) to unoccupied cell, mark the corner
with minus sign (-) and continue down the column or row to unoccupied cell and mark
the corner with plus sign and minus sign alternatively. Close the path back to the selected
unoccupied cell.
Step5. Select the smallest quantity amongst the cells marked with minus sign on the
corners of closed path. Allocate this value to the selected unoccupied cell and add it to
other occupied cells marked with plus signs and subtract it from the occupied cells
marked with minus sign
Step6. Obtain a new improved solution by allocating units to the unoccupied cell
according to step 5 and calculate the new total transportation cost.
Step7. Test the revised solution further for optimality. The procedure terminates when all
dij > 0 for unoccupied cell.
Remark:
The path starts and ends at the selected unoccupied cell. It consists of successive
horizontal and vertical lines whose end points must be occupied cells, except for
the end point associated with entering unoccupied cell.
It is immaterial whether the loop (path) is traced in a clockwise or anti-clockwise
direction. However, for a given solution only one loop can be constructed for each
unoccupied cell.
Any two adjacent cells of the ordered set lie either in the same row or in the same
column.
Closed loops (path) may or may not be square in shape.
All cells that have a plus or minus sign, except the starting unoccupied cell, must
be occupied cell.
Ever loop has an even number of cells and at least four.
Example6
For better understanding turn back to Example 5 with the initial feasible solution found
using LCM method. The following table is the initial feasible solution found using LCM
(see Example5)
A B C Supply Ui
S1 5 7 9 500 U1=0
500
S2 6 2 4 600 U2=-5
500 100
S3 3 9 5 800 U3=-4
600 200
Solution
To test the initial feasible solution of LCM by using MODI methods followed the
following steps.
Step1. Calculate the index for Ui (row), Vj(column) for occupied cell.
S1, C 9 = U1 + V3 U1=0
9 = 0 + V3 V3=9
S2, B 2 = U2 +V2
2 = -5 + V2 V2=7
S2, C 4 = U2 +V3
4 = U2 +9 U2=-5
S3, A 3 = U3+V1
3 = -4+V1 V1=7
S3, C 5 = U3+V3
5 = U3 + 9 U3=-4
Step2 calculate the opportunity cost for unoccupied cell using the relation:
S1,A 5– 0– 7 -2
S1,B 7- 0– 7 0
S2,A 6 – (-5) - 7 4
S3,B 9- (-4)- 7 6
if one or more dij values < 0, then improved solution can be obtained by entering
unoccupied cell (i,j) having the largest negative value.
Based on the above data, cell (S1, A) has a negative index (dij).
This indicates that, if we occupied cell (S1, A), there will be birr two reduction of
transportation cost per unit.
Thus, total reduction of cost = 2Q (where Q is number of unit assigned to this cell)
Step4. Construct a closed path for the unoccupied cell with largest negative opportunity
cost starting the closed path with the selected unoccupied cell (i.e cell S1,A).
A B C Supply Ui
(+)500
S2 6 2 4 600 U2=-5
500 100
Step5. Select the smallest quantity amongst the cells marked with minus sign on the
corners of closed path and allocate this value to cell (S1, A). Then, add it to cell (S3, C)
and subtract it from cell (S3, A and S1, C).
Based on the above table, 500 and 600 are quantities found in cell marked with minus
sign. So, select 500 units and adjust the demand and supply as follows
Exercise
a) find the initial feasible solution using NWCM and VAM method
b) test the optimality of the initial feasible solution found in a above using stepping stone
methods and MODI method
In this section we will consider some special cases in solving transportation problems.
The most common special cases includes the following
For a feasible solution to exist, it is necessary that the total supply must equal total
demand. Unbalanced case may be of two types: supply> demand or demand>supply.
If total demand exceeds total supply, a dummy row (called a dummy supply
centre) can be added to the transportation table to account for the excess demand.
The unit transportation cost for the cells in the dummy row is set equal to zero.
Example 1
Three cement factories are located at A, B, C are to be supplied with carbonates stones
from mines X, Y, and Z. any mine can supply to any factory. The only condition is to
minimize the overall transportation cost. The total requirements (demand) and capacities
(supply of mines) are indicated in the table below, along with cost of transportation in
birr per ton.
i. ii. iii
Mines capacity per day factory requirements per cost of transpo
In Lorries loads in tones day in lorry loads in tons ration in birr
A B C
X 20 A 40 X 2 4 4
Y 50 B 25 Y 8 12 8
Z 30 C 85 Z 2 8 12
100 150
Required: formulate the transportation table and find the initial feasible solution using
NWCM.
Since demand and supply are not equal (i.e. demand exceeds supply) feasible solution
will not exist. This problem can be solved by a balanced situation and proceed as per
method explained earlier.
In order to balance supply and demand, create a dummy source of mine (row) having a
capacity by which supply is less than demand in the initial solution. In this case it is 50
(i.e.150-100=50). Dummy source is included in the table by adding one extra row.
Dummy source will have distribution cost =0. Then find the initial feasible solution
using NWCM. See table below.
To A B C Supply
From
X 2 4 4 20
20
Y 8 12 8 50
20 25 5
Z 2 8 12 30
30
Dummy 0 0 0 50
The initial feasible
solution using NWCM is:
50
Route unit
cost Demand 40 25 85 150 quantity shipped
cost
150
X, A 2
20 40
Y, A 8 20 160
Y, B 12 25 300
Y, C 8 5 40
Z, C 12 30 360
Dummy, C 0 50 0
Now proceed as above to till optimum solution is obtained using either stepping stone
method or MODI method.
If total supply exceeds total demand, an additional column (called dummy demand
centre) can be added to the transportation table to absorb the excess supply. The unit
transportation cost for these cells in this column is set equal to zero because these
represent product items that are not being made and not being sent.
Example 2
In the above example given in Example1 of this section, let as have the following supply
and demand. Unit cost is assumed to be the same.
i. ii. iii
A B C
X 20 A 40 X 2 4 4
Y 50 B 25 Y 8 12 8
Z 80 C 35 Z 2 8 12
150 100
Required: formulate the transportation table and find the initial feasible solution using
NWCM.
Sine total supply is greater than total demand by 50 units, add a dummy column with 0
unit cost and 50 units of demand. See the following table.
To A B C Dummy Supply
From
X 2 4 4 0 20
20
Y 8 12 8 0 50
20 25 5
Z 2 8 12 0 80
30 50
Demand 40 25 35 50 150
150
Z, C 12 30 360
Z, dummy 0 50 0
Total cost 900
Now proceed as above to till optimum solution is obtained using either stepping stone
method or MODI method.
b) Degeneracy
Degeneracy occurs when the number of occupied cell contains is less than (m+n-1). In
such cases the current solution cannot be improved because it is not possible to draw a
closed path for every occupied cell. Also the values dual variables Ui and Vj which are
used to test the optimality ca not be computed. Thus, we need to remove the degeneracy
to improve the given solution.
The degeneracy can arise either in developing the initial solution or at some intermediate
solution.
ε Δ
Add an artificial solution ‘ ’ (Greek letter epsilon) or (delta) to one or more of the
Δ
unoccupied cell to get ‘m + n -1’ number of occupied cells. The cell with is
considered as an occupied cell and is modified in the same manner as other occupied
cells.
ε Δ
‘ or is assumed to be a very small quantity close to zero
Δ
Generally, in a minimization transportation problem, it is better to allocate to
unoccupied cells that have lowest transportation costs. In maximization problems it
should be allocated to a cell that has a high pay off value.
Δ
The placement of Delta ( ) cannot be in just any unoccupied cells. It has to be
placed in either the column or the row associated with the cell which simultaneously
full filled the supply and demand constraints.
Δ
The quantity of is considered to be so small that if it is transferred to an occupied
cell it does not change the quantity of allocation. Hence it does not affect the total
Δ
transportation cost. Finally we can eliminate from the table.
Example. A manufacturer wants to ship 200 loads of his products from two factories
(factiry1 and factory2) to three warehouses (warehouse A, B, and C) as shown below.
To A B C Supply
From
X 5 10 8 75
75 Δ
Y 4 15 12 125
60 65
Demand 75 60 65 200
200
Solution
ii. To test the optimality of the initial solution found above, let as use stepping stone
method.
To A B C Supply
From
X (-)5 (+)10 8 75
75 Δ
60 65
Demand 75 60 65 200
200
Unoccupied cell Closed path with alternate ‘+’ and ‘–’ sign
X,C +8 – 12 + 15 – 10 = 1
Y, A +4 – 15 +10– 5 = -6
Since cell (Y,A) has a negative index, we must occupy this cell. This enables as to save
birr6 per unit of assignment at this cell. Therefore, total cost reduction=6Q.
To A B C Supply
From
X 5 10 8 75
75 Δ
Y 4 15 12 125
60 65
Demand 75 60 65 200
200
The quantities in the cell with minus sign are 75 and 60. Since the smallest is 60, assign
this units at cell (Y, A) and adjust the respective demand and supply amount.
To resolve degeneracy which occurs during optimality test , the quantity may be allocated
to one or more cells which have become unoccupied recently to have m+n-1 number of
occupied cells in the new solution.
Example
Goods have to be transported from sources S1, S2, and S3 to destination D1, D2, and D3.
The transportation cost per unit, capacity of the sources and requirements of the
destinations are given in the following table.
To D1 D2 D3 Supply
From
S1 8 5 6 120
S2 15 10 12 80
S3 3 9 10 80
Solution
Using the NWCM, the non degenerate initial feasible solution is given in the following
table:
To D1 D2 D3 Supply
From
S1 8 5 6 120
120
S2 15 10 12 80
30 50
S3 3 9 10 80
30 50
Total cost= (120 X8) + (15X30) + (10X50) + (9X30) + (10X50) = birr 2680
To test the optimality of the above solution found in NWCM using MODI method,
calculate Ui, Vj and dij. As usual as shown below.
To D1 D2 D3 Supply Ui
From
S1 8 5 6 120 U1=0
120
S2 15 (+)10 12 80 U2=7
30 (-) 50
30 50
S1, D1 8 = U1 + V1 U1=0
8 = 0 + V1 V1=8
S2, D1 15 = U2 +V1
15 = U2 + 8 U2=7
S2, D2 10 = U2 +V2
10 = 7 +V2 V2=3
S3, D2 9 = U3+V2
9 = U3+3 U3=6
S3, D3 10 = U3+V3
10 = 6+ V3 V3=4
Then calculate the opportunity cost for unoccupied cell using the relation:
S1,D2 5– 0– 3 2
S1,D3 6- 0– 4 2
S2, D3 12 – 7 -4 1
S3,D1 3- 6- 8 -11
Since the unoccupied cell (S3,D1) has the largest negative opportunity cost of -11,there
fore, cell (S3, D1)is entered in to the new solution mix. The closed path for (S3.D1) is
shown in the above table. The maximum allocation to this cell is 30. However, when this
amount is allocated to (S3,D1) both cells with minus sign i.e. (S2,D1) and(S3,D2)
become unoccupied because these two have same allocation(30 units each).thus, the
number of positive allocations became less than the required number(i.e. m+n-1)
To D1 D2 D3 Supply
From
S1 8 5 6 120
120
S2 15 10 12 80
80
S3 3 9 10 80
30 Δ 50
Δ
To remove the degeneracy a quantity is assigned to one of the cells that has become
Δ
occupied so that there are m+n-1 occupied cells. Assign to either (S2, D1) OR (S3,
D2) and proceed with the usual solution procedure.
Example
A company has received a contract to supply gravels from its three gravel pits (located in
towns W, X and Y) to three new construction projects (located in town A, B and C). The
transportation cost per unit, capacity of the sources (gravel pits) and requirements of the
destinations (construction projects) are given in the following table.
To A B C Dummy Supply
From
W 4 8 8 0 76
X 16 24 16 0 82
Y 8 16 24 0 77
(NB: Since demand is less than supply, add dummy column with a demand of 20 units.)
Required: determine the transportation Schedule so that cost is minimized.
Solution
First let as find the initial feasible solution using VAM
(Note that if the column difference and row differences are equal; select the cell in which
you can assign the maximum amount with the minimum transportation cost)
W 4 8 8 0 76 4 4 0 -
76
X 16 24 16 0 82 16 8 8 8
21 41 20
Y 8 16 24 0 77 8 8 8 8
72 5
Column 4 8 8 0
difference
8 8 10
4 8 8 10
First calculate the index for Ui (row), Vj(column) for occupied cell.
A B C dummy Supply Ui
NB: the
index of W 4 8 8 0 76 U1=0 U1=0
The 76 initial
feasible
X 16(+) 24 (-) 16 0 82 21 U2=16
21 41 20
72 5
8 = 8 +V1 V1 = 0
Y, B 16 = U3 + V2
16 = U3 + 8 U3= 8
Then, calculate the opportunity cost for unoccupied cell using the relation:
The opportunity costs in all unoccupied cells are positive except for cell (X,A) which has
a zero opportunity cost. This means, if (X, A) is occupied, no change in the transportation
cost would occur. To determine the alternative solution, form a closed path for cell(X,A)
as shown above. Quantity of 21 and 72 are in the cell with minus sign. So, the maximum
quantity which can be allocated to cell (X, A) is 21.after this change, the new solution is
shown in the following table.
A B C Dummy Supply Ui
W 4 8 8 0 76 U1=0
76
X 16 24 16 0 82 21 U2=16
21 41 20
Y 8 16 24 0 77 5 U3=8
51 26
Since all dij values are positive or zero, this solution is optimal with a minimum total
transportation cost of birr 2424which is the same as in the previous solution.
d) Unacceptable route
In general, transportation model is used for cost minimization problems. However, some
transportation-type problems concern profits or revenues rather than costs. In this case
the objective is to maximize rather than to minimize.
The algorithm for solving maximization problem is the same as that for the minimization
problems except that it adds one additional step at the start:
First identify the cell with the largest profit and subtract all the other cell profits
from that value. These values reflect the opportunity costs that would be incurred
by using routes with unit profits that are less than the largest unit profit.
Then replace the original unit profits with these opportunity costs and solve in the
usual way for the minimum opportunity cost solution. This will be identical to
maximizing the total profit.
When the optimal distribution plan has been identified, use the original cell
values (i.e. profits) to compute the total profit for that plan.
Example
A company has three manufacturing plants and three warehouses. These data along with
the unit profit are given below.
NB: the steps for developing an initial feasible solution, evaluation of optimality and
reallocation are identical to those used for cost minimization.
Solution
To D1 D2 D3 Supply
From
S1 4 2 8 100
S2 5 1 9 200
S3 7 6 3 200
ii. To find the initial feasible solution, first identify the largest unit profit and the subtract
all cells from that value. In this case the largest unit profit is 15.
Then subtract each cell from 15 to get the opportunity cost as shown below.
D1 D2 D3 D1 D2 D3
S1 9- 4=5 9- 2=7 9- 8 =1 S1 5 7 1
S2 9-5=4 9- 1= 8 9- 9=0 S2 4 8 0
S3 9-7 =2 9- 6 =3 9- 3=6 S3 2 3 6
To D1 D2 D3 Supply
From
S1 5 7 1 100
S2 4 8 0 200
S3 2 3 6 200
Find the initial feasible solution using VAM, LCM or NWCM method. Let as use
NWCM.
To D1 D2 D3 Supply
From
S1 5 7 1 100
50 50
S2 4 8 (+)0 200
100(-) 100
200
500
As usual let as use the stepping stone method to test for optimality.
The index for each unoccupied cell can be determined by tracing a closed path.
Unoccupied cell Closed path with alternate ‘+’ and ‘–’ sign
S1, D3 +1 – 0 + 8 – 7 = 2
S1, D1 +4 – 8 +7 – 5 = -2
S3, D1 +2 – 6 +0 – 8+7-5 = 0
The index for cell (S1, D1) and (S3, D2) has a negative index which indicates that further
cost reduction (in this case profit improvement) is possible. Since cell (S3, D2) has the
largest number with mines sign (-11), occupied this cell. Thus the total reduction of cost
=11Q (where Q is the amount allocated to cell (S3, D2)
For the optimum assignment see the clothed path in the above table. The quantity 200
and 100 are assigned in the cell with minus sign. Therefore, select the smallest (i.e.100)
to be assigned at cell (S3, D2). The revised allocation is shown below.
To D1 D2 D3 Supply
From
S1 5 7 1 100
50 50
S2 4 8 0 200
200
The optimum total
S3 2 3 6 200
transportation
cost can 100 100 be calculated as
follows:
Demand 50 150 300 500
Route unit cost
quantity shipped cost
S1, D1 5 50 250
S1, D2 7 50 350
S2, D 3 9 200 0
We stated that the total transportation cost reduction is 11Q. Since Q=100, the total cost
reduction is birr 1100.
Using this optimal allocation the total profit can be computed by using the original cell
values (unit profits).
The assignment problem is a special case of Transportation problem in which the number
of sources and destinations are the same, and the objective is to assign the given
job(tasks) to most appropriate machine(person)so as to optimize the objective function
like minimizing cost or maximizing profit.
The problem of assignment arise because available resources such as men, machine etc
have varying degree of efficiency for performing different activities. Therefore cost,
profit, or time of performing the different activities is different. Thus, the problem is:
‘How the assignment is made so as to optimize the given objective’.
Assumptions
number of jobs is equal to number of machines or persons
each man or machine is loaded with one and only one job
each man or machine is independently capable of handling any of the job being presented
Loading criteria must be clearly stated such as “minimizing operating time”, or
“maximizing profit”, or “minimizing production costs” or “minimizing production cycle
time” etc.
the data matrix of the assignment model is similar to transportation model except that the
supply of each resource and demand at each of the destinations is taken to be one i.e.
assignment is made on the bases of 1:1
If the numbers of rows are not equal to the number of columns and vise-versa, a dummy
row or dummy column must be added. The assignment costs for dummy cells are always
zero.
Step2. Find the opportunity cost table using row reduction and column reduction.
The technique for determining whether an optimal assignment is possible at this stage
consists of drawing straight lines (vertically and horizontally) through the opportunity
cost table in such a way as to minimize the number of lines necessary to cover all zero
entries (i.e. draw a line through a column or row which has the largest number of zeros)
If the number of lines equals the number of either the number of rows or columns, an
optimal assignment is made. If the number of lines less than the number of rows or
column, an optimal assignment cannot be determined, and the opportunity cost table must
be revised.
From among the cells not covered by any line, chose the smallest elements. Call
this value ‘k’.
Subtract k from every element in the cell not covered by a line
Add k to every element in the cell covered by the two lines i.e. intersection of two
lines.
Element in cells covered by one line remain unchanged.
Step5. You should repeat the process of finding the minimum number of covering lines.
Repeat steps 3 and 4 until number of lines equal number of rows or column.
If number of line is equal with either number of columns or rows, the optimum
assignment is possible. Here, the question is “how to assign?”
To avoid confusion in making the assignments, begin by identifying a row or column that
has only one zero, and make an assignment to that cell (box that cell). Since no other
assignments can be made in this row or column (row or column which contain the
assigned cell), we cross them off. With crossed off column and row, look for a row or
column in which there is only one zero and assigned it. Continue until all zeros are either
crossed or assigned.
Example1
The district court assigns four judges (A, B, C, and D) to four court dockets (i, ii, iii, and
iv). Estimated time (in days) each court will take to clear each docket are given below.
The objective of the assignment is to minimize the total time required to complete all of
the cases on the four dockets. Based on the given data, determine the number of days
each judge would take to clear each docket.
Docket
judge i ii iii iv
A 14 13 17 14
B 16 15 16 15
C 18 14 20 17
D 20 13 15 18
Solution
The number of rows = number of column=4. so no need of adding dummy row or column
Row reduction: the minimum time (day) in row A, B, C, and D are: 13, 15 14, and 13
respectively. Subtract these elements in their respective row. This shows the additional
time that would be incurred if the lowest time assignment is not made either in terms of
rows or columns. For example, judge A can handle docket ii with the lowest processing
time of 13 days. Therefore, if judge A were assigned docket i, the additional time would
be, 14-13=1 day. Similarly, if judge A were assigned docket iii and iv, the additional time
would be 4(17-13=4) and 1(14-13=1). We can perform similar calculation for the other
rows.
Original estimated days to clear dockets estimated days after row reduction
Docket Docket
judge i ii iii iv judge i ii iii iv
A 1 0 4 1 A 0 0 3 1
B 1 0 1 0 B 0 0 0 0
C 4 0 6 3 C 3 0 5 3
D 7 0 2 5 D 6 0 1 5
Column
Minimum 1 0 1 0
docket
i ii iii iv
A 0 0 3 1
B 0 0 0 0
C 3 0 5 3
D 6 0 1 5
Now compare the number of line with the number of columns or rows. Since the
minimum number of covering lines (i.e. 3) is not equals with the number of rows or
columns (i.e.4), at this stage optimum assignment is not possible.
(NB: note that there is another way of covering the first column by vertical line instead
of the line drawn through the first row. This will not change the final result. Here the
point is that, only three lines drawn would be needed.)
Step4. Develop the new revised the opportunity cost table
Since number of line is not equal with number of rows/columns, we must revise the
opportunity time. See the following:
From among the cells not covered by any lines (i.e.3, 6, 5, 13, & 5), chose the
smallest element (i.e. 1). So k =1
Subtract 1 from every element in the cell not covered by a line
Add k to every element in the cell covered by the intersection of two lines i.e.
cell(A, ii)=0+1 and cell (B, ii)=0+1
Element in cells covered by one line remain unchanged: i.e. cell (A, i) ,(A, iii),
(A,vi), (B,i) ,(B, iii) (B,iv), (C, ii), and (D, ii)
The revised opportunity cost table is shown below
docket
i ii iii vi
A 0 1 3 1
B 0 1 0 0
C 2 0 4 2
D 5 0 0 4
Step5. You should repeat the process of finding the minimum number of covering lines.
Draw the line and compare it with the number of rows or columns.
Now the number of line is 4 which are equal with the number of row or columns. There
fore the next step is making the assignment.
The procedure for identifying the assignment cell is to select the row or column in which
there is only one zero. In this example, the candidates are: row A, row C & column iv.
Then cell at row B column iv can be assigned, and cross row B and column iv. Next
assign cell row D column iii. Finally assign row C column ii.
docket
i ii iii vi
0
A 1 3 1
0
0
B 0 1 0
0
C 2 4 2
D 5 0 4
The total days associated with this solution is obtained by adding original day in the
occupied cells as shown below.
Judge A to docket i 14
Judge B to docket iv 15
Judge C to docket ii 14
Example 2
A computer centre has three expert programmers. The centre wants three application
programs to be developed. The head of the computer centre, after studying carefully the
programs to be developed, estimate the computer time in minutes required by the expert
for the application programs as follows:
Programmers
A B C
1 120 100 80
Programs 2 80 90 110
Required: assign the programmers to the programs in such a way that the total computer
time is minimum.
Solution
Step1. Since numbers of rows are equal with number of column, dummy is not needed.
A B C minimum
A B C
Programs 2 80 90 110 80 2 0 10 30
Programmers’ programmers
A B C A B C
Programs1 40 20 0 Programs 1 40 10 0
2 0 10 30 2 0 0 30
3 0 30 10 3 0 20 10
Column
Minimum 0 10 0
Programmers
A B C
Programs 1 40 10 0
2 0 0 30
3 0 20 10
At this stage optimal assignment is possible because the numbers of lines are equal with
the number of row (columns)
Programmers
A B C
0
Programs 1 40 10
0
2 0 30
0
3 20 10
Program 1 to programmer C 80
Program 2 to programmer B 90
Total 280
While making an assignment in the reduced assignment matrix, it is possible to have two
or more ways to strike off a certain number of zeros. Such a situation indicates multiple
optimal solutions with the same optimal value of objective function. In such cases the
more suitable situation may be considered by the decision maker.
There may arise situations when the assignment problem calls for maximization of profit,
revenue, sales etc. as the objective function. In this case, one extra step must be added to
start of the process: identify the largest value (pay off element) in the assignment table
and then subtract all the elements of the table from the largest element.
Example
A marketing manager has three sales men and three sales districts. Considering the
capabilities of the sales men and the nature of districts, the marketing manager estimates
that cell per month (in ‘000’birrs) for each sales man in each district would be as follows:
Districts
1 2 3
Sales men A 14 22 30
B 20 18 40
C 11 12 50
Required: find the assignment of salesmen to districts that will result maximum sales.
Solution
The new cost data (which is equivalent with maximization of the original data) is given
below:
Districts
1 2 3
Sales men A 36 28 20
B 30 32 10
C 39 38 0
* First reduce the row *Table after row r deduction
District District
1 2 3 1 2 3
sale A 16 8 0 A 0 0 0
men B 20 22 0 B 4 14 0
C 39 38 0 C 23 30 0
Column
Minimum 16 8 0
Then as usual, cover all zeros (from the table found after reducing the table) by lines as
shown below.
District
1 2 3
sale A 0 0 0
men B 4 14 0
C 23 30 0
Since the number of table is less than the number of rows or columns, develop a revised
opportunity table. To develop the revised table: add 4 to cell (A,3) and deduct 4 from un
covered elements i.e. cell (B,1), (B,2) (C,1) and (C, 2) and cover zeros by line.
District
1 2 3 since the number of lines (3) is equal with the
Sale A0 0 4 number of column or rows (3) optimum
Men B0 10 0 assignment can be possible at this stage.
C 19 26 0
District
1 2 3 * there is only one zero in row C. Thus assign this
Sale A 0 00 4 cell and cross row C and column3.
Men 0
B 10 00 * now row B has only one zero. Assign this cell and
C 19 26 0 cross off row B and column 1.
* Finally assign row cell (A, 2)
Optimum assignments sales (in’000’)
Salesman A to district 2 22
Salesman B to district 1 20
Salesman C to district 3 50
Total 280 Thus, total sale = 280,000 birr.
Any assignment problem can be solved if and only if the opportunity cost matrix is a
square matrix i.e. the number of column must be equal with the number of rows. That is
why solving the assignment problem begins with checking whether the matrix is square
or ton.
If the given matrix is not a square matrix, the assignment problem is called unbalanced
problem. In such cases, a dummy row(s) or column(s) are added in the matrix (with zero
as a cost elements) to make it a square matrix. After making the given assignment matrix
a square matrix, similar procedures can be used to solve the problem. See the Example
given under next section (section 4)
3. Restriction on assignment
Example
In the modification of a plant lay out of a factory 4 new machines M1; M2, M3 and M4
are to be installed in a machine shop. There are 5 vacant places A, B, C, D and E
available. Because of limited space, machine M2 cannot be placed at C and M3 cannot be
placed at A. The cost of locating a machine at a place (in ‘00’s of birr) is as follows:
Location
A B C D E
M1
9 11 15 10 11
Machine M2 12 9 -- 10 9
M3 -- 11 14 11 7
14 8 12 7 8
M4
Required: find the optimal assignment schedule and the minimum cost.
Solution
Since the number of rows is not equal with number of column, add one dummy row to
get a square cost matrix.
* Second assign a high cost, denoted by M to the pair ( M2, C) and ( M3,A). The cost
matrix so obtained is shown below.
Location
A B C D E
M1
9 11 15 10 11
M2 12 9 M 10 9
M3 M 11 14 11 7
14 8 12 7 8
M4
0 0 0 0 0
M5
Now apply the assignment procedures to find the optimum assignment i.e. first reduce the
row and then the column. Next draw the line to cover zeros. If all zeros are covered by a
line and the numbers of lines are equal with the number of rows or columns, assign the
cell. See the following.
Location location
A B C D E A B C D E
M1 9 11 15 10 11 M1 0 2 6 1 2
3 0 M 1 0
12 9 M 10 9
M2 M2 M 4 7 4 0
7 1 5 0 1
M 11 14 11 7 0 0 0 0 0
M3 M3
14 8 12 7 8
M4 M4
0 0 0 0 0
M5 M5
Since the number of lines (5) is equal with the number of rows or columns (5), the
optimal solution can be obtained at this stage.
Location
A B C D E
0 2 6 1 2
M1 The optimum assignment and the total mini
3 M 1 0
0
M2 mum costs are:
M 4 7 4
M4 0 0 0 00 MachineM1 to location A 9
M5 MachineM2 to location B 9
0
MachineM3 to location E 7
MachineM4 to location D 7
MachineM5 to location C 0
We examine all the rows starting from row ‘A’ one-by-one until a row containing only a
single zero element is located (assigned). Here rows ‘C’ and ‘D’ have only one zero
element in the cell (C ,ii) and (D, ii). Assignment is made first at cell (C, ii) and these
cells. All zeros in the assigned columns are now crossed off as shown below
We now examine each column starting from column1 (start from the first column).there
is one zero in column iii and column iv in the cell (A, iii) and (A, iv). Assignment is
made in these cells. Thus, cell
If the numbers of rows are not equal to the number of columns and vise-versa, a dummy
row or dummy column must be added. The assignment costs for dummy cells are always
zero.
Step2. Find the opportunity cost table using row reduction and column reduction.
Examine the rows successively until a row with exactly one unmarked zero is
obtained. Make an assignment to this single zero by making square () around
it.
For each zero value that becomes assigned, eliminate (strike off) all other zeros
in the same row and /or column
Repeat these two steps for each column and also with exactly single zero value
cell that has not been assigned
If a row and/or column has two or more unmarked zeros and one cannot be
chosen by inspection, then choose the assigned zero cell arbitrarily
Continues these process until all zeros in rows/columns are either enclosed
(assigned) or struck off with ‘X’ mark.
If the number of assigned cells is equal to the number of rows /columns, then it is an
optimal solution. The total cost associated with this solution is obtained by adding
original cost figures in the occupied cells. If a zero cell was chosen arbitrarily in step 3,
there exists an alternative optimal solution. But if no optimal solution is found, then go to
step 5.
Draw a set of horizontal and vertical lines to cover all the zeros in the revised cost table
obtained from step 3 by using the following procurers:
For each row in which no assignment was made, mark a tick ()
Examine the marked rows. If any zero cell occurs in those rows, mark a to
the respective columns that contain those zeros.
Examine the marked columns. If any assigned zero occurs in those columns, tick
the respective rows that contain those assigned zeros.
If the number of lines drawn (or total assignment) is equal to the number of rows (or
rows) the current solution is the optimal solution, otherwise, go to step 6.
From among the cells not covered by any line, chose the smallest elements. Call
this value ‘k’.
Subtract k from every element in the cell not covered by a line
Add k to every element in the cell covered by the two lines i.e. intersection of
two lines.
Element in cells covered by one line remain unchanged.
Example1
The district court assigns four judges (A, B, C, and D) to four court dockets (i, ii, iii, and
iv). Estimated time (in days) each court will take to clear each docket are given below.
The objective of the assignment is to minimize the total time required to complete all of
the cases on the four dockets. Based on the given data, determine the number of days
each judge would take to clear each docket.
Docket
judge i ii iii iv
A 14 13 17 14
B 16 15 16 15
C 18 14 20 17
D 20 13 15 18
Solution
Row reduction: the minimum time (day) in row A, B, C, and D are: 13, 15 14, and 13
respectively. Subtract these elements in their respective row. This shows the additional
time that would be incurred if the lowest time assignment is not made either in terms of
rows or columns. For example, judge A can handle docket ii with the lowest processing
time of 13 days. Therefore, if judge A were assigned docket i, the additional time would
be, 14-13=1 day. Similarly, if judge A were assigned docket iii and iv, the additional time
would be 4(17-13=4) and 1(14-13=1). We can perform similar calculation for the other
rows.
Original estimated days to clear dockets estimated days after row reduction
A 14 13 17 14 13 A 1 0 4 1
B 16 15 16 15 15 B 1 0 1 0
C 18 14 20 17 14 C 4 0 6 3
D 20 13 15 18 13 D 7 0 2 5
Column reduction: the procedure of column reduction is similar with row reduction.
Docket Docket
A 1 0 4 1 A 0 0 3 1
B 1 0 1 0 B 0 0 0 0
C 4 0 6 3 C 3 0 5 3
D 7 0 2 5 D 6 0 1 5
Column
Minimum 1 0 1 0
We examine all the rows starting from row ‘A’ one-by-one until a row containing only a
single zero element is located (assigned). Here rows ‘C’ and ‘D’ have only one zero
element in the cell (C, ii) and (D, ii). Assignment is made first at cell (C, ii) and, all zeros
in the assigned column and now must be crossed off.
We now examine each column starting from column1 (start from the first column and
then to the next). In this example, there is one zero in column iii and column iv in the cell
(A, iii) and (A, iv). Assignment is made in these columns (first column iii as per their
order and then the next column). Then, all zeros in the assigned cells column and row
must be crossed off. See the following.
Docket
i ii iii vi
A 0 0 3 1
B 0 0 0 0
C 3 0 5 3
D 6 0 1 5
Chapter Summery
In order to find initial solution of a transportation methods the following methods used:
The north –west corner method is a systematic approach for developing an initial basic
feasible solution. Its chief advantages are that it is simple to use and easy to understand.
List cost approach, also called the minimum –cost method, uses lowest cell as the basis
for selecting route.
The assignment problem is a special case of Transportation problem in which the number
of sources and destinations are the same, and the objective is to assign the given job
(tasks) to most appropriate machine (person) so as to optimize the objective function like
minimizing cost or maximizing profit.
Review exercise
1. Suppose that a firm has three factories / sources of supply /& four warehouses/point of
demand/ .The firm's production capacity at the three factories, the demand for the four
destination centers located at various regions & the cost of shipping each unit from the
factories to the warehouses through each route is given as follows:
Destinations
Factory
W1 W2 W3 W4 Capacity
F1 3 2 7 6 5000
F2 7 5 2 3 6000
F3 2 5 4 5 2500
Required:
a. Develop an initial feasible solution using NWCM & Compute the total cost
b. Develop an initial feasible solution using least-cost method & compute the total cost.
[
2. Three garment plants are available for monthly education of four styles of men's shirts. The
capacities of the three plants are 45,000, 93,000 and 60,000 shirts. The number of shirts required
in style "a" through "d" is 28,000, 65,000, 35,000 & 70,000, respectively. The profits, in $ per
shirt, at each plant for each style are shown below.
STYLE a b c D
PLANT
1 8 12 -2 6
2 13 4 3 10
3 0 7 11 8
How many shirts of each type to produce in each plant so that profit is maximized?
3. A car rental company has one car at each of five depots a, b, c, d and e. A customer in
each of the five towns A, B, C, D and E requires a car. The distance in (in kilometers)
between the depots and towns where the customers are, is given in the following distance
matrix:
Depots
a B c d e
D 50 50 90 80 110
E 55 35 70 80 105
How should the cars be assigned to the customers so as to minimize the distance traveled?
Chapter 4
Decision theory
Learning objectives:
Describe and give example of decision making under the three conditions
(certainty, risk and complete uncertainty)
Construct and understand the pay off table
Make conclusion by using maximin, maximax, minimax regret, inseaficent reason
and expected value criteria.
Find out the expected value of perfect information
4.1. Introduction
Making appropriate decision is the most vital aspects in management. Infect certain
authors have defined management as decision making. Decision making is an action.
Every action has a reaction. Some decision initiates a set of activities; some put an end a
certain activities. Some decision can be withdrawn without any consequential actions and
losses, but majority of the decisions are not. Therefore, the success or failure of an
individual or organization experiences, depends to large extent on the ability of making
appropriate decisions
Decision theory provides an analytical and systematic approach to the study of decision
making. Decision models useful in helping decision makers make the best possible
decisions are classified according to the degree of certainty. The scale of certainty can
range from complete certainty to complete uncertainty. The region which falls between
these two extreme points corresponds to the decision making under risk (probabilistic
problems).
Irrespective of the types of decision model, there are certain essential terms which are
common in decision making. Let as defines theses terms:
Decision alternatives: these are act, action, option, or strategies available for decision
maker before any decision is made. Decision alternatives are under the control and
known to decision maker.
State of nature: a possible future condition (consequence or event) resulting from the
choices of a decision alternatives depends up on certain factors beyond the control of
decision maker. These factors are depends on the environment. Decision maker can
at the most attach subjective probability of occurrence of such events.
Payoff: a numerical value resulting from each possible combination of alternatives
and state of nature is payoff. The payoff might be profits, revenues costs, or other
measures of values. The payoff values are always conditional values because of
unknown state of nature.
Degree of certainty. There can be different degree of certainty. One extreme is
complete certainty and the other is complete uncertainty. The latter exists when the
likely hood of the various state of are unknown. Between theses two extremes is risk,
a term that implies that probabilities are known for the state of nature. We will see
these in detail in the next section.
Step2. List all possible future events (state of nature) which are beyond the control of
the decision maker
Step3. Identify all courses of actions (alternatives) which are under the control of the
decision maker
Step4. Express the payoff resulting from each pair of course of action and state of nature.
Normally payoffs are expressed in a monetary value.
Step5. Apply appropriate decision theory model to select the best courses of action.
A tabular arrangement conditional outcome (payoff) values is known as payoff matrix (or
decision table or payoff table). It includes a list of the alternatives, the possible state of
nature, and the payoff associated with each of the alternatives-statue of nature
combinations. Conventionally, rows represent the state of nature and columns represent
decision alternatives to be considered. Please not that pay off can be profits, revenues
costs, or other measures of values.
Example 1
A firm manufactures three types of products (A, B, C). The fixed cost and the variable
cost of producing each product are given below.
Product A 25000 12
Product B 35,000 9
Product C 53,000 7
In this case the payoff is expressed in profit. So first calculate the profit for each
alternative. There are three products and three demands. Therefore, there are 9(3x3=9)
pairings or pay off cells.
Let D1, D2, and D3 be poor, moderate, and high demand respectively, then payoff is
given by:
The calculation for payoff (in ‘000’s) for each pair of alternative demand (course of
action) and the types of product (state of nature) are shown below:
Payoff =( p q) – [(v q)+ FC ]
D 1 , A=(3000 X 25)−¿
Thus, the pay off values in the payoff table (in ‘000’s) is shown below.
(State of nature) D1 D2 D3
Product A 14 66 118
Product B 13 77 141
Product C 1 73 145
Knowledge of the likely hood of each of the state of nature can play an important role in
selecting a course of action. Generally there are three decision making environments:
certainty, uncertainty and risk.
Decision under certainty is the case where the decision maker has a perfect knowledge
(information) about the state of nature. This situation does not need any techniques or
managerial tool for the manager to take his decision. In this case, the decision maker will
simply select an alternative that yields the largest return (payoff) for the known future
(state of nature). However, in practice, decision under certainty is very rare. So this
situation is not of our interest and we left it for the reader.
This is a situation in which decision maker has neither previous knowledge of the
outcomes of his decisions, nor is in a position to attribute a probability of occurrence of
such outcomes. Still various options are opened to you to make decision in the absence of
perfect knowledge which depends on a variety of criteria. In this section we will
introduce you only four approaches for decision making under uncertainty.
In this section, the decision maker ensures that he should not miss the opportunity to
achieve the largest possible profit (maxmaxi) or lowest possible cost (minimini). Thus, he
selects the alternatives that represent the maximum of the maxima (or minimum of the
minima) payoffs (outcomes).
Locate the maximum (or minimum) payoff values corresponding to each alternative.
Then,
Select an alternative with best anticipated payoff value(maximum for profit and
minimum for cost)
Example 2
To illustrate the maximaxi criteria, let as refer the example given under Example1.
Solution
In order to determine the strategy (alternative action) first determine the maximum value
for each column. Then, the column with the largest value is considered as the best
strategy.
(State of nature) D1 D2 D3
Product A 14 66 118
Product B 13 77 141
Product C 1 73 145
Maximaxi
The maximum of column maxima is 145. Hence the manufacturing company should
adopt alternative D3 (sale in the high demand).
In this criterion, the decision makers ensures that he would earn no less or (pay no more)
tan some specified amount. Thus he selects the alternative that represents the maximum
of the minima (in case of profit) and the minimum of the maximum (in case of loss). The
application is simple. First select the minimum from each column and select the largest
values from the minimum column rows.
Example3
(State of nature) D1 D2 D3
Product A 14 66 118
Product B 13 77 141
Product C 1 73 145
maxmin
Since the probabilities of state of nature are not known, it is assumed that all state of
nature will occur with equal probability, i.e. each state of nature is assigned an equal
probability.
compute the expected(average) value payoff for each alternative by adding all the
pay offs and dividing by the number of possible states of nature
Select the best expected payoff value (maximum for profit or minimum for cost).
a) minimax regret
This criteria also called opportunity loss decision criteria or savage criterion because
decision makers feels regret after adopting a wrong course of action (alternative)
resulting an opportunity loss of payoff. Thus, the decision maker always intends to
minimize this regret. The work method is as follows:
from the given pay off matrix, develop an opportunity loss(or regret)matrix as
follows:
For each course of action(alternative) identify the worst (maximum regret value)
Select the course of action with the smallest anticipated opportunity loss value. For
better understanding of the concept let as see the following Example
Example
Given the following payoff tables, determine which alternatives would be chosen using
each of these decision criteria.
State of alternative
nature a b c d
1 12 17 22 14
2 18 10 16 14
3 15 14 10 14
Solution
i. The maximaxi approach to select the best alternative using maximaxi approach; first
select the maximum values from each column. Then, the best strategy is the
minimum column row with the highest values. See the following computation.
Strategy
State of (alternative)
nature a b c d
1 12 17 22 14
2 18 10 16 14
3 15 14 10 14
Maximum
mximaxi
The maximum value from the maximum column row is 22. Thus select strategy C
To select the best alternative using maximini approach, first select the minimum values
from each column. Then, the best strategy is the column with the highest values. See the
following computation.
Strategy
State of (alternative)
nature a b c d
1 12 17 22 14
2 18 10 16 14
3 15 14 10 14
Column 12 10 10 14----maximin column row
minimum
Mximin
The maximum value from the maximum column row is 14. Thus select strategy d.
In order to the maximin regret decision, we must first obtain the opportunity loss or regret
table. To get this table, first identify the best payoff in each row (i.e. 22, 18, and 15 for
row 1, 2, and 3 respectively). Then subtract every payoff in each row from the best pay
off rows. See the following.
nature a b c d
Strategy
State of (alternative)
nature a b c d
1 10 5 0 8
2 0 8 2 4
3 0 1 5 1
Maximum)
Minimum regret
The minimum value from the worst regret row is 4. Thus select strategy C.
To select the best alternative using Laplace approach, first we must determine the average
payoff for each alternative. Then select the minimum values from each column. Then, the
best strategy is the column with the highest values. To illustrate the laplace criterion
application, let as refer the above Example once again.
Strategy
State of (alternative)
nature a b c d
1 12 17 22 14
2 18 10 16 14
3 15 14 10 14
Best payoff
Since the largest payoff is obtained from alternative C, select this strategy.
In this case also the decision maker has to make decision when outcomes are not certain.
However, unlike the previous case, he has sufficient information to assign probability
values to the likely occurrence of each of the states. Knowing the probability distribution
of the state of nature, the best decision is to select that course of action (decision
alternatives) which has the largest expected payoff value.
The most widely used criterion for evaluating various courses of action (alternatives)
under risk includes: expected monetary value (EMV), expected opportunity loss (EOL),
expected value of perfect information (EVPI).
Expected monetary value is the sum of the payoffs for each course of action multiplied
by the probabilities associated with each state of nature.
Construct a payoff matrix listing all possible courses of action and state of nature.
Enter the conditional payoff values associated with each possible combination of
courses of action and state of nature along with the probabilities of the occurrences
of each state of nature.
Calculate the EMV for each course of action by multiplying the conditional payoffs
by the associated probabilities and add these weighted values for each course of
action.
Select the course of action that yields the optimal values.
Example
Mr. Abebe quite often flies from Adis Ababa to Debrezeyt. He can use the air port bus
but it costs him birr 25 and there is a 0.8 chance that he will miss the bus. He can also
stay in hotel but it costs him birr 270 with a 0.96cahnces of being on time for the flight.
Still he can use Taxi which costs him birr 350 with a 0.99 chance of being on time for
flight. If Mr. A bebe catches the plane, he will conclude a business transaction which will
produce a profit of birr 10,000 otherwise, he will lose it.
Required: using EMV method, which mode of transportation should Mr.Abebe use?
Solution
Computations of EMV for the three courses of action are shown below.
Bus
Stay in hotel
Taxi
Thus, comparing the EMV of the three courses of action, the alternative with the highest
EMV is the third alternative i.e. using Taxi. So Mr. Abebe must use taxi to maximize his
expected monetary value.
Example
A company manufactures goods for a market in which the technology of the product is
changed very rapidly. The R&D department has produced a new product which appears
to have potential for commercial exploitation. A further birr 60,000 is required for
development testing. The company has 100 customers and each customer might purchase
at the most one unit of the product. Market research suggest that a selling price of birr
6000 for each unit with total variable costs and selling estimates are birr 2000 for each
unit.
From the previous experience, it has been possible to drive a probability distribution
relating to the proportion of customers who will buy the product as follows:
a) show the expected opportunity loss value using the opportunity loss table
b) state whether or not the company should develop the product
Solution
=400000p- 60,000
The conditional profit values (payoffs) for each pair of state of nature (proportion of
customers who will by the product) and alternative (to develop or not to develop the new
product) is shown in the following table.
S1 S2
0.04 -44,000 0
0.08 -28000 0
0.12 -12000 0
0.16 4000 0
0.20 20,000 0
State of nature
(Proportion of customers, p)
c) Using the given estimate of probabilities of associated with each state of nature, the
expected opportunity loss (EOL) for each alternative is given below:
Since the company wants to minimize the expected opportunity loss, it should select
courses of action S2 (do not develop the product) with minimum EOL.
In decision making under risk each state of nature is associated with the probability of its
occurrence. However, if the decision maker can acquire perfect(complete and
accurate)information, about the occurrence of various state of nature , then he will be able
to select a course of action that yields the desired payoff for whatever state of nature that
actual occurs.
EVPI represents the maximum amount of money the decision maker has to pay to get
information about the occurrence of various state of nature before a decision has to be
made.
Example
A company needs to increase its production beyond its existing capacity. It has narrowed
the alternatives to two approaches to increase the production capacity: expansion, at a
cost of birr 8 million or modernization at a cost of birr5 million. Both approaches would
require the same amount of time for implementation. Management believe that over the
required payback period, demand will either be high or moderate. Since high demand is
considered to be somewhat less likely than moderate demand, the probability of high
demand has been set at 0.35. If the demand is high, expansion would gross an estimated
additional birr 12 million but modernization only an additional birr 6 million, due to
lower maximum production capability. On the other hand, if the demand is moderate, the
comparable figure would be birr 7 million for expansion and birr 5 million for
modernization.
d) Construct the conditional opportunity loss table and also calculate EOL.
Solution
Given: - state of nature: high demand (N1) and moderate demand (N2)
NB: since the probability that the demand is high estimated at 0.35, the probability of
moderate demand must be (1-0.35) = 0.6
b) The calculation of EMV for each courses of action S1 and S2 is given below:
To maximize EMV, the company must choose courses of action S1 (expand). Thus, the
EMV=0.75million
c) To calculate EVPI, we shall first calculate expected payoff with perfect information
(EPPI). To calculate EPPI, we chose optimal courses of action for each state of nature ,
multiply its conditional profit by the given probability to get weighted profit, and then
sum these weights as shown below:
EPPI =1.40
The optimal EMV is 0.75million corresponding to the courses of action S1. Then,
= 1.40-0.75
(Demand) probability S1 S2 S1 S2
Moderate DD 0.65 -1 0 1 0
Since decision makers seek to minimize the EOL, he must select alternativeS1.
Decision threes
In the discussion of decision problems until now our concern has been with the single
stage a problem where in the decision maker has to select the best course of section on
the basis of whatever information is achievable at appoint of time. We shall now consider
the decision situations that involve multiple stages. Also called the sequential decision
problems, they are characterized by a sequence of decisions in which following each
decision, a chance event occurs which in turn affects the next decision.
In analyzing multiple stage decision situations we have to evaluate the decision procedure
in a backward manner by evaluating the best course of action at the later stages to decide
the best action at the earlier stages. For this purpose the decision tree or the decision flow
diagram as it is sometimes called is a very effective device.
Example1.
1. An oil company has recently acquired rights in a certain area to conduct surveys and test
drillings to lead to lifting oil if it is found in commercially exploitive quantities. The area
is considered to have good potential for finding oil in commercial quantities. At the
outset the company has the choice to conduct further geological tests or to Cary out a
drilling program immediately. On the known conditions, the company estimates that
there is a 70:30 chance of future tests showing a success.
Whether a tests show the possibility of ultimate success or not or even if no test are
undertaken at all, the company could still pursue its drilling program or alternatively
consider selling its rights to drill in the area. Therefore, However, if it carries out the
drilling program, the likely hood of final success or failure is considered depends on the
forgoing stages. Thus;
If successful tests have been carried out the expectation of success in drilling is
given as 80:20.
If the test indicates failure, then the expectation of success in drilling is given as
20:80.
If no tests have been carried out at all the expectation of success in drilling is given
as 55:45.
Costs and revenues have been estimated for all possible outcomes and the net present
value of each is as follows:
Success
Failure
A. Draw up a decision (probability) tree diagram to represent the above information: and
B. Evaluate the tree in order to advice the management of the company on its best course of
action.
-40 65
-50
Failure (0.45)
Sell l Failure (0.2)
Index
Decision nods
Chance nods
As you observe that the tree has several branches which originate from squares or from
circles. A squire represents a decision node or decision fork at which the decision maker
has to take a decision, while a circle represents a chance node or chance fork at which
events (i.e. the state of nature) are branched out.
Total (20)
AMU, CBE, Department of Management December 2010 Page 205
2. Sell 1.0 15 15
Operations Research
The expected value of the drilling option is a loss of Birr 20 million while if the site is
sold, we can get birr 15 million. On the basis of the criterion of maximization of the
expected profit, our decision would be to sell the site.
Total 70
2. Sell 1.0 65 15
The two alternative course of action, selling and drilling have expected value equal to 65
million and birr 70 million, respectively. Obviously, therefore, provided that a positive
result is indicated by the test, the best course would be to go in for oil drilling.
Now at this stage we have two conditional decisions, sell the site if negative is obtained
on the test and drill in case the test indicates a positive. At each node the branch at which
we have not to move. Next, we move to the decision node where a decision has to be
taken whether to drill at the outset.
Total 48
Positive 0.7 70 49
Total 53.5
3. Sell 1.0 45 45
The expected value of the alternative of carrying out a test is birr 53.5 million, which is
the highest of the three. Therefore, it is better to test before carried out. If it proves
negative, the right should be sold to give a return of birr 15 million. To proceed with
relining, if that happens, would expectedly lead to a loss. However, if the test proves
positive the drilling should be under taken.
Chapter Summary
Knowledge of the likely hood of each of the state of nature can play an important role in
selecting a course of action. There are three decision making environments: certainty,
uncertainty and risk.
Decision under certainty is the case where the decision maker has a perfect knowledge
(information) about the state of nature
Decision making under uncertainty is a situation in which decision maker has neither
previous knowledge of the outcomes of his decisions nor is in a position to attribute a
probability of occurrence of such outcomes.
The most widely used criterion for evaluating various courses of action (alternatives)
under risk includes: expected monetary value (EMV), expected opportunity loss (EOL),
expected value of perfect information (EVPI),
Expected monetary value is the sum of the payoffs for each course of action multiplied
by the probabilities associated with each state of nature.
Expected value of perfect information (EVPI) represents the maximum amount of money
the decision maker has to pay to get information about the occurrence of various state of
nature before a decision has to be made.
Review question
a. Certainty c. Risk
b. Uncertainty d. All
3. Which of the following criteria is not applicable to decision making under risk
4. Which of the following criteria is not used for decision making under uncertainty
a. Maximin c. Minimax
b. Maximax d. Minimize expected loss
Regular sharing 2 5 8
Risky sharing -5 7 15
Property -10 10 20
Over the past 300 days, 150 days have been medium market conditions and 60 days have
had high market increases.
On the basis of these data, state the optimal investment strategy for the investor.
According to the f given information the probabilities of low medium and high market
conditions would be 90/300 0r 0.30, 150/300 or 0.50 and 60/300 or 0.20, respectively.
Answer b, regular share 0.047, b. risky shares 0.050 and property 0.060. c. since the
expected of 6% highest for property, the investor should invest in this alternative.
CHAPTER 5
NETWORK MODEL
Learning objectives
5.1. Introduction
Network model is a technique used for planning, and scheduling large projects in the
field of construction, maintenance, fabrication, purchasing computer system, distribution,
transportation, facility allocation, and many others. There are several reasons for the wide
range of applications. First many problems that can be formulated as LP models can also
be visualized with the context of networks. Second network models of problems are
often easier to formulate than alternative LP models of the same problem. Third the
pictorial representation of network models makes them relatively easy to understand,
particularly by the nontechnical decision makers. Finally the advances in network
optimization procedure have resulted in some amazing results in solving network model.
Network modeling: Techniques that enable complex projects to be scheduled, taking into
account the precedence of each activity.
The different systems satisfying the definition of network in the physical world
PROJECTS
Examples of a project:
Characteristics of a project:
i. Identify various activities (tasks or work elements) to be performed in the project, that is,
develop a breakdown structure (WBS).
ii. Determine requirement of resources such as man, material, machines, money, e.t.c. for
carrying out activities listed above.
iii. Assign responsibility for each work package. The work package corresponds to the
smallest work effort defined in a project and from the set of elements which are the base
for planning, scheduling and controlling the project
iv. Allocate resources to work package.
v. Estimate cost and time at various level of project completion.
vi. Develop work performance criteria
vii. Establish control channels for project personnel.
v. Based on the time estimate, calculate the total project duration, identify critical path;
calculate floats; carry out resource smoothing (leveling) exercise for critical (or scarce)
resources taking in to account resource constraints if any.
There are two basic planning and control technique that utilize network to complete
predetermined Project or schedule. These are Program Evaluation and Review Technique
(PERT) and Critical Path Method.
PERT is an acronym for "Program Evaluation and Review Technique". This was created
as a means to plan and accelerate the development of the Polaris Ballistic Missile. In
USA the defense department developed a nuclear missile to be launched from beneath the
ocean's surface by a mobile submarine, which would be an effective deterrent against
aggression by an enemy. This paved way to plan how to design, develop and plan the
different stages in the production of a missile and how quickly this task could be
completed. A planning and scheduling technique named PERT gave the answer to these
questions.
In any new venture, uncertainties are bound to creep in. PERT incorporated these
uncertainties into a model, which provides a reasonable answer to these uncertainties.
There are certain statistical aspects scheduling large projects consisting of numerous
activities whose completion times are uncertain and are independent of one another.
PERT is an event-oriented technique. By 'event' we mean reaching a certain stage of
completion of the project.
If we employ more persons or skilled people or given overtime to the workers, the
activity could be completed in a reduced duration known as crash duration. But this
involves an increased cost in the form of additional resources. With CPM the amount
needed to complete the various activities is assumed to be known with certainty. So, the
direct costs for the activities increase and hence the cost of the project also increases. By
reducing the activity duration of some or all possible completed ahead of the schedule.
This will naturally reduce the overhead cost for the entire project. On one hand the direct
expenses increase, if we shorten the activity duration, but, the indirect expenses for the
project are reduced. We have to strike a balance or an optimum time schedule, or at least
cost schedule is to be obtained. This is the purpose of the Critical Path Method. Thus
CPM is not concerned with uncertain job times as in PERT. PERT is useful in research
and developmental projects, whereas CPM is mostly used in construction projects, or in
situations already handled, so that the details like the normal completion time, crash
duration and cost of crashing are already known.
The following are the suggested applications when PERT or CPM is found useful.
1. It assumes a probability distribution for the duration of each activity. Thus completion
time estimates for all of the activities are needed
2. To perform PERT analysis on a project, the emphasis is given on the completion of a task
rather than the activities required to perform to reach a particular event or task. Thus, it is
also called event-oriented technique.
3. It is used for one time project involving activities of none-repetitive nature (i.e. activities
which may never have been performed before) in which time estimates are uncertain,
such as redesigning an assembly line or installing a new information system
4. It helps in identifying critical areas in a project so that necessary adjustment can be made
to meet the schedule completion date of the project.
CPM
1. This technique was developed in connection with a construction and maintenance project
in which duration of each activity was known with certainty.
2. It is suitable for establishing a trade-off for optimum balancing between schedule time
and cost of the project.
3. It is used for completion of project involving activities of repetitive nature.
4. Activity times may be estimated either as a single point estimate (CPM) or as a 3- point
estimate (PERT)
CPM- Critical Path Method
PERT- Program Evaluation and Review Technique
5. The duration of the activities is independent of each other.
Basic terms
Events: Events of the network represent protect milestone, such as the start or the
completion of an activity (task) or activities, and occur at a particular instant of time at
which some specific part of the project has been or is to be achieved. Event An event is a
point in time that marks the beginning or ending of activity. Events are commonly
represented by circles (nodes) in the network diagram.
The events can be further classified into the following two categories (types):
i. Merge event
An event which represents the joint completion of more than one activity is known as
merge event.
Event
Events in the network diagram are identified by numbers. Each event should be identified
by a number higher than that allotted to its immediately preceding event to indicate
progress of work. The numbering of events in the network diagram must start from left
(start of the project) to right (completion of the project) and top to the bottom. Care
should be taken that there is no duplication in the numbering.
Activities: Activity is a time-consuming job or task that is a key subpart of the total
project. Activities of the net work represent project operations or tasks to be conducted.
As such each activity except dummy consumes time and resources and incurs costs. An
arrow is commonly used to represent an activity with its head indicating the direction of
progress in the project.
Activities are identified by the numbers of their starting (tail or initial) event and ending
(head or terminals) event. An arrow (i , j) extended between two events; the tail event i
represents the start of the activity and the head event j, represents the completion of the
activity.
Activity
A dummy activity in the network is added only to represent the given precedence
relationships among activities of the project and is needed when
a) Two or more parallel activities in a project have same head and tail events, or
b) Two or more activities have some (but not all) of their immediate predecessor
activities in common.
A B C
Activity Activity Activity
1 2 4
D Dummy Activity
B is the predecessor of C and D Activity
3
Network model use the following two type of network to show precedence requirements
of the activities in the project.
Activity-on-Node (AON) Network in this type of precedence network each node (or
circle) represent a specific task while the arcs represent the ordering between tasks. AON
network diagrams place the activities when the nodes and the arrows are used to indicate
sequencing requirement. Generally, these diagrams have no particular starting and ending
nods for the whole project. The lack of dummy activities in this diagram always makes
them easier to draw and to interpret.
Activity-on-Arrow (AOR) network: in this type of precedence network at the each end
of the activity arrow is a node (or circle). This node represents points in time or instance,
when an activity is starting or ending. The arrow itself represents the passage of time
required for that activity to be performed.
These diagrams have a single beginning node from which all activities with no
predecessor may start. The diagram then works its way from left to right, ends with a
single ending node, where all activities with no followers come together. Three important
advantage of using AOA are:
1. Activity A A
A
A B
2. B must follow A A B
B
A
C
3 C must follow B
A&B C
A
4 C must follow A, A C
and D must follow
A&B
B D
1. In network diagram, arrows represent activities and circles the events. The length of an
arrow is of no significance.
2. Each activity should be represented by one arrow and must start and end in a circle
called event. The tail of an activity represents the start and head the completion of work.
3. The event numbered 1 denotes start of the project and is called initial event. All activities
emerging (or taking off) from event 1 should not be preceded by any other activity or
activities. Event carrying the highest number denotes the completion events. A network
should have only one initial event and only one terminal event.
4. The general rule of numbering the event is that the number at an activity’s head should
always be larger than that at its tail. That is, events should be numbered such that for
each activity (i, j), i< j.
5. An activity must be uniquely identified by its starting and completion event which
implies that
a. An event number should not get repeated or duplicated
b. Two activities should not be identified by the same completion event
c. Activities must be represented either by their symbols or by the
corresponding ordered pair of starting-completion events.
6. The logical sequence ( or interrelationship ) between activities must follow the following
rules :
a. An event cannot occur until all the incoming activities into it have been
completed.
b. An activity cannot start unless all the preceding activities, on which it depends,
have been completed.
c. Though a dummy activity does not consume ether any resource or time, even then
it has to follow the rules 6(a) and 6(b).
1. Looping
A case of endless loop in a network which is also known as looping is shown in figure
below, where activities A, B, and C from a cycle.
4
3
C
1
Looping
B
A
2
In this case it is difficult to number three events associated with activity A, B, and C so as
to satisfy rule 6 of constructing the network.
2. Dangling
A case of disconnect activity before the completion of all activities which is also known
as dangling.
1 A 2 B 4
≈
C Dangling
In this case, activity C does not give any result as per the rules of the network. The
dangling may be avoided by adopting rule 5 of constructing the network.
When two or more parallel activities in a project have the same head and tail events:
i.e: two events are connected with more than one arrow.
In figure a, activities B and C have a common predecessor activity A. At the same time,
they have activity D as a common successor. To get the correct network a dummy
activity for the ending event of B to show that D may not start before both B and C are
completed is shown in figure b). A dummy which is used in such a case to establish
relationship proper logical relationships is also known as logic dummy activity.
B 2
A 1 B 3 D
A 2 D B Dummy
1
C Dummy A 1 C 3 D
2
C
C
Example
1. An assembly is to be made from two parts X and Y. Both parts must be turned on
a lathe and Y need not be polished. The sequence of activities together with their
predecessors is given below:
D Turn X on Lathe B
E Turn Y on Lathe B, C
F Polish Y E
G Assemble X and Y D, F
H Pack G
Solution
3
B D
1 2
A
Dummy 6 G 7 H 8
C
F
4 5
E
2. Listed in the table are the activities and sequencing necessary for a maintenance job on
the heat exchangers in a refinery.
_______________________________________________________________________
D Clean bolts B
G Clean shell C
Required:
Solution
5
F D2
C 4
A 2 B 3 G 6
1 H 8 J 10
D I 9
E
7 D1
A -
B -
C -
D A
E C
F A
G D, B, E
Solution:
4 5
1
Example A project has the following activities. The relationships among the activities
are given below. Construct the network.
D, E and F follow B.
G follows E.
I and J succeed G.
F and J precede K.
H and I precede L.
M succeeds L and K.
H
D
L
I
E
G J K M
B F N
A
C
The duration of individual activities may be uniquely determined (in case of CPM) or
may involve the three time estimates (in case of PERT) out of which the expected
duration of an activity is computing. Having done this, the following factors should be
known to prepare project scheduling.
In a network diagram, there should only be one initial event & one final event, while
other events are numbered consecutively with integers 1, 2… n such that i< j , for any
two event i&j connected by an activity which starts at i & finish at j
For calculating the above mentioned times, there are two methods; namely;
1) Set the earliest occurrence time of initial event 1 to zero, i.e. E1 ¿ 0 , i=1
2) Calculate earliest start time for each activity that begins at event i(=1) This is equal to
earliest occurrence time of event, i (tail event). That is ESij =E j , for all activities (i, j )
starting at event i.
3) Calculate the earliest finish time of each activity that begins at event i. This equal to the
earliest start time of the activity plus the duration of the activity. That is,
EFij =ES ij +t ij =Ei +t ij for all activities ( i j ) beginning at event
4) Proceed to the next event, say j; j > j
5) Calculate the earliest occurrence time for the event j. This is the maximum of the
earliest finish times of all activities ending into that event, i.e.
Ej=max { EF ij }=max { E i+ t ij } , for all immediate predecessor activities.
6) If j = N (final event ), then earliest finish time for the project, i.e. the earliest
occurrence time E N for the final event is given by
E N =max { EF ij }=¿, for all terminal activities.
1. Set the latest occurrence of last event N equal to its earliest occurrence tome ( known from
forward pass method), that is
LN =E N , j=N
2. Calculate latest finish time of each activity which ends at event j. This is equal to latest
occurrence time of final event N. That is,
LF ij = L j , for all activities (i, j) ending at even j.
3. Calculate the latest start times of all activities ending at j. I t is obtained by subtracting the
duration of the activity from the latest finish time of the activity, That is,
LSij =LF ij – t ij = L j –t ij, for each activity (i, j) ending at Event j,
6. If j = 1 (initial event), then the latest finish time for project, i.e. latest occurrence time L1
for the initial event is given by.
Li=min {LS ij } For all immediate successor activities
1. Event Float
The float (also sometimes called ‘slack’) of an event is the difference between its latest
time (Li) and its earliest time ( E i).That is
Event Float=Li−E i
It is a measure of how much later than expected a particular event could occur without
delaying the completion of the entire project.
2. Activity Float
a. Total Float
It is the difference between the time available to perform the activity (measured from the
earliest start time to the latest finish time) and the expected completion time of the
activity. That is, for activity (i, j), the total float is given by
¿ LSij – ESij
¿ LF ij – EF ij
b. Free Float
For calculating the total float, only a particular activity was considered with respect to its
tail and head event times or by considering latest start and finish time of an activity with
respect to its earliest start and finish time. However, it may be needed to know that how
much an activity’s scheduling time can be increased or decreased without changing its
immediate successor activities.
The time by which the completion of an activity can be delayed beyond its earliest finish
time without affecting the earliest start time of a succeeding activity is called free float.
Remarks:
1. Since latest event time is always greater than or equal to the earliest event time (i.e
Li> ¿Ei¿ ,(TFij )> (FF〗 ij).
This implies the value of free float may range from zero to total float and can never
exceed total float value.
2. The calculation of various floats can help the decision-maker in identifying the
underutilized resources, flexibility in the total schedule and possibilities of
redeployment of resources.
3. Once the float of an activity is disturbed, float of all other activities of the project is
changed and should be recalculated.
-resources are not adequate and activity may not finish in time
-Induce extra resources or certain activities (also called critical activities) need crashing
to reduce negative float value.
-Resources are just sufficient to complete the activity. Any delay in activity execution
will necessarily increase the project cost.
-Surplus resources can be deployed elsewhere or execution of the activity can be delayed.
Critical Path
Path: A path is unbroken chain of activities that connects form the start node to the end
node (terminal node).
Certain activities in a network diagram of a project are called critical activities because
delay in their execution will cause further delay in the project completion time. Thus, all
activities having zero total float value are identified as critical activities.
The critical path is the continuous chain of critical activities in a network diagram. It is
the longest path starting from first to the last event & is shown by the thick line or double
lines in the network diagram most of the time.
The length of the critical path is the sum of the individual times of the critical activities
living on it & defines the minimum time required to complete the project.
i. For all activities (i, j) lying on the critical path the E –values & L-values for tail &
head events are equal That is, E j=L j ∧Ei=Li
ii. On critical path E j−Ei =L j−Li=Lij
Finding the critical path is important for directing decision-makers attention & effort to
critical activities where improvement will pay the largest dividend.
Example
1. An established company has decided to add a new product to its line. It will buy the
product form manufacturing concern, package is & sell it to a number of distributes
selected on a geographical basis. Market research has identified the volume expected &
the size of sales force required. The steps
B Hire Salesman 4
C Train salesman 7
F Conduct 10
G Design package 2
K Select distributors 9
Required:
L Sell to distributors 3
a) Draw the
M Ship stocks 5
an arrow
diagram
b) Indicate the critical path; the critical activities & the total project time
c) For each non-critical activity, find the total float
Solution: (a). the arrow diagram for the given project along with E-values and L-values
is shown in the following diagram.
E4=13
L4=14 4 I (6)
)
J (13 H (10)
M (5)
G (2) 9 10
1 3
6 E9=20
E1=0 A (6) E3=2 E10=25
L9=20
L1=0 L3=4 K (9) L10=25
AMU, CBE, Department of Management December 2010 E6=17 Page 234
2 C (7) L6=17 F (10)
5
B (4) E5=10
E2=6 L5=10
Operations Research
Determine the earliest start time E i and latest finish Lj for each event by processing as
follows.
E9=Max {Ei+ti, 9} =Max {E4+t4, 9;E6+t6, 9} =Max {E8+t8, 10; E9+t9, 10}
J=5, 6, 7 j=2, 3, 4
=min {L5, - t2, 5; L6, t2, 6;L7-t2, 7} =min{L2,- t1, 2; L3, t1, 3;L4-t1, 4}
B. The vertical path in the network diagram above has been shown by thick lines by joining
all those events where two values Ei& Lj are equal. The critical path of the project is: 1-
2-5-6-9-10 the critical activities are A, B, C, L&M. the total project time is 25 weeks
C. For each non-critical activity, the total float is:
1-3 2 0 2 2 4 2 0
1-4 13 0 13 1 14 1 0
2-6 9 6 15 8 17 2 2
2-7 2 6 8 9 11 3 0
3-4 10 2 12 4 14 2 1
4-9 6 13 19 14 20 1 1
7-8 4 8 13 11 15 3 0
8-10 10 12 22 15 25 3 3
Remark:
Example2. The following table gives the activities of a construction project and
duration.
Duration (days) 20 25 10 12 6 10
(iii) Find the total, free and independent floats of each activity.
Solution: The first step is to draw the network and fix early start and early finish
schedule and then late start-late finish schedule.
E1=24
E1=0 20 2 12
L1=36
L1=20 E1=36
E1=20 L1=46 5
4
1 15 L1=30 10
6
25
E1=5 3 E1=30
L1=30 L1=36
AMU, CBE, Department of Management December 2010 Page 237
Operations Research
1-2 0 0 0
1-3 5 5 5
2-3 0 0 0
2-4 4 4 4
3-4 0 0 0
4-5 0 0 0
To find the critical path, connect activities with 0 total slack and we get 1-2-3-4-5 as the
critical path.
1-2-4-5 = 42
1. Optimistic time (to or a): This is the shortest possible time to perform an activity,
assuming that everything goes well.
2. Pessimistic time (tp or b): This is the maximum (longest) time that is required to perform
an activity, under extremely bad conditions. However, such conditions do not include
acts of God like earthquakes, flood, etc.
3. Most likely time (tm or m): This is the most realistic time to complete the activity.
Statistically, it is the modal values of duration of the activity i.e. consider the time
estimate under the normal situation.
The β distribution is not necessarily systematic; the degree of skewness depends on the
location of t m ¿ t o andt p. Thus the range specified by the optimistic time (t0) and
pessimistic time (tp) estimates is assumed to enclose every possible estimates of the
duration of the activity. The most likely time (tm) estimate may not coincide with the
midpoint ( 〖 to+ tp)/2 and may occur to its left or its right as shown in the following
figure.
Pr
ob
ab
ilit
y
Activity duration
In beta distribution the midpoint (t o +t p )/2 is assumed to weight half as much as the most
likely point (t m). Thus the expected or mean (te) or μ) value of the captivity duration can
be approximated as the arithmetic mean of (to+ tp)/2 and 2tm. That is
6 σ ≡ tp−¿∨σ=tp −¿/6
As we are expecting variability in the activity duration, the total project may not be
completed exactly in time. Thus, it is necessary to calculate the probability of actually
meeting the schedule time to the project as well as activities.
The probability distribution of time for completing an event can be approximated by the
normal distribution due to central limit theorem. Thus, the probability of completing the
project by schedule time (T s) is given by
Prob(Z= te/σe )
Z= number so standard deviations the scheduled time or target date lies away
from the mean or expected date.
In order to find out the probability of completing the project in some given time, we shall
consider only the expected length of the critical path and its variance. The expected time
of the project can be calculated by adding the expected time of each activity lying on the
critical path can be known by adding variance of critical activity.
Example1.
1. A project is represented by a network shown below & has the following data
Task: A B C D E F G H I
Optimistic time 5 18 26 16 15 6 7 7 3
Pessimistic time: 10 22 40 20 25 12 12 9 5
F
I
3 6 7
B
E H
A 2 D 5
1
G
C
4
d. The probability of an event occurring at the expected completion date if the original
scheduled time of completing the project is 41.5 weeks.
e. The duration of the project that will have 95% chance of being completed.
Solution (a) using the following formula, the expected activity time (te) or μ) and
variance (σ^2) is given in the following table.
(b) The earliest and latest estimated time presented in the following graph.
E6=29
E3=20
L6=38.8
L3=29.8
F (9) E7=42.8
I (4)
3 6 7
E5=25.8 L7=42.8
E (20) H (8)
B (20)
E2=7.8 L5=34.8
A (7.8) 2 D (18) 5
1
L2=16.8
G (9.8)
E1=0
L1=0 C (33) 4
E4=33
L4=33
(c). the critical path is shown by thick line in the above figure where E values and L
values are the same. The critical path is 1-4-7 and the critical completion time for the
project is 42.8 week.
(d) the last event 7 will occur only after 42.8 weeks. For this, we require only the duration
of critical activity. This will help us in calculating the standard of the duration of the last
event.
distribution table
Thus, the probability that the project can be completed in less than or equal to 41.5 weeks
is 0.30. In other words, the probability that the project well gate delayed beyond 41.5
weeks is 70%
Project crashing
Crashing is employed to reduce the project quickly when time is reduced. Similarly,
completion time by spending extra beyond point B, the time increases while the
resource (cost). However, as shown in the cost decrease. Since for technical reasons
figure beyond point A, cost increases more time may not be reduced indefinitely,
therefore, we call this limit as crash point. provides the starting point for crashing
There is also a cost efficient duration called analysis.
normal point. Thus extending the activity b) When all critical activities are crashed.
duration beyond normal point may increase This provides the stopping point for
costs. crashing analysis.
Step2. Identify critical activity and compute
For simplicity the relationship between
the cost slope for each of these by using the
normal time and cost as well as crash time
relationship
and cost for an activity is assumed to be
linear instead of being concave and/or Cost slope=(crash cost-normal cost)/(crash
discrete. Thus the crash cost per unit of time time-normal time)
can be estimated by computing the relative
The values of cost slop for critical activities
change in cost (cost slope) per unit change
indicate the direct extra cost required to
in time.
execute an activity per unit of time.
Remark: crashing an activity means
performing it in the shortest technically
possible time by allocating to it necessary
resources.
Activity time
Operations Research
Example1. The following table gives data on normal time, and cost and crash time and cost for a
project
2-3 3 30 3 30
4-5 0 0 0 0
a) Draw the network and identify the critical path with a double line.
b) What are the normal project duration and associated cost?
c) Find out the total float associated with each activity.
d) Crash the relevant activities systematically and determine the optimal project completion
time and cost.
Solution (a) the network for normal activity times indicates a project completion time of 32
weeks with the critical path: 1- 2- 5- 6 – 7 -8, as shown in the following figure.
E4=10
E7=22
L4=12
L7=22
4 E5=12 E6=18 7
0
7 L5=12 L6=18
AMU, CBE, Department of Management December 2010 Page 246
Operations Research
1 3 9 6 10
2 5 6
E1=0 E2=3
3 13
5
L1=0 L2=3
8
3
E3=6 E8=32
L3=7 L8=32
(b). normal project duration is 32 weeks and the associated cost is as follow:
(c). Calculation for total float associated with each activity are shown in the following table.
(Li−Ei)−tij
(d). for critical
1-2 (3−0)−3=0
activities, crash cost
slope is 2-3 (7−3)−3=1 given.
2-4 (12−3)−7=2
E4=10
E7=21
2-5 (12−3)−9=0
L4=11
4 L7=21
E6=17 7
3-5 0(12−6)−5=1
7 E5=11
L6=17 4
3 4-5 8 (12−10)−0=2
L5=11 10
1 2
× 5
6 6
7-8 (32−22)−10=0
The minimum value of crash cost per week is for activity 2-5 and 5-6. Hence, crash activity 2-5
by 2days from 9 weeks to 7 weeks. But the time should be reduced by 1 week only otherwise
another path 1-2-3-4-6-7-8 becomes a parallel path. Net work as shown in the above figure is
developed when it is observed that the new project time is 31 weeks and the critical paths are 1-
2-5-6-7-8 and 1-2-3-5-6-7-8. The following table shows the crash cost per week (in birr)
1-2 400−300
=100
3−2
2-5 810−720
=45
9−7
5-6 410−320
=45
6−4
6-7 470−400
=70
4−3
7-8 1200−1000
=200
10−9
With crashing of activity 2-5 the new total cost involved can be calculated as follows;
New Total cost =Total Direct no rmal cost+increased direct cost due ¿ crashing of activity
Dear learners now with respect to the new network given, the new possibilities for crashing in
the critical paths are listed in the following table.
1-2 100
2-5 × (Crashed)
6-7 50
7-8 45
1-2 70
2-5 200
The minimum value of crashed cost slope is for activity 5-6. Hence, crashing it by 2 weeks from
6 weeks to 4 weeks. The new network diagram is presented as follow. E7=19
E4=10
4 E6=15 L7=19
7
L4=11 0
7 E5=11 L6=15 4
3 8 L5=11 4 10
1
E2=3
2
× 5 × 6
E1=0 3 13
5
L2=3 E8=29
E3=6 8
L1=0
3
L3=6 L8=29
As you see in the above figure the critical paths are not changed because activity 5-6 is common
for the two critical paths. But with crashing of activity 5-6 by 2 weeks, the new cost involved is:
New Total cost =Total Direct normal cost +increased direct cost due ¿
With respect to the new network depicted above the new possibility for crashing in the critical
paths are listed below.
1-2 100
2-5 × (Crashed)
5-6 × (Crashed)
6-7 70
7-8 200
The further crashing of 6-7 activity time from 4 weeks to 3 weeks will results in increased direct
cost than the gain due to reduction in project time. Consequently, here we must stop further
crashing. The optimal project duration is 29 weeks with associated cost of birr 5,805 as shown in
the table below.
Summery
Network model is a technique used for planning, and scheduling large projects in the
field of construction, maintenance, fabrication, purchasing computer system, distribution,
transportation, facility allocation, and many others.
A project is a set of proposals for the investment of resources in to a clearly identified set
of actions that are expected to produce future benefits that exceeds costs.
There are two basic planning and control technique that utilize network to complete
predetermined Project or schedule. These are Program Evaluation and Review Technique
(PERT) and Critical Path Method.
Critical path analysis helps to estimate the total project duration & to assign starting &
finishing times to all activities involved in the project. This helps in checking actual
progress against the scheduled duration of the project.
List cost approach, also called the minimum –cost method, uses lowest cell as the basis
for selecting route.
The assignment problem is a special case of Transportation problem in which the number
of sources and destinations are the same, and the objective is to assign the given job
(tasks) to most appropriate machine (person)so as to optimize the objective function like
minimizing cost or maximizing profit.
Review Exercise
6. Listed in the table are the activities & sequencing necessary for the completion of a
recruitment procedure for management trainees (MT) in a firm.
A -
B A
C A
D C
E B
F C
G D, E
H G
I H
J F
K I, J
L K
M L
N M, N.
7. The owner of a chain of fast food restaurants is considering a new computer system for
accounting and inventory control. A computer company sent the following information
F Develop input-output C 8 9 16
routines
G Create database E 4 8 12
Required:
A. Construct an arrow diagram for this problem
B. Determine the critical path and compute the expected completion time
C. Determine the probability of completing the project in 55 days
Chapter Six
Queuing Theory
Objectives:
6.1. Introduction
The queuing theory also called the waiting line theory, owes its development to A K
Erlang’s effort to analyze telephone traffic congestion with a view to satisfying the
randomly arising demand for the services of the Copenhagen automatic telephone system,
in year 1909. The theory is applicable to situations where the “customers” arrive at some
service stations for some service, wait (occasionally not): and then leave the system after
getting the service.
Queues or waiting lines arise when the demand for a service facility exceeds the capacity
of that facility, that is, the customers do not get service immediately upon request but
must wait, or the service facilities stand idle and wait for customers.
Some customers wait when the total number of customers requiring service exceeds the
number of service facilities, some service facilities stand idle when the total number of
service facilities exceeds the number of customers requiring service.
Queuing theory can be applied to a wide variety of operational situations where imperfect
matching between the customers and the service facilities is caused by one’s inability to
predict accurately the arrival and service time of customers. In particular, it can be used
to determine the level of service (either the service rate or the number of service
facilities) that balances the following two conflicting costs:
An increase in the existing service facilities would reduce the customer’s waiting time.
Conversely, decreasing the level of service should result in long queue(s). This means an
increase (decrease) in the level of service increases (decreases) the cost of operating
service facilities but decreases (increases) the cost of waiting time.
Figure1.1. Illustrates both types of costs as a function of level of service. The optimum
service level is one that minimizes the sum of the two costs.
Cost of operation
Total expected cost
Level of service
Of the two types of costs mentioned above, the cost of waiting is difficult to estimate.
However, it is measured in terms of loss of sales or goodwill when the customer is a
human being who has no sympathy with the service. But, if the customer is a machine
waiting for repair, then cost of waiting is measured in terms of the cost of lost production.
The waiting line is developed because the service to a customer may not be rendered
immediately as the customer reaches to the service facility. Thus lack of adequate service
facility would cause waiting lines of customers to be formed. The only way that the
service demand can be met with ease is to increase the service capacity (raising the
efficiency of the existing capacity if possible) to a higher level. A manager has to decide
an appropriate level of service which is neither to low nor too high.
1. Input source
2. Queuing process
3. Queue discipline
4. Service process(or mechanism)
An input source need not be homogeneous population but may consist of several sub-
populations.
Behavior of arrivals: if a customer, on arriving at the service system stays in the system
until served, no matter how much he has to wait for the service is called a patient
customer. Whereas the customer, who waits for a certain time in the queue and leaves the
service system without getting service due to certain reasons such as a long queue in front
of him is called an impatient customer.
Pattern of arrivals: the arrival process (or pattern) to the service system is classified into
two categories: static and dynamic. These two are further classified based on the nature
of arrival rate and the control which can be exercised on the arrival process.
In static arrival process, the control depends on the nature of arrival rate (random or
constant). Random arrivals are either at a constant rate or varying with time. Thus to
analyze the queuing system, it is necessary to attempt to describe the probability
distribution of arrivals. From such distributions we obtain average time between
successive arrivals, also called inter-arrival time (time between two consecutive arrivals),
and the average arrival rate (i.e. number of customers arriving per unit of time at the
service system).
The dynamic arrival process is controlled by both service facility and customers. The
service facility adjusts its capacity to match changes in the demand intensity, by either
varying the staffing levels at different timings of service, varying service charges (such as
telephone call charges at different hours of the day or week) at different timings, or
allowing entry with appointments.
The variation in the demand intensity also affects the customer’s behavior. They either
balk or renege from the service system when confronted with a long slow moving waiting
line.
The arrival time distribution can be approximated by one of the following probability
distributions;
a) Poisson distribution
b) Exponential distribution
c) Erlang distribution
Frequently in queuing problems, the number of arrivals per unit of time can be estimated
by a probability distribution known as the Poisson distribution, as it adequately supports
many real world situations.
λt
P (×=n/Pn= λt) =( λt ) e −¿ ¿ , n=0, 1 , 2, … .
n
n!
Now, if we are interested in the inter arrival time probability distribution, then the time
interval from 0 to t, the probability of no arrival is given by
λt
e−¿ n e λt −¿ ¿
0
λt
Let us define the random variable T as the time between successive arrivals. Since a
customer can arrive at any time, T must be a continuous random variable. The probability
of no arrival in the time interval from 0 to t will be equal to the probability that T exceeds
t, so we have
λt
The probability that there is an arrival in the time interval from 0 to t is given by
λt
This probability is also called the cumulative probability distribution function, F (t) of T.
Also the distribution of the random variable T is referred to as the exponential
distribution. Its probability density function (pdf) is
f (t)= { λ e 0 ; t <0 }
λt
−¿ ;t ≥0 ¿
F(t)
f (t) λ=3.0
t
0 1 2 3 4
The mean of exponential distribution is the expected or average time between arrivals E
1
(T). Thus with λ arrivals per unit time, E(T )=
λ
Hence from the above discussion we may conclude that the Poisson distribution of
arrivals with arrival rate λ is equal to the negative exponential distribution of inter-arrival
1
times with mean inter arrival times .
λ
The probability density function of exponential distribution can also be used to compute
the probability that the next customer arrives within time, T of the previous arrival. This
means that if a customer has already arrived at the service system, then probability of
arriving for the next customer can be determined by using this equation.
= 1-0.165-0.835
Also, if customers arrive at the service system at an average rate of 24 customers per
hour, and a customer has already arrived, then probability of another customer arriving in
the next 5 minutes (i.e. t=1/12 hr) is
Further if λ =24 customers per hour, t= 1/12 hour and λ t=24 ×(1 /12)=2, the probability
of n customers arriving within the next 10 minutes is:
( 121 )
( ) =0.27
2
−24
24
e
12
P(x=2 / p n=2)=
2!
Queuing process
The queuing process refers to the number of queues, and their respective lengths. The
number of queues depends upon the layout of the service system. Thus there may be
single or multiple queues. Certain service systems adopt “take a number” policy to avoid
a formal queue to form; the length (or size) of the queue depends upon the operational
situation such as physical space, legal restrictions, and attitude of the customers.
In certain cases, a service system is unable to accommodate more than the required
number of customers at a time. No further customers are allowed to enter until space
becomes available to accommodate new customers. Such type of situations are referred to
us finite (or limited) source queue. On the other hand, if the service system is able to
accommodate any number of customers at a time, then it is referred to as infinite (or
unlimited) source queue.
In many other situations, when arriving customers experience long queue(s) in front of a
service facility, they often do not enter the service system even though additional waiting
space is available. The queue length in such cases depends upon the attitude of the
customers. In some finite source queue systems, the maximum permissible queue is of
zero length.
A multiple queue system at a service system can also be finite or infinite. But this has
certain advantages such as:
Queue discipline
In the queue structure, the important thing to know is the queue discipline. The queue
discipline is the order or manner in which customers from the queue are selected for
service.
There are a number of ways in which customers in the queue are served. Some of these
are:
a) Static queue disciplines: are based on the individual customer's status in the queue. Few
of such disciplines are:
i. If the customers are served in the order of their arrival, then this is known as the first-
come, first-served (FCFS) service discipline. Prepaid taxi queue at airports where a taxi
is engaged on a first-come, first-served basis is an example of this discipline.
ii. Last-come-first-served (LCFS) sometimes, the customers are serviced in the reverse
order of their entry so that the ones who join the last are served first. For example,
assume that letters to be typed, or order forms to be processed accumulate in a pile, each
new addition being put on the top of them. The typist or the clerk might process these
letters or orders by taking each new task from the top of the pile. Thus, a just arriving
task would be the next to be serviced provided that no fresh task arrives before it is
picked up. Similarly, the people who join an elevator last are the first ones to leave it.
b) Dynamic queue disciplines are based on the individual customer attributes in the queue.
Few of such disciplines are:
i. Service in Random Order (SIRO) -- Under this rule customers are selected for service
at random, irrespective of their arrivals in the service system. In this every customer in
the queue is equally likely to be selected. The time of arrival of the customers is,
therefore, of no relevance in such a case.
ii. Priority Service-- Under this rule customers are grouped in priority classes on the basis
of some attributes such as service time or urgency or according to some identifiable
characteristic, and FCFS rule is used within each class to provide service. Treatment of
VIPs in preference to other patients in a hospital is an example of priority service.
iii. Pri-emptive priority- under this rule, the highest priority customer is allowed to enter
into the service immediately after entering into the system even if a customer with lower
priority is already in service. That is lower priority customer’s service is interpreted (pre-
empted) to start service for a special customer. This interrupted service resumed again
after the highest priority customer is served or according to the pre-emptive rules.
iv. Non-pre-emptive priority- in this case highest priority customer goes ahead in the
queue, but service is started immediately on completion of the current service.
For the queuing models that we shall consider, the assumption would be that the
customers are serviced on the first-come-first-served basis.
The service process is concerned with the manner in which customers are serviced and
leave the service system. It is characterized by
The series arrangement consists of a sequence of a number of service facilities such that a
customer must go through one facility after another in a particular sequence before the
whole service is completed. Each service facility may however, work independently of
the others, having its own rules of service.
Service Facility
Customers
Served
Customers
Customers
Served
Customers
b) Single queue
In the parallel arrangement the incoming customer may join the queue of his choice in
front of service facilities or may be served by any service facilities.
Service facility
Served customers
Customers
Served customers
Arrangement of service facilities in parallel
Service facilities
Customer
Served
Customers
Distribution of service time-the time taken by the server from the commencement of
service to the completion of service for a customer is known as the service time.
If μ is the average service rate, then the expected number of customers served during time
interval t will be μt. Thus, if we consider zero time at the start of service, the probability
that services is not completed by time t is given by
−μ t
p(x=0/ p n=μ t)e
If the random variable T represents the service time, then the probability of service
completion within time t is given by
−μt
p ( T ≤ t ) =1−e ,t ≥0
It is assumed that customers joining a queuing system arrive in random manner and
follow a Poisson distribution or equivalently the inter-arrival times follow exponential
distribution.
In most of the cases, service times are also assumed to be exponentially distributed. It
implies that the probability of service completion in any short-time period is constant and
independent of the length of time that the service has been in progress. The basic reason
for assuming exponential service is that it helps in formulating simple mathematical
models which ultimately help in analyzing a number of aspects of queuing problems.
The number of arrivals and departures (those served) during an interval of time in a
queuing system is controlled by the following assumptions (also called axioms).
a. The probability of an event (arrival or departure) occurring during the time interval (t,
t+Δt) depends on the length of time interval Δt. That is, probability of the event does not
depend either on number of events that occur up to time t or the specific value of t,
meaning that the events that occur in non- overlapping time are statistically independent.
b. The probability of more than one event occurring during the time interval (t, t+ Δt) is
negligible. It is denoted by 0(Δt).
c. At most one event (arrival or departure) can occur during a small time interval Δt. the
probability of an arrival during the time interval ( t, t + Δt) is given by
The arrival process assumes that the customers arrive at the queuing at the queuing
system and never leave it. Such a process is called pure birth process. The aim is to
derive an expression for the probability P n (t) of n arrivals during time interval (t, t+ Δt).
The terms commonly used in the development of various queuing models are the
following:
Δt = a time interval so small that the probability of more than one customer’s arrival is
negligible, i.e. during any given small interval of time Δt only one customer can arrive.
λΔt = probability that a customer will arrive in the system during time Δt.
1- λΔt = probability that no customer will arrive in the system during time Δt.
If the arrivals are completely random, then the probability distribution of a number of
arrivals in a fixed time interval follows Poisson distribution.
( λ t )n −λt
pn ( t ) = e ,n=0 , 1 ,2 , … .
n!
The Markovian property of inter arrival times states that the probability that a customer
currently in service is completed at some time t is independent of how long he has
already been in service. That is
The departure process assumes that no customer joins the system while service is
continued for those who are already in the system. Let, at time t =0 (starting time) there
be N ≥ 1 customers in the system. Since service is being provided at the rate of μ.
Therefore, customers leave the system at the rate μ after being serviced. Such a process is
called pure death process
Basic axioms
t + Δ t is negligible.
The following terms are used in the development of various queuing models.
μ Δ t = probability that a customer in-services at time t will complete service during time
Δ t.
1 - μ Δ t = probability that the customer in-service at time t will not complete service
during time Δ t.
This shows that service times follows negative exponential distribution with mean
1/μ
2.
and variance 1 / μ
The service time of short duration have the highest probability of occurrence. As the
service time increases, the probability of occurrence tails off exponentially towards zero
probability.
The area under the negative exponential distribution curve is determined as:
{
F ( T )= t μ e
−μt t
−μ t dt = [− μ e 0 ]
0
−μt 0 −μ t
¿−e +e −¿ 1−e
− μt
F (T )=f (t ≤ T )=1−e
Where F(T) is the area under the curve to the left of T. Thus
−μ t
1−F (T )=f ( t >T )=e is the area under the curve to the right of T.
The derivation of this model is based on certain assumptions about the queuing system:
The following events (possibilities) may occur during a small interval of time, Δt just
before time t. it is assumed that the system is in state n (number of customers) at time t.
1) The system is in state n (number of customers) and no arrival and no departure, leaving
the total to n customer.
2) The system is in state n+1 (number of customers) and no arrival and one departure,
reducing the total to n customers.
3) The system is in state n-1(number of customers) and one arrival and no departure,
bringing the total to n customers.
We will begin with determining probability, Pn of n customers in the system at time t and
value of its various operating characteristics is summarized as follows:
0 1 2 n-1-- n n =1
1
If pn (t ) is the probability of a customers at time t in the system, then the probability that
the system will contain n customers at time ( t+ ∆ t ) can be expressed as the sum of the
joint probabilities of the three mutually exclusive and collectively exhaustive cases as
mentioned above. That is
For n ≥ 1∧t ≥ 0.
pn ( t ) { 1− λ ∆ t ¿ ( 1−μ ∆ t ) + pn+1 ( t )¿ ¿
p'
n (t )=
d { p n(t ) } =λ p n−1(t )+ μ pn +1(t ) – ( λ+ μ ) pn (t ) ;n ≥1 (1)
dt
p n(t +∆ t )− p 0( t )
=μ p 1 ( t )−λ p 0 (∆ t)
∆t 0 (t )+
∆t
p'
o ( t )=
d { p 0(t ) } =μ p 1(t)−¿ λ p (t ) ;n=0 ¿
0 (2)
dt
lim pn (t )= pn
n→∞
lim d
And n →∞
{ pn (t )} =0 ; n=0 , 1 ,2 … …
dt
Thus these equations constitute the system of steady-state difference equations. The
solution of these equations can be obtained either by (i) iterative method, (ii) generating
functions method, or (iii) linear operator method. Here we shall be using the iterative
method to find the values of P1, P2 ….in terms of P0, λ and μ .
λ p0 =μ p1 ∨p 1= ( μλ ) p ; n=0 0
0=−( λ+ μ ) p1+ ¿ λ p + μ p ¿ 0 2
p p p
Or 2=(
λ+μ
μ
) 1−
λ
μ
0=¿ ( λ+μ μ ) μλ p − λμ p
0 0=¿ ¿ ¿
p λ+μ p λ p1
3=( ) 2−( )
μ μ
=¿ ( λ +μ μ )( λμ ) 2 p ( μλ ) p λ
0−( )
μ
λ
0=¿( ) p0; n=2 ¿
μ
2
¿¿
To obtain the value of p0, we make use of the fact that (sum of all possibilities)
∞
∑ p n= p 0+ p1+….+ p
n=0
n+….=1
Or
∑( )p
∞ n
λ
0=1∨ p <1
n=0 μ 0= ∞
1
∑ ( μλ )
n
;
λ
μ
n=0
The denominator of this expression is an infinite geometric series, whose sum is:
∞ n
∑ ( μλ )=
1
1−( λ/ μ)
n=0
p
And n= ( λμ ) (1− λμ )= p (1− p ); p<1 ,n =0 , 1 ,2 … …
n
n
This expression gives the required probability distribution for exactly n customers in the
queuing system. Hence p0 denotes the probability of the system being empty (no
customer).
Step 4 obtain probability density function of waiting time excluding service time
distribution
The waiting time distribution of each customer in the steady-state is same, and it is a
continuous random variable except that there is a non-zero probability that the delay will
be zero probability that the delay will be zero, i.e. waiting time is zero. Let w be the time
required by the server to serve all the customers present in the system at a particular time
in the steady state. Let ∅ w (t) be the probability distribution function of w, i.e.
∅ w ( t )= p( w ≤t ), 0≤ t ≤ ∞ .
If an arriving customer find n(≥1)customers already in the system, then in order for a
customer to get service at time between 0 and t, all the customers must have been
serviced by time t. let s1 , s2…………… sn denote service of n customers, respectively.
Thus
{
0 ; n=0
n
∑ s i ;n ≥1
i=1
The distribution function of waiting time w for a customer who has to wait is given by
{{
p 0=1− p ;n=0 ,t =0
p { w ≤ t }= n
p ∑ s i ≤ t }; n≥ 1 ,t >0
i=1
Since the service time for each customer is independent and identically distributed,
therefore its probability function density is φ s (t )=μe −μt
, t> 0 , where μ the mean service rate is.
Thus
n
∅ n ( t ) =∑ [ pn prob { ( n−1 ) customers got service at time t } ]
i=1
( )
n
λ λn
¿ ∑ 1− ( )¿
i=1 μ μ
{
∞ 1− p;t =0
∅ w (t)= p { w ≤ t }= ∑ t
n=1 pn ∫ ∅n ( t ) dt
0
¿¿
¿¿
{
1−p ;t=0
¿ t
p (1− p)∫ μe
− μ ( 1− p) t dt
This shows that the waiting time distribution is discontinuous at t=0 and continuous in
the range 0< t<∞ . thus expression for ∅ w (t) may also be written as
∅ w ( 0 )=1− p ; t=0
( λμ ) e
¿ λ 1− − ( μ−λ ) t
dt ; t >0
For the busy period distribution, let the random variable w denote the total time (waiting
and service) that a customer spends in the system. Then the probability density function
for the distribution is given by
∅ (w ) d
∅ ( w : w>0 ) = ; ∅ (w)= { ∅ w (t ) }
prob ( w> 0 ) di
λ −(μ− λ)t
λ(1− )e
∅ (w) μ
¿ ∞
=∞
∫ ∅ ( t ) dt ∫ λ(1− μλ ) e−(μ− λ ) t dt
0 o
λ
λ(1− )e−(μ− λ)t
μ
=( μ−λ ) e−( μ− λ ) t ; t=0
λ
μ
∞ ∞
∫ ∅ ( w :w >0)dt =∫ ( μ−λ ) e ❑
−( μ−λ ) tdt
=1
Now
0 0
1. a) Expected number of customers in the system (customers in the line plus the customer
being served)
∞ ∞
Ls =∑ nPn=∑ n (1− p ) p n , 0< p< 1
n=0 n=0
∞ ∞
¿ ( 1− p ) ∑ np = p ( 1− p ) ∑ np n−1
n
n=0 n=1
¿ p ( 1− p ) { 1+2 p+ 3 p + … } 2
¿ ( 1− p )
{ p
( 1−p )2 }[∑ of an artithmatico−geometric series ]
0r Ls= p
=
λ
; p=
λ
1− p μ −λ μ
L ∞ ∞ ∞
q=¿ ∑ (n−1 ) pn=∑ nPn− ∑ p n¿
n=1 n=1 n=1
∞
¿ ∑ nPn−¿ ¿
n=0
λ λ λ
Ls – ( 1− p0 )= − ; 1− p0 =
λ−μ μ μ
2
λ
Lq =
μ( μ−λ)
∞
¿ ∫ t . λ (1− p)e
−μ ( 1− p ) t
dt
0
( λμ ) (μ−λ1 )
¿ λ 1− 2
λ L
Or W q= = q
μ (μ− λ) λ
b) Expected waiting time for a customer in the system (waiting and service)
1 λ 1
¿ W q+ = +
μ μ(μ−λ) μ
1 Ls
0r ¿ =
μ−λ λ
3. probability that the number of customers in the system is greater than or equal to k
p ( n ≥ k )=¿
4. the variance (fluctuation) of queue length
∞
Var ( n )=∑ n2 p n – ¿
n=1
∞
¿ ∑ n 2 pn−¿
n=1
¿ ( 1− p ) [ 1. p + 22 . p2+ 32 . p3 +… . ] −¿
2
p
Var ( n )=
¿¿
5. probability that the queue is non-empty
p ( n>1 )=1− p0− p1
¿ 1− 1− ( λμ )−(1− μλ )( μλ )=¿
6. probability of k or more customers in the system
∞ ∞
p ( n ≥ k )=∑ p k =∑ (1− p) p k
n=k n=k
∞
¿ ( 1− p ) pk ∑ pn−k
n=k
(1− p) pk
¿ ( 1− p ) pk [ 1+ p + p2 +… ]= = pk
(1− p)
Or p ( n ≥ k )=¿
¿¿
Example
A television repairman finds that the time spent on his jobs has an exponential
distribution with a mean of 30 minutes. If he repairs sets in the order in which they come
in, and if the arrival of sets follows a Poisson distribution approximately with an average
rate of 10 per 8-hour day, what is the repairman’s expected idle time each day? How
many jobs are headed of the average set just brought in?
Solution
Number of hours for which the repairman remains busy in an 8 hour day (traffic
intensity) is given by
Example 2
a) What is the probability that a person arriving at the booth will have to wait?
b) The telephone department will install a second booth when convinced that an arrival
would expect waiting for at least 3 minutes for a phone call. By how much should the
flow of arrivals increase in order to justify a second booth?
c) What is the average length of queue that forms from time to time?
d) What is the probability that it will take him more than 10 minutes altogether to wait for
the phone and complete his call?
Solution
λ=¿1/10= 0.10 person per minute and μ =1/3=0.33 person per minute
Where W q=¿ given) and λ=λ ' (say) for the second booth
Hence, the increase in the arrival rate is 0.60-0.10=0.06 arrivals per minute
0 μ
[[ ]
∞ −0.23 t
e
¿ ∫ ( 0.3 ) ( 0.23 ) e
0.23 t
dt=0.069
10 −0.23
This shows that 3 percent of the arrivals on average will have to wait for 10 minutes or
more before they can use the phone.
Model II ¿
This model is identical to the model I with a difference only in queue discipline. Since
the derivation of Pn is independent of any specific queue discipline, therefore in this
model also we have
n
Pn= (1− p ) p ; n=1 ,2 , … … .
Consequently, other results will also remain unchanged in any queuing system as long as
Pn remain unchanged.
Model III{ (M /M /1):(N / FCFS) } exponential service, Finite (or limited) Queue
The model is different from model I in respect of the capacity of the system. Suppose that
no more than N customers can be accommodated at any time in the system due to certain
reasons. Thus any customer arriving when the system is already contains N customers
dies not enter the system and is lost.
The difference equations derived in model I will also be same for this model as long as
n< N . The systems of steady-state difference equations for this model are:
λ p0 =μ p1 ; n=0
λ p N−1=μ p N ; n=N
Assumptions of this model are same as that of Model I except that the length of the queue
is limited. Moreover, in this case the service rate does not have to exceed arrival rate (
μ> λ ) in order to obtain steady state conditions.
Using the usual procedure from first two difference equations the probability of a
customer in the system for n= 0, 1. 2 ….N are obtained as follows.
Pn=¿
N
¿ p0 ∑ ¿ ¿
n=0
¿ p0 [ 1+ p + p 2+ … p N ]= p
[ ]
N +1
1− p
0
1− p
1− p
And consequently, p0= N +1
; p ≠1∧p <1; p=λ /μ
1− p
{( 1−1−p p ) p ; n ≤ N ; p ≠ 1(λ ≠ μ)
n
N +1
1
; p=1( λ=μ)
N +1
The steady-state solution in this case exists even for p>1. this is due to the limited
capacity of the system which controls the arrivals by the queue length ( ¿ N−1 ) not by the
same as in model I.
{
N+ 1
p (N +1) p
− ; p ≠ 1( λ ≠ μ)
1−p 1− p N +1
N
; p=1 (λ=μ)
2
Example
Consider a single server queuing system with Poisson input, exponential service times.
Suppose the mean arrival rate is 3 calling units per hour, the expected service time is 0.25
hour and the maximum permissible calling units in the system is two. Derive the steady-
state probability distribution of the distribution of the number of calling units in the
system, and then calculate the expected number in the system.
Solution
¿ ( 1−0.75 ) ¿ ¿
p (1− p) 1−0.75
0= N +1
= ¿
1− p 1−¿¿
2
¿ 0.43 ∑ n ¿ ¿
n=1
Chapter summary
The queuing theory also called the waiting line theory, owes its development to A K
Erlang’s effort to analyze telephone traffic congestion with a view to satisfying the
randomly arising demand for the services of the Copenhagen automatic telephone system,
in year 1909.
It is assumed that customers joining a queuing system arrive in random manner and
follow a Poisson distribution or equivalently the inter-arrival times follow exponential
distribution.
Review Question
2. Customers arrive at a one window drive according to a poison distribution with mean of 10
minutes and service time per customer is Solve exponential with mean of 6 minutes. The
space in front of the window can accommodate only 3 vehicles including the serviced
one. Other vehicles have to wait outside this space. Calculate
a. probably that an arriving customer can drive directly to the space in front of the window
b. Probability that an arriving customer will have to wait outside the directed space.
c. how long an arriving customer is expected to wait before gating the service
3. A branch of Dashen Bank has only one typist. Since the typing work varies in length
(number of pages to the typed) the typing rate is randomly distributed approximating a
poison distribution with mean service rate of 8 letters per hour. The letter arrives at a rate
of 5 per hour during the entire 8 hour work day. If the type writer is values at birr 1.50
per hour, determine
a. Equipment utilization
b. Average system time
c. The percent time that an arriving letter has to wait
d. Average cost due to waiting on the part of type writer, i.e. it remaining idle.
Chapter one
1. A, 2. A, 3. D, 4. A, 5. A
Chapter 2
1. A 7. D 13. B
2. D 8. F 14. A
3. F 9. A 15. A
4. D 10. B 16. B
5. B 11. B 17. A
6. D 12. C 18. D
Part II
1. Answer: X1=4, X2 =2 2. Answer:X1=1,X2=5
and Maximize Z=230 and Minimize Z=13
3. Answer: 10 units of model X, 2 units of model Y and the maximum profit is $
560.
4. Answer: X1=187.5, X2 =125 and Maximize Z=21,875
5. Minimize 18Y1 + 14Y2
Subject to 2Y1 + 2Y2 ≥ 12
3Y1 + Y2 ≥ 10
Y1 >= 0, Y2 ≥ 0
Chapter three
Initial feasible solution
W1 W2 W3 W4 Factory Capacity
F1 3 2 7 6 5000
Factory 5000
F2 7 5 2 3 6000
F3 2 5 4 5 2500
1000 1500
F1 W1 5000 3 $ 15000
F2 W1 1000 7 7000
F2 W2 4000 5 20000
F2 W3 1000 2 2000
F3 W3 1000 4 4000
F3 1500 5 7500
W4
Total transportation cost =$55,500
Least- Cost method is better than the NWCM because it considers cost factories
3.
A e
B c
C b
D d
Chapter 4
1. D 3. B 5. A
2. C 4. D
6. Answer b, regular share 0.047, b. risky shares 0.050 and property 0.060. c.
since the expected of 6% highest for property, the investor should invest in
this alternative.
Chapter five
1. A 3. A 5. D
2. A 4. B
Chapter six