Powerco has three electric power plants that supply the needs of four cities. The numbers of kilowatt-hours (kwh) of electricity that can be supplied by power plants 1, 2, and 3 are 35 million, 50 million, and 40 million, respectively. The peak power demands in the cities, which occur at the same time of day (2 pm), are as follows (in kwh).
City 1: 45 million; City 2: 20 million; City 3: 30 million; City 4: 30 million. The costs (in dollars) of sending one million kwh of electricity from each plant to each city depend on the distance the electricity must travel. These costs are shown in the following table.
To:
City 1
City 2
City 3
City 4
From:
Plant 1
8
6
10
9
Plant 2
9
12
13
7
Plant 3
14
9
16
5
Powerco wants to determine how much electricity to send from each plant to each city to minimize the total cost
of meeting each city's peak power demand.
a.
Use the northwest corner method to obtain an initial feasible solution to Powerco's transportation problem.
b.
Starting with the initial feasible solution from part a. find an optimal solution to Powerco's problem.
c.
If the cost of sending one million kwh of electricity from Plant 1 to City 1 changes from $8 to $12,
does the optimal solution found in part a. change? Explain.
2.
a. State the Complementary Slackness Theorem for a pair of linear programming
problems of the following form.
P.minimizectx
subjectto Ax
³ b
x
³ 0
Q.maximizebty
subjectto Aty
£ c
y
³ 0
where A is an m×n matrix.
b.
Use the Complementary Slackness Theorem to prove or disprove that (5,0,1) is an optimal solution of the following linear programming problem.
minimize 2x1+3x2+4x3
subjectto x1- x2+3x3
³ 5
2x1+ x2- x3
³ 4
x1+2x2+ x3
³ 6
xj
³ 0, j=1,2,3
3.
Let problems P and Q be as follows:
P.minimizectx
subjectto Ax
³ b
x
³ 0
Q.maximizebty
subjectto Aty
£ c
y
³ 0
where A is an m×n matrix.
a.
Prove: If [`(x)] and [`(y)] are feasible vectors for P and Q, respectively, then ct[`(x)] ³ bt[`(y)].
b.
Consider the following problem.
Q:maximize b1y1-4y2+ y3
subjectto y1- y2+ y3
£ 1
-2y1+ y2+2y3
£ 1
yj
³ 0, j=1,2,3
For what values of b1 is the objective function of this problem unbounded above over the feasible region? For these values of b1, what conclusion can you draw about the dual of problem Q?
Justify your answers.
4.
Consider the following problem.
P:minimize z=6x1+2x2+x3-3x4
subjectto x1-2x2- x3+3x4
£ 2
2x1-4x2- x3+3x4
=3
3x1-6x2+ x3+6x4
=6
xj
³ 0, j=1,2,3,4
a.
Solve problem P using the two phase simplex method.
b.
Let c2 be the coefficient of x2 in the objective function of problem P, and let b1 be the number on the right side of the first constraint. Use sensitivity analysis to determine:
i.) the range of values for q such that replacing c2=2 by c2=2+q does not
change the optimal basis found in part a.
ii.) the range of values for r such that replacing b1=2 by b1=2+r does not
a. Prove that the equation ex=2+x+[(x2)/2] has exactly one real solution.
b.
Let a be the solution in part a. Find an approximation b of a such that
|a-b| < 10-6.
c.
Prove that the inequality in b. holds.
2.
Recall that the truncation errors ETn and ESm in approximating òabf using the (extended) trapezoidal rule with n subintervals and Simpson's rule with m subintervals are given by the following formulas.
ETn
=-
(b-a)3
12n2
f¢¢(xT), some xT Î (a,b), and
ESm
=-
(b-a)5
180m4
fiv(xS), some xS Î (a,b)
Suppose one approximates òp/6p/3sinq dq using the trapezoidal rule and Simpson's rule. Show that Simpson's rule with m=6 yields a more accurate approximation of this integral than the trapezoidal rule with n=100; i.e., for this function |ES6| < |ET100|.
3.
Let A be a n×n matrix of the following form.
A=
æ ç ç ç ç ç ç ç ç
ç ç ç ç ç ç ç è
a1
d1
0
0
0
0
¼
0
b1
a2
d2
0
0
0
¼
0
c1
b2
a3
d3
0
0
¼
0
0
c2
b3
a4
d4
0
¼
0
:
···
:
0
¼
¼
cn-4
bn-3
an-2
dn-2
0
0
¼
¼
0
cn-3
bn-2
an-1
dn-1
0
¼
¼
0
0
cn-2
bn-1
an
ö ÷ ÷ ÷ ÷ ÷ ÷ ÷ ÷
÷ ÷ ÷ ÷ ÷ ÷ ÷ ø
a.
Construct an algorithm which generates the upper-triangular matrix obtained by applying Gaussian elimination to A. That is, show how to construct the upper-triangular matrix U so that A=LU where L is lower-triangular and U has only two significant diagonals; i.e.,
U=
æ ç ç ç ç ç ç
ç ç ç ç ç è
u1
v1
0
0
¼
0
0
u2
v2
0
¼
0
0
0
u3
v3
¼
0
:
···
:
0
0
0
¼
un-1
vn-1
0
0
0
¼
0
un
ö ÷ ÷ ÷ ÷ ÷ ÷
÷ ÷ ÷ ÷ ÷ ø
.
You may assume that pivoting is unnecessary. Exploit the sparsity pattern by writing everything in terms of the given vectors. (Do not use double subscripts.)
b.
Determine the number of additions/subtractions, multiplications, and divisions in your algorithm.
4.
Let
A=
æ ç ç ç
ç ç è
3
-1
0
0
-5
2
0
1
4
ö ÷ ÷ ÷
÷ ÷ ø
.
Suppose Gauss-Seidel is used to solve the linear system Ax=b.
a.
Construct the matrix TGS such that Gauss-Seidel may be described as
x(i+1)=TGSx(i)+c, i=0,1,2,3,¼ .
b.
Determine the minimum k, (0 < k < 1) where it is guaranteed that
According to a survey conducted by the American Bar Association 1 in every 410 Americans is a lawyer, but 1 in every 64 residents of Washington, D.C. is a lawyer.
a.
If you select a random sample of 1500 Americans, what is the approximate probability that the sample contains at least one lawyer?
b.
If you select the same size random sample from among the residents of Washington, D.C., what is the approximate probability that the sample contains at least thirty lawyers?
c.
What is the expected number of lawyers in America from the sample of 1500 persons? What is the expected number of lawyers in Washington, D.C. from the sample of 1500 persons?
d.
If you stand on a Washington, D.C. street corner and interview the first 1000 persons who walk by, and thirty say that they are lawyers, does this suggest that the frequency of lawyers passing the corner exceeds the density within the city? Justify and explain your answer.
2.
Suppose the joint probability mass function (p.m.f.) of X and Y is given by the following table.
\BeginTable \BeginFormat | l | c | c | c | c | \EndFormat " " \use4 x " "y | 0 " 1 " 2 " 3 " +20 _ " 1 | 1/8 " 1/16 " 3/16 " 1/8 " " 2 | 1/16 " 1/16 " 1/8 " 1/4 " \EndTable
a.
Determine the marginal p.m.f. of X, fX(x), and the marginal p.m.f. of Y, fY(y).
Sketch these p.m.f.s and describe their shapes.
b.
Calculate E[X], E[Y], and E[X | Y=y], for y=1,2.
c.
Calculate E[E[X | Y]]. Compare E[X] and E[E[X | Y]].
d.
Calculate Var(X) and Var(Y).
3.
A test correctly identifies a disease D with probability 0.95 and wrongly diagnoses D with probability 0.01. From past experience it is known that the disease occurs in a targeted population with frequency 0.2%. An individual is chosen at random from this population and is given the test. Calculate the probability that
a.
the test is positive, P(+).
b.
the individual actually suffers from the disease D if the test turns out to be positive, P(D|+).
c.
the individual actually does not suffer from the disease Dc if the test turns out to be negative, P(Dc|+).
4.
If X1,X2 are independent Gamma(ai,1), i=1,2, where
f(xi)=
1
G(ai)
xiai-1e-xi 0 < xi < ¥
determine the following.
a.
Determine the joint density of X1 and X2, f(x1,x2).
b.
Let
Y1=
X1
X1+X2
and let
Y2=X1+X2.
Determine the joint density of Y1 and Y2, f(y1,y2).
c.
Calculate the marginal densities of Y1, Y2, fY1(y1), and fY2(y2).
d.
Is Y2 independent of Y1?
(Note: The distribution of Y1 is called the Beta distribution.)