Visit https://testbankdeal.

com to download the full version and

explore more testbank or solution manual

Applied Numerical Methods WMATLAB for Engineers

and Scientists 3rd Edition Chapra Solutions
Manual

_____ Click the link below to download _____

https://testbankdeal.com/product/applied-numerical-
methods-wmatlab-for-engineers-and-scientists-3rd-
edition-chapra-solutions-manual/

Explore and download more testbank at testbankdeal.com

Here are some suggested products you might be interested in.
Click the link to download

Applied Numerical Methods with MATLAB for Engineers and

Scientists 2nd Edition Steven Chapra Solutions Manual

https://testbankdeal.com/product/applied-numerical-methods-with-
matlab-for-engineers-and-scientists-2nd-edition-steven-chapra-
solutions-manual/

Numerical Methods for Engineers 7th Edition Chapra

Solutions Manual

https://testbankdeal.com/product/numerical-methods-for-engineers-7th-
edition-chapra-solutions-manual/

Numerical Methods for Engineers 6th Edition Chapra

Solutions Manual

https://testbankdeal.com/product/numerical-methods-for-engineers-6th-
edition-chapra-solutions-manual/

Digital Control System Analysis and Design 4th Edition

Phillips Solutions Manual

https://testbankdeal.com/product/digital-control-system-analysis-and-
design-4th-edition-phillips-solutions-manual/
Network Management Principles and Practices 2nd Edition
Subramanian Solutions Manual

https://testbankdeal.com/product/network-management-principles-and-
practices-2nd-edition-subramanian-solutions-manual/

Managerial Economics 3rd Edition Froeb Solutions Manual

https://testbankdeal.com/product/managerial-economics-3rd-edition-
froeb-solutions-manual/

Strategic Management 3rd Edition Rothaermel Test Bank

https://testbankdeal.com/product/strategic-management-3rd-edition-
rothaermel-test-bank/

Social Statistics for a Diverse Society 8th Edition

Frankfort Nachmias Test Bank

https://testbankdeal.com/product/social-statistics-for-a-diverse-
society-8th-edition-frankfort-nachmias-test-bank/

Introduction to Audiology 13th Edition Martin Solutions

Manual

https://testbankdeal.com/product/introduction-to-audiology-13th-
edition-martin-solutions-manual/
Contemporary Business Mathematics Canadian 10th Edition
Hummelbrunner Solutions Manual

https://testbankdeal.com/product/contemporary-business-mathematics-
canadian-10th-edition-hummelbrunner-solutions-manual/
 1

CHAPTER 9
9.1 The flop counts for the tridiagonal algorithm in Fig. 9.6 can be summarized as

Mult/Div Add/Subtr Total

Forward elimination 3(n – 1) 2(n – 1) 5(n – 1)
Back substitution 2n – 1 n–1 3n – 2
Total 5n – 4 3n – 3 8n – 7

Thus, as n increases, the effort is much, much less than for a full matrix solved with Gauss elimination
which is proportional to n3.

9.2 The equations can be expressed in a format that is compatible with graphing x2 versus x1:

18  2 x1
x2   3  0.333333 x1
6
40  x1
x2   5  0.125 x1
8

which can be plotted as

a11=2;a12=-6;b1=-18;
a21=-1;a22=8;b2=40;
x1=[0:20];x21=(b1-a11*x1)/a12;x22=(b2-a21*x1)/a22;
plot(x1,x21,x1,x22,'--'),grid
xlabel('x_1'),ylabel('x_2')

10

8
x2

2
0 5 10 15 20
x1

Thus, the solution is x1 = 9.6, x2 = 6.2. The solution can be checked by substituting it back into the
equations to give

2(9.6)  6(6.2)  19.2  37.2  18

9.6  8(6.2)  9.6  49.6  40

9.3 (a) The equations can be rearranged into a format for plotting x2 versus x1:

x2  14.25  0.77 x1

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 2

1.2
x2  11.76471  x1
1.7
20

0
0 20 40 60 80
-20

-40

-60

If you zoom in, it appears that there is a root at about (38.7, –15.6).
38 38.5 39 39.5

-15

-15.4

-15.8

-16.2

The results can be checked by substituting them back into the original equations:

0.77(38.7)  15.6  14.2  14.25

1.2(38.7)  1.7(15.6)  19.92  20

(b) The plot suggests that the system may be ill-conditioned because the slopes are so similar.

(c) The determinant can be computed as

D  0.77(1.7)  1(1.2)  0.11

which is relatively small. Note that if the system is normalized first by dividing each equation by the largest
coefficient,

0.77 x1  x2  14.25
0.705882 x1  x2  11.76471

the determinant is even smaller

D  0.77(1)  1(0.705882)  0.064

9.4 (a) The determinant can be computed as:

A1  1 1  1(0)  1(1)  1
1 0

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 3

A2  2 1  2(0)  1(3)  3
3 0

A3  2 1  2(1)  1(3)  1
3 1
D  0(1)  2(3)  5(1)  1

(b) Cramer’s rule

1 2 5
1 1 1
2 1 0 2
x1    2
D 1
0 1 5
2 1 1
3 2 0 8
x2   8
D 1
0 2 1
2 1 1
3 1 2 3
x3    3
D 1

(c) Pivoting is necessary, so switch the first and third rows,

3x1  x2 2
2 x1  x2  x3  1
2 x2  5 x3  1

Multiply pivot row 1 by 2/3 and subtract the result from the second row to eliminate the a21 term. Note that
because a31 = 0, it does not have to be eliminated

3x1  x2 2
0.33333x2  x3  0.33333
2 x2  5 x3  1

Pivoting is necessary so switch the second and third row,

3x1  x2 2
2 x2  5 x3  1
0.33333x2  x3  0.33333

Multiply pivot row 2 by 0.33333/2 and subtract the result from the third row to eliminate the a32 term.

3x1  x2 2
2 x2  5 x3  1
 0.16667 x3  0.5

Note that, at this point, the determinant can be computed as the product of the diagonal elements

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 4

D  3  2  0.16667  ( 1) 2  1

The solution can then be obtained by back substitution

0.5
x3   3
0.16667
1  5(3)
x2  8
2
2  0(3)  1(8)
x1   2
3

The results can be checked by substituting them back into the original equations:

2(8)  5(3)  1
2(2)  8  3  1
3(2)  8  2

9.5 (a) The equations can be expressed in a format that is compatible with graphing x2 versus x1:

x2  0.5 x1  9.5
x2  0.51x1  9.4

The resulting plot indicates that the intersection of the lines is difficult to detect:

22

20

18

16

14

12

10
5 10 15 20

Only when the plot is zoomed is it at all possible to discern that solution seems to lie at about x1 = 14.5 and
x2 = 10.

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 5

14.7
14.65
14.6
14.55
14.5
14.45
14.4
14.35
14.3
9.75 10 10.25

(b) The determinant can be computed as

0.5 1
 0.5(2)  (1)(1.02)  0.02
1.02 2

which is close to zero.

(c) Because the lines have very similar slopes and the determinant is so small, you would expect that the
system would be ill-conditioned

(d) Multiply the first equation by 1.02/0.5 and subtract the result from the second equation to eliminate the
x1 term from the second equation,

0.5 x1  x2  9.5
0.04 x2  0.58

The second equation can be solved for

0.58
x2   14.5
0.04

This result can be substituted into the first equation which can be solved for

9.5  14.5
x1   10
0.5

(e) Multiply the first equation by 1.02/0.52 and subtract the result from the second equation to eliminate the
x1 term from the second equation,

0.52 x1  x2  9.5
 0.03846 x2  0.16538

The second equation can be solved for

0.16538
x2   4.3
0.03846

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 6

This result can be substituted into the first equation which can be solved for

9.5  4.3
x1   10
0.52

Interpretation: The fact that a slight change in one of the coefficients results in a radically different solution
illustrates that this system is very ill-conditioned.

9.6 (a) Multiply the first equation by –3/10 and subtract the result from the second equation to eliminate the
x1 term from the second equation. Then, multiply the first equation by 1/10 and subtract the result from the
third equation to eliminate the x1 term from the third equation.

10 x1  2 x2  x3  27
 4.4 x2  1.7 x3  53.4
0.8 x2  6.1x3  24.2

Multiply the second equation by 0.8/(–4.4) and subtract the result from the third equation to eliminate the x2
term from the third equation,

10 x1  2 x2  x3  27
 4.4 x2  1.7 x3  53.4
6.409091x3  33.9091

Back substitution can then be used to determine the unknowns

33.9091
x3   5.29078
6.409091
(53.4  1.7(5.29078))
x2   10.0922
4.4
(27  5.29078  2(10.0922))
x1   0.152482
10

(b) Check:

10(0.152482)  2(10.0922)  (5.29078)  27

3(0.152482)  5(10.0922)  2(5.29078)  61.5
0.152482  10.0922  5( 5.29078)  21.5

9.7 (a) Pivoting is necessary, so switch the first and third rows,

8 x1  x2  2 x3  20
3x1  x2  7 x3  34
2 x1  6 x2  x3  38

Multiply the first equation by –3/(–8) and subtract the result from the second equation to eliminate the a21
term from the second equation. Then, multiply the first equation by 2/(–8) and subtract the result from the
third equation to eliminate the a31 term from the third equation.

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 7

8 x1  x2  2 x3  20
 1.375 x2  7.75 x3  26.5
 5.75 x2  1.5 x3  43

Pivoting is necessary so switch the second and third row,

8 x1  x2  2 x3  20
 5.75 x2  1.5 x3  43
 1.375 x2  7.75 x3  26.5

Multiply pivot row 2 by –1.375/(–5.75) and subtract the result from the third row to eliminate the a32 term.

8 x1  x2  2 x3  20
 5.75 x2  1.5 x3  43
8.108696x3  16.21739

At this point, the determinant can be computed as

D  8  5.75  8.108696  ( 1) 2  373

The solution can then be obtained by back substitution

16.21739
x3   2
8.108696
43  1.5( 2)
x2  8
5.75
20  2(2)  1(8)
x1  4
8

(b) Check:

2(4)  6(8)  (2)  38

3(4)  (8)  7(2)  34
8(4)  (8)  2(2)  20

9.8 Multiply the first equation by –0.4/0.8 and subtract the result from the second equation to eliminate the
x1 term from the second equation.

 0.8 0.4   x1   41 
    x   45.5
 0.6 0.4  2  
 0.4 0.8   x3   105 

Multiply pivot row 2 by –0.4/0.6 and subtract the result from the third row to eliminate the x2 term.

 0.8 0.4   x1   41 
    
 0.6 0.4   x2    45.5 
 0.533333  x3  135.3333

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 8

The solution can then be obtained by back substitution

135.3333
x3   253.75
0.533333
45.5  (0.4)253.75
x2   245
0.6
41  (0.4)245
x1   173.75
0.8

(b) Check:

0.8(173.75)  0.4(245)  41
0.4(173.75)  0.8(245)  0.4(253.75)  25
0.4(245)  0.8(253.75)  105

9.9 Mass balances can be written for each of the reactors as

200  Q13c1  Q12 c1  Q21c2  0

Q12 c1  Q21c2  Q23c2  0
500  Q13c1  Q23c2  Q33c3  0

Values for the flows can be substituted and the system of equations can be written in matrix form as

130 30 0   c1  200 

 90 90    
 0  c2    0 
   
 40 60 120   c3  500 

The solution can then be developed using MATLAB,

>> A=[130 -30 0;-90 90 0;-40 -60 120];

>> B=[200;0;500];
>> C=A\B

C =
2.0000
2.0000
5.8333

9.10 Let xi = the volume taken from pit i. Therefore, the following system of equations must hold

0.52 x1  0.20 x2  0.25 x3  4800

0.30 x1  0.50 x2  0.20 x3  5800
0.18 x1  0.30 x2  0.55 x3  5700

MATLAB can be used to solve this system of equations for

>> A=[0.55 0.2 0.25;0.3 0.5 0.2;0.18 0.3 0.55];

>> b=[4800;5800;5700];
>> x=A\b

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 9

x =
1.0e+003 *
4.0058
7.1314
5.1628

Therefore, we take x1 = 4005.8, x2 = 7131.4, and x3 = 5162.8 m3 from pits 1, 2 and 3 respectively.

9.11 Let ci = component i. Therefore, the following system of equations must hold

15c1  17c2  19c3  2120

0.25c1  0.33c2  0.42c3  43.4
1.0c1  1.2c2  1.6c3  164

The solution can be developed with MATLAB:

A=[15 17 19;0.25 0.33 0.42;1 1.2 1.6];

b=[2120;43.4;164];
c=A\b

c =
20.0000
40.0000
60.0000

Therefore, c1 = 20, c2 = 40, and c3 = 60.

9.12 Centered differences (recall Chap. 4) can be substituted for the derivatives to give

ci 1  2ci  ci 1 c c
0D  U i 1 i 1  kci
x 2 2x

collecting terms yields

– (D + 0.5Ux)ci–1 + (2D + kx2)ci – (D – 0.5Ux)ci+1 = 0

Assuming x = 1 and substituting the parameters gives

2.5ci 1  4.2ci  1.5ci 1  0

For the first interior node (i = 1),

4.2c1  1.5c2  200

For the last interior node (i = 9)

2.5c8  4.2c9  15

These and the equations for the other interior nodes can be assembled in matrix form as

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 10

 4.2 1.5 0 0 0 0 0 0 0   c1   200 

 2.5 4.2 1.5
 0 0 0 0 0 0  c2   0 
 0 2.5 4.2 1.5 0 0 0 0 0   c3   0 
    
 0 0 2.5 4.2 1.5 0 0 0 0  c4   0 
 0    
0 0 2.5 4.2 1.5 0 0 0  c5    0 
   
 0 0 0 0 2.5 4.2 1.5 0 0  c6 0 
   
 0 0 0 0 0 2.5 4.2 1.5 0  c7   0 
    
 0 0 0 0 0 0 2.5 4.2 1.5 c8   0 
 
 0 0 0 0 0 0 0 2.5 4.2  c9   15 

The following script generates the solution with the Tridiag function from p. 247 and develops a plot:

clear,clc,clf
D=2;U=1;k=0.2;c0=80;c10=10;dx=1;
diag=(2*D+k*dx^2);
super=-(D-0.5*U*dx);
sub=-(D+0.5*U*dx);
r1=-sub*c0; rn=-super*c10;
n=9;
e=ones(n,1)*sub;f=ones(n,1)*diag;g=ones(n,1)*super;
r=zeros(n,1);r(1)=r1;r(n)=rn;
c=Tridiag(e,f,g,r)
c=[80 c 10]; x=0:1:10;
plot(x,c) ylim([0 90])

Alternatively, as in the following script, the solution can be generated directly with MATLAB left division:

clear,clc,clf
A=[4.2 -1.5 0 0 0 0 0 0 0
-2.5 4.2 -1.5 0 0 0 0 0 0
0 -2.5 4.2 -1.5 0 0 0 0 0
0 0 -2.5 4.2 -1.5 0 0 0 0
0 0 0 -2.5 4.2 -1.5 0 0 0
0 0 0 0 -2.5 4.2 -1.5 0 0
0 0 0 0 0 -2.5 4.2 -1.5 0
0 0 0 0 0 0 -2.5 4.2 -1.5
0 0 0 0 0 0 0 -2.5 4.2];
b=[200 0 0 0 0 0 0 0 15]';
c=(A\b)'
c=[80 c 10]; x=0:1:10;
plot(x,c), ylim([0 90])

In either case, the results are:

c =
68.6613 58.9183 50.5357 43.3029 37.0219 31.4898 26.4684 21.6284 16.4454

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
Visit https://testbankdead.com
now to explore a rich
collection of testbank,
solution manual and enjoy
exciting offers!
 11

80

60

40

20

0
0 2 4 6 8 10

9.13 For the first stage, the mass balance can be written as

F1 yin  F2 x2  F2 x1  F1 x1

Substituting x = Ky and rearranging gives

 F  F
  1  2 K  y1  2 Ky2   yin
 F1  F1

Using a similar approach, the equation for the last stage is

 F  F
y4   1  2 K  y5   2 xin
 F1  F1

For interior stages,

 F  F
yi-1   1  2 K  yi  2 Kyi 1  0
 F1  F1

These equations can be used to develop the following system,

11 10 0 0 0   y1  0.1

 1 11 10 0 0   y2   0 
 0 1 11 10 0   y3    0 
0 0 1 11 10   y4   0 
 0 0 0 1 11   y5   0 

The solution can be developed in a number of ways. For example, using MATLAB,

>> format short g

>> A=[11 -10 0 0 0;-1 11 -10 0 0;0 -1 11 -10 0;0 0 -1 11 -10;0 0 0 -1 11];
>> B=[0.1;0;0;0;0];
>> Y=A\B
Y =
0.0099999
0.0009999
9.99e-005

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 12

9.9e-006
9e-007

Note that the corresponding values of X can be computed as

>> X=5*Y
X =
0.05
0.0049995
0.0004995
4.95e-005
4.5e-006

Therefore, yout = 0.0000009 and xout = 0.05.

9.14 Assuming a unit flow for Q1, the simultaneous equations can be written in matrix form as

 2 1 2 00  Q2  0 
0
0
 0 2 10  Q3  0 
2
0 0 0 03  Q4  0 
2
     
1 1 0 00  Q5  1 
0
0 1 1 1 0 0  Q6  0 
    
 0 0 0 1 1 1 Q7  0 

These equations can then be solved with MATLAB,

>> A=[-2 1 2 0 0 0;
0 0 -2 1 2 0;
0 0 0 0 -2 3;
1 1 0 0 0 0;
0 1 -1 -1 0 0;
0 0 0 1 -1 -1];
>> B=[0 0 0 1 0 0 ]';
>> Q=A\B

Q =
0.5059
0.4941
0.2588
0.2353
0.1412
0.0941

9.15 The solution can be generated with MATLAB,

>> A=[1 0 0 0 0 0 0 0 1 0;
0 0 1 0 0 0 0 1 0 0;
0 1 0 3/5 0 0 0 0 0 0;
-1 0 0 -4/5 0 0 0 0 0 0;
0 -1 0 0 0 0 3/5 0 0 0;
0 0 0 0 -1 0 -4/5 0 0 0;
0 0 -1 -3/5 0 1 0 0 0 0;
0 0 0 4/5 1 0 0 0 0 0;
0 0 0 0 0 -1 -3/5 0 0 0;
0 0 0 0 0 0 4/5 0 0 1];
>> B=[0 0 -74 0 0 24 0 0 0 0]';
>> x=A\B

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 13

x =
37.3333
-46.0000
74.0000
-46.6667
37.3333
46.0000
-76.6667
-74.0000
-37.3333
61.3333

Therefore, in kN

AB = 37.3333 BC = 46 AD = 74 BD = 46.6667 CD = 37.3333

DE = 46 CE = 76.6667 Ax = 74 Ay = 37.33333 Ey = 61.3333

9.16
function x=pentasol(A,b)
% pentasol: pentadiagonal system solver banded system
% x=pentasol(A,b):
% Solve a pentadiagonal system Ax=b
% input:
% A = pentadiagonal matrix
% b = right hand side vector
% output:
% x = solution vector

% Error checks
[m,n]=size(A);
if m~=n,error('Matrix must be square');end
if length(b)~=m,error('Matrix and vector must have the same number of
rows');end
x=zeros(n,1);

% Extract bands
d=[0;0;diag(A,-2)];
e=[0;diag(A,-1)];
f=diag(A);
g=diag(A,1);
h=diag(A,2);
delta=zeros(n,1);
epsilon=zeros(n-1,1);
gamma=zeros(n-2,1);
alpha=zeros(n,1);
c=zeros(n,1);
z=zeros(n,1);

% Decomposition
delta(1)=f(1);
epsilon(1)=g(1)/delta(1);
gamma(1)=h(1)/delta(1);
alpha(2)=e(2);
delta(2)=f(2)-alpha(2)*epsilon(1);
epsilon(2)=(g(2)-alpha(2)*gamma(1))/delta(2);
gamma(2)=h(2)/delta(2);
for k=3:n-2
alpha(k)=e(k)-d(k)*epsilon(k-2);
delta(k)=f(k)-d(k)*gamma(k-2)-alpha(k)*epsilon(k-1);

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 14

epsilon(k)=(g(k)-alpha(k)*gamma(k-1))/delta(k);
gamma(k)=h(k)/delta(k);
end
alpha(n-1)=e(n-1)-d(n-1)*epsilon(n-3);
delta(n-1)=f(n-1)-d(n-1)*gamma(n-3)-alpha(n-1)*epsilon(n-2);
epsilon(n-1)=(g(n-1)-alpha(n-1)*gamma(n-2))/delta(n-1);
alpha(n)=e(n)-d(n)*epsilon(n-2);
delta(n)=f(n)-d(n)*gamma(n-2)-alpha(n)*epsilon(n-1);
% Forward substitution
c(1)=b(1)/delta(1);
c(2)=(b(2)-alpha(2)*c(1))/delta(2);
for k=3:n
c(k)=(b(k)-d(k)*c(k-2)-alpha(k)*c(k-1))/delta(k);
end
% Back substitution
x(n)=c(n);
x(n-1)=c(n-1)-epsilon(n-1)*x(n);
for k=n-2:-1:1
x(k)=c(k)-epsilon(k)*x(k+1)-gamma(k)*x(k+2);
end

A script to test the function can be developed as:

clear,clc
A=[8 -2 -1 0 0;-2 9 -4 -1 0;-1 -3 7 -1 -2;0 -4 -2 12 -5;0 0 -7 -3 15];
b=[5 2 1 1 5]';
x=pentasol(A,b)'

x =
1.0825 1.1759 1.3082 1.1854 1.1809

9.17 Here is the M-file function based on Fig. 9.5 to implement Gauss elimination with partial pivoting

function [x, D] = GaussPivotNew(A, b, tol)

% GaussPivotNew: Gauss elimination pivoting
% [x, D] = GaussPivotNew(A,b,tol): Gauss elimination with pivoting.
% input:
% A = coefficient matrix
% b = right hand side vector
% tol = tolerance for detecting "near zero"
% output:
% x = solution vector
% D = determinant

[m,n]=size(A);
if m~=n, error('Matrix A must be square'); end
nb=n+1;
Aug=[A b];
npiv=0;
% forward elimination
for k = 1:n-1
% partial pivoting
[big,i]=max(abs(Aug(k:n,k)));
ipr=i+k-1;
if ipr~=k
npiv=npiv+1;
Aug([k,ipr],:)=Aug([ipr,k],:);
end
absakk=abs(Aug(k,k));
if abs(Aug(k,k))<=tol
D=0;

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
 15

error('Singular or near singular system')

end
for i = k+1:n
factor=Aug(i,k)/Aug(k,k);
Aug(i,k:nb)=Aug(i,k:nb)-factor*Aug(k,k:nb);
end
end
for i = 1:n
if abs(Aug(i,i))<=tol
D=0;
error('Singular or near singular system')
end
end
% back substitution
x=zeros(n,1);
x(n)=Aug(n,nb)/Aug(n,n);
for i = n-1:-1:1
x(i)=(Aug(i,nb)-Aug(i,i+1:n)*x(i+1:n))/Aug(i,i);
end
D=(-1)^npiv;
for i=1:n
D=D*Aug(i,i);
end

Here is a script to solve Prob. 9.5 for the two cases of tol:

clear; clc; format short g

A=[0.5 -1;1.02 -2];
b=[-9.5;-18.8];
disp('Solution and determinant calculated with GaussPivotNew:')
[x, D] = GaussPivotNew(A,b,1e-5)
disp('Determinant calculated with det:')
D=det(A)

The resulting output is

Solution and determinant calculated with GaussPivotNew:

x =
10
14.5
D =
0.02

Determinant calculated with det:

D =
0.02

PROPRIETARY MATERIAL. © The McGraw-Hill Companies, Inc. All rights reserved. No part of this Manual
may be displayed, reproduced or distributed in any form or by any means, without the prior written permission of the
publisher, or used beyond the limited distribution to teachers and educators permitted by McGraw-Hill for their
individual course preparation. If you are a student using this Manual, you are using it without permission.
Exploring the Variety of Random
Documents with Different Content
 "A good business!" declared Captain Stanley emphatically.

"Aye. Now, as to the Mary—— She's moored, I warrant you, so that

she's as safe as if she was beached. I'll leave you five of my men under the
bo'sun to guard her and her cargo until the court makes the order to have her
fetched to St. John's."

"Five will be enough. I'm not very doubtful of the temper of the men
here. They're cowed, and I think that now Pierre and his fellows are locked
up they'll lose any initiative they ever had. Still, we won't take risks, for the
Mary is a prize of considerable value as she stands."

"That's so. And speaking of prizes, I shall recommend that man of yours
for a good competence. It's impossible to over-estimate the value of his and
Dare's work. My word, Stanley, that boy of yours is a good plucked one!"

Captain Stanley flushed with pleasure and looked in Dare's direction.

Dare had heard his name pronounced and had turned inquiringly. His father
beckoned him to approach.

"Well, Dare, my boy, we've settled up the odds and ends of this business.
It's been the most complete success, thanks to you and Ben. You took risks
that I could never approve of, but the results have been so splendid that I've
had no difficulty in promising Captain McDonnell to overlook that part of
the affair. You did splendidly, my boy, splendidly. But I'll spare your
blushes. Besides, if I'm not mistaken, you'll hear more of this from another
and a higher quarter."

"That's so," interpolated Captain McDonnell. "The Government will

learn of your services, my lad, both through the official report and the
medium of your humble servant. And as you've saved them some thousands
in revenue and gained them a great deal more in seizures, you can count on
them doing the right thing."

"But I don't want them——" began Dare, considerably abashed by the

turn the conversation had taken, though he could not help feeling delight in
having earned the praise of his superiors.
 "Of course you don't, boy," interrupted Captain McDonnell, "but that's
neither here nor there. You've been of service, and as it's a Government
affair things must take their proper course. Now, as to the present—— But
you'd better break this to him, Stanley."

Dare looked at his father questioningly. Captain Stanley returned the

look, smiling gravely.

"I've decided, Dare, and Captain McDonnell supports my decision, that

it's best for you to leave Saltern now that our object in coming here has
been attained. The temper of the villagers is uncertain. They're disappointed
and scared, and at such times people are apt to be excessive in their
demonstrations of emotion. It's not that there's any great danger, but they
know of the part you played in the cleaning up of the gang and they don't
feel very friendly towards you, to say the least; and under the circumstances
I'd rather that you left here as soon as possible.

"I'll be following you shortly myself. As soon as the Mary is taken to St.
John's, someone will be sent to relieve me and in time a permanent official
will be appointed. Then we'll do some hunting and fishing in the Humber
Valley. In the meantime I hope you won't mind obliging me by leaving here
alone. I won't order you to go; you've earned the right to decide for
yourself, but I own I'll be considerably relieved if you'll consent to follow
my advice."

Dare flushed.

"Of course, dad," he burst out impulsively, "whatever you wish——"

"But where am I to go?" he asked, when his father had placed his hand
on his shoulder to show his approval.

"Ah! that will interest you, I think. Captain McDonnell has offered to
take you cruising in the Drake for a month."

"Dad!"
 Both Captain Stanley and Captain McDonnell smiled at that enthusiastic,
forceful exclamation.

"Appeals to you, eh?" chuckled McDonnell.

"Rather!" ejaculated Dare. "There's nothing I'd like better, seeing I can't
stay on here."

"Then be on board by five o'clock."

*****

At half-past five the Drake broke out her anchor and, dipping her flag to
the Customs House ensign, slowly got under way. When she reached the
Oven she slackened speed, and a gun was trained on the former harbour of
the smugglers. The shell expelled from it struck the face of the cliff just
above the narrow opening. There came a report as though the cliff itself had
split in twain, then hundreds of tons of loosened rock fell to form a barrier
for all time to the entrance to the cave.

Dare, who was with Captain McDonnell on the bridge, witnessed the
result with considerable satisfaction.

"Well, that's the end of the Oven," he said.

"And a jolly good thing too," said Captain McDonnell. Then he reached
out a hand and rang "Full speed ahead" to the engine-room.

And the Drake, shuddering from stem to stern at the sudden revolutions
of her propeller, leapt forward like a greyhound, and with a white wave at
her prow headed jauntily for the open sea.

The Mayflower Press, Plymouth. William Brendon & Son, Ltd.

 *** END OF THE PROJECT GUTENBERG EBOOK CONTRABAND: A
TALE OF MODERN SMUGGLERS ***

Updated editions will replace the previous one—the old editions will
be renamed.

Creating the works from print editions not protected by U.S.

copyright law means that no one owns a United States copyright in
these works, so the Foundation (and you!) can copy and distribute it
in the United States without permission and without paying
copyright royalties. Special rules, set forth in the General Terms of
Use part of this license, apply to copying and distributing Project
Gutenberg™ electronic works to protect the PROJECT GUTENBERG™
concept and trademark. Project Gutenberg is a registered trademark,
and may not be used if you charge for an eBook, except by following
the terms of the trademark license, including paying royalties for use
of the Project Gutenberg trademark. If you do not charge anything
for copies of this eBook, complying with the trademark license is
very easy. You may use this eBook for nearly any purpose such as
creation of derivative works, reports, performances and research.
Project Gutenberg eBooks may be modified and printed and given
away—you may do practically ANYTHING in the United States with
eBooks not protected by U.S. copyright law. Redistribution is subject
to the trademark license, especially commercial redistribution.

START: FULL LICENSE

THE FULL PROJECT GUTENBERG LICENSE
 PLEASE READ THIS BEFORE YOU DISTRIBUTE OR USE THIS WORK

To protect the Project Gutenberg™ mission of promoting the free

distribution of electronic works, by using or distributing this work (or
any other work associated in any way with the phrase “Project
Gutenberg”), you agree to comply with all the terms of the Full
Project Gutenberg™ License available with this file or online at
www.gutenberg.org/license.

Section 1. General Terms of Use and

Redistributing Project Gutenberg™
electronic works
1.A. By reading or using any part of this Project Gutenberg™
electronic work, you indicate that you have read, understand, agree
to and accept all the terms of this license and intellectual property
(trademark/copyright) agreement. If you do not agree to abide by all
the terms of this agreement, you must cease using and return or
destroy all copies of Project Gutenberg™ electronic works in your
possession. If you paid a fee for obtaining a copy of or access to a
Project Gutenberg™ electronic work and you do not agree to be
bound by the terms of this agreement, you may obtain a refund
from the person or entity to whom you paid the fee as set forth in
paragraph 1.E.8.

1.B. “Project Gutenberg” is a registered trademark. It may only be

used on or associated in any way with an electronic work by people
who agree to be bound by the terms of this agreement. There are a
few things that you can do with most Project Gutenberg™ electronic
works even without complying with the full terms of this agreement.
See paragraph 1.C below. There are a lot of things you can do with
Project Gutenberg™ electronic works if you follow the terms of this
agreement and help preserve free future access to Project
Gutenberg™ electronic works. See paragraph 1.E below.
1.C. The Project Gutenberg Literary Archive Foundation (“the
Foundation” or PGLAF), owns a compilation copyright in the
collection of Project Gutenberg™ electronic works. Nearly all the
individual works in the collection are in the public domain in the
United States. If an individual work is unprotected by copyright law
in the United States and you are located in the United States, we do
not claim a right to prevent you from copying, distributing,
performing, displaying or creating derivative works based on the
work as long as all references to Project Gutenberg are removed. Of
course, we hope that you will support the Project Gutenberg™
mission of promoting free access to electronic works by freely
sharing Project Gutenberg™ works in compliance with the terms of
this agreement for keeping the Project Gutenberg™ name associated
with the work. You can easily comply with the terms of this
agreement by keeping this work in the same format with its attached
full Project Gutenberg™ License when you share it without charge
with others.

1.D. The copyright laws of the place where you are located also
govern what you can do with this work. Copyright laws in most
countries are in a constant state of change. If you are outside the
United States, check the laws of your country in addition to the
terms of this agreement before downloading, copying, displaying,
performing, distributing or creating derivative works based on this
work or any other Project Gutenberg™ work. The Foundation makes
no representations concerning the copyright status of any work in
any country other than the United States.

1.E. Unless you have removed all references to Project Gutenberg:

1.E.1. The following sentence, with active links to, or other

immediate access to, the full Project Gutenberg™ License must
appear prominently whenever any copy of a Project Gutenberg™
work (any work on which the phrase “Project Gutenberg” appears,
or with which the phrase “Project Gutenberg” is associated) is
accessed, displayed, performed, viewed, copied or distributed:
 This eBook is for the use of anyone anywhere in the United
States and most other parts of the world at no cost and with
almost no restrictions whatsoever. You may copy it, give it away
or re-use it under the terms of the Project Gutenberg License
included with this eBook or online at www.gutenberg.org. If you
are not located in the United States, you will have to check the
laws of the country where you are located before using this
eBook.

1.E.2. If an individual Project Gutenberg™ electronic work is derived

from texts not protected by U.S. copyright law (does not contain a
notice indicating that it is posted with permission of the copyright
holder), the work can be copied and distributed to anyone in the
United States without paying any fees or charges. If you are
redistributing or providing access to a work with the phrase “Project
Gutenberg” associated with or appearing on the work, you must
comply either with the requirements of paragraphs 1.E.1 through
1.E.7 or obtain permission for the use of the work and the Project
Gutenberg™ trademark as set forth in paragraphs 1.E.8 or 1.E.9.

1.E.3. If an individual Project Gutenberg™ electronic work is posted

with the permission of the copyright holder, your use and distribution
must comply with both paragraphs 1.E.1 through 1.E.7 and any
additional terms imposed by the copyright holder. Additional terms
will be linked to the Project Gutenberg™ License for all works posted
with the permission of the copyright holder found at the beginning
of this work.

1.E.4. Do not unlink or detach or remove the full Project

Gutenberg™ License terms from this work, or any files containing a
part of this work or any other work associated with Project
Gutenberg™.

1.E.5. Do not copy, display, perform, distribute or redistribute this

electronic work, or any part of this electronic work, without
prominently displaying the sentence set forth in paragraph 1.E.1
with active links or immediate access to the full terms of the Project
Gutenberg™ License.

1.E.6. You may convert to and distribute this work in any binary,
compressed, marked up, nonproprietary or proprietary form,
including any word processing or hypertext form. However, if you
provide access to or distribute copies of a Project Gutenberg™ work
in a format other than “Plain Vanilla ASCII” or other format used in
the official version posted on the official Project Gutenberg™ website
(www.gutenberg.org), you must, at no additional cost, fee or
expense to the user, provide a copy, a means of exporting a copy, or
a means of obtaining a copy upon request, of the work in its original
“Plain Vanilla ASCII” or other form. Any alternate format must
include the full Project Gutenberg™ License as specified in
paragraph 1.E.1.

1.E.7. Do not charge a fee for access to, viewing, displaying,

performing, copying or distributing any Project Gutenberg™ works
unless you comply with paragraph 1.E.8 or 1.E.9.

1.E.8. You may charge a reasonable fee for copies of or providing

access to or distributing Project Gutenberg™ electronic works
provided that:

• You pay a royalty fee of 20% of the gross profits you derive
from the use of Project Gutenberg™ works calculated using the
method you already use to calculate your applicable taxes. The
fee is owed to the owner of the Project Gutenberg™ trademark,
but he has agreed to donate royalties under this paragraph to
the Project Gutenberg Literary Archive Foundation. Royalty
payments must be paid within 60 days following each date on
which you prepare (or are legally required to prepare) your
periodic tax returns. Royalty payments should be clearly marked
as such and sent to the Project Gutenberg Literary Archive
Foundation at the address specified in Section 4, “Information
 about donations to the Project Gutenberg Literary Archive
Foundation.”

• You provide a full refund of any money paid by a user who

notifies you in writing (or by e-mail) within 30 days of receipt
that s/he does not agree to the terms of the full Project
Gutenberg™ License. You must require such a user to return or
destroy all copies of the works possessed in a physical medium
and discontinue all use of and all access to other copies of
Project Gutenberg™ works.

• You provide, in accordance with paragraph 1.F.3, a full refund of

any money paid for a work or a replacement copy, if a defect in
the electronic work is discovered and reported to you within 90
days of receipt of the work.

• You comply with all other terms of this agreement for free
distribution of Project Gutenberg™ works.

1.E.9. If you wish to charge a fee or distribute a Project Gutenberg™

electronic work or group of works on different terms than are set
forth in this agreement, you must obtain permission in writing from
the Project Gutenberg Literary Archive Foundation, the manager of
the Project Gutenberg™ trademark. Contact the Foundation as set
forth in Section 3 below.

1.F.

1.F.1. Project Gutenberg volunteers and employees expend

considerable effort to identify, do copyright research on, transcribe
and proofread works not protected by U.S. copyright law in creating
the Project Gutenberg™ collection. Despite these efforts, Project
Gutenberg™ electronic works, and the medium on which they may
be stored, may contain “Defects,” such as, but not limited to,
incomplete, inaccurate or corrupt data, transcription errors, a
copyright or other intellectual property infringement, a defective or
damaged disk or other medium, a computer virus, or computer
codes that damage or cannot be read by your equipment.

1.F.2. LIMITED WARRANTY, DISCLAIMER OF DAMAGES - Except for

the “Right of Replacement or Refund” described in paragraph 1.F.3,
the Project Gutenberg Literary Archive Foundation, the owner of the
Project Gutenberg™ trademark, and any other party distributing a
Project Gutenberg™ electronic work under this agreement, disclaim
all liability to you for damages, costs and expenses, including legal
fees. YOU AGREE THAT YOU HAVE NO REMEDIES FOR
NEGLIGENCE, STRICT LIABILITY, BREACH OF WARRANTY OR
BREACH OF CONTRACT EXCEPT THOSE PROVIDED IN PARAGRAPH
1.F.3. YOU AGREE THAT THE FOUNDATION, THE TRADEMARK
OWNER, AND ANY DISTRIBUTOR UNDER THIS AGREEMENT WILL
NOT BE LIABLE TO YOU FOR ACTUAL, DIRECT, INDIRECT,
CONSEQUENTIAL, PUNITIVE OR INCIDENTAL DAMAGES EVEN IF
YOU GIVE NOTICE OF THE POSSIBILITY OF SUCH DAMAGE.

1.F.3. LIMITED RIGHT OF REPLACEMENT OR REFUND - If you

discover a defect in this electronic work within 90 days of receiving
it, you can receive a refund of the money (if any) you paid for it by
sending a written explanation to the person you received the work
from. If you received the work on a physical medium, you must
return the medium with your written explanation. The person or
entity that provided you with the defective work may elect to provide
a replacement copy in lieu of a refund. If you received the work
electronically, the person or entity providing it to you may choose to
give you a second opportunity to receive the work electronically in
lieu of a refund. If the second copy is also defective, you may
demand a refund in writing without further opportunities to fix the
problem.

1.F.4. Except for the limited right of replacement or refund set forth
in paragraph 1.F.3, this work is provided to you ‘AS-IS’, WITH NO
OTHER WARRANTIES OF ANY KIND, EXPRESS OR IMPLIED,
INCLUDING BUT NOT LIMITED TO WARRANTIES OF
MERCHANTABILITY OR FITNESS FOR ANY PURPOSE.

1.F.5. Some states do not allow disclaimers of certain implied

warranties or the exclusion or limitation of certain types of damages.
If any disclaimer or limitation set forth in this agreement violates the
law of the state applicable to this agreement, the agreement shall be
interpreted to make the maximum disclaimer or limitation permitted
by the applicable state law. The invalidity or unenforceability of any
provision of this agreement shall not void the remaining provisions.

1.F.6. INDEMNITY - You agree to indemnify and hold the Foundation,

the trademark owner, any agent or employee of the Foundation,
anyone providing copies of Project Gutenberg™ electronic works in
accordance with this agreement, and any volunteers associated with
the production, promotion and distribution of Project Gutenberg™
electronic works, harmless from all liability, costs and expenses,
including legal fees, that arise directly or indirectly from any of the
following which you do or cause to occur: (a) distribution of this or
any Project Gutenberg™ work, (b) alteration, modification, or
additions or deletions to any Project Gutenberg™ work, and (c) any
Defect you cause.

Section 2. Information about the Mission

of Project Gutenberg™
Project Gutenberg™ is synonymous with the free distribution of
electronic works in formats readable by the widest variety of
computers including obsolete, old, middle-aged and new computers.
It exists because of the efforts of hundreds of volunteers and
donations from people in all walks of life.

Volunteers and financial support to provide volunteers with the

assistance they need are critical to reaching Project Gutenberg™’s
goals and ensuring that the Project Gutenberg™ collection will
remain freely available for generations to come. In 2001, the Project
Gutenberg Literary Archive Foundation was created to provide a
secure and permanent future for Project Gutenberg™ and future
generations. To learn more about the Project Gutenberg Literary
Archive Foundation and how your efforts and donations can help,
see Sections 3 and 4 and the Foundation information page at
www.gutenberg.org.

Section 3. Information about the Project

Gutenberg Literary Archive Foundation
The Project Gutenberg Literary Archive Foundation is a non-profit
501(c)(3) educational corporation organized under the laws of the
state of Mississippi and granted tax exempt status by the Internal
Revenue Service. The Foundation’s EIN or federal tax identification
number is 64-6221541. Contributions to the Project Gutenberg
Literary Archive Foundation are tax deductible to the full extent
permitted by U.S. federal laws and your state’s laws.

The Foundation’s business office is located at 809 North 1500 West,

Salt Lake City, UT 84116, (801) 596-1887. Email contact links and up
to date contact information can be found at the Foundation’s website
and official page at www.gutenberg.org/contact

Section 4. Information about Donations to

the Project Gutenberg Literary Archive
Foundation
Project Gutenberg™ depends upon and cannot survive without
widespread public support and donations to carry out its mission of
increasing the number of public domain and licensed works that can
be freely distributed in machine-readable form accessible by the
widest array of equipment including outdated equipment. Many
small donations ($1 to $5,000) are particularly important to
maintaining tax exempt status with the IRS.

The Foundation is committed to complying with the laws regulating

charities and charitable donations in all 50 states of the United
States. Compliance requirements are not uniform and it takes a
considerable effort, much paperwork and many fees to meet and
keep up with these requirements. We do not solicit donations in
locations where we have not received written confirmation of
compliance. To SEND DONATIONS or determine the status of
compliance for any particular state visit www.gutenberg.org/donate.

While we cannot and do not solicit contributions from states where

we have not met the solicitation requirements, we know of no
prohibition against accepting unsolicited donations from donors in
such states who approach us with offers to donate.

International donations are gratefully accepted, but we cannot make

any statements concerning tax treatment of donations received from
outside the United States. U.S. laws alone swamp our small staff.

Please check the Project Gutenberg web pages for current donation
methods and addresses. Donations are accepted in a number of
other ways including checks, online payments and credit card
donations. To donate, please visit: www.gutenberg.org/donate.

Section 5. General Information About

Project Gutenberg™ electronic works
Professor Michael S. Hart was the originator of the Project
Gutenberg™ concept of a library of electronic works that could be
freely shared with anyone. For forty years, he produced and
distributed Project Gutenberg™ eBooks with only a loose network of
volunteer support.
Project Gutenberg™ eBooks are often created from several printed
editions, all of which are confirmed as not protected by copyright in
the U.S. unless a copyright notice is included. Thus, we do not
necessarily keep eBooks in compliance with any particular paper
edition.

Most people start at our website which has the main PG search
facility: www.gutenberg.org.

This website includes information about Project Gutenberg™,

including how to make donations to the Project Gutenberg Literary
Archive Foundation, how to help produce our new eBooks, and how
to subscribe to our email newsletter to hear about new eBooks.
Welcome to our website – the ideal destination for book lovers and
knowledge seekers. With a mission to inspire endlessly, we offer a
vast collection of books, ranging from classic literary works to
specialized publications, self-development books, and children's
literature. Each book is a new journey of discovery, expanding
knowledge and enriching the soul of the reade

Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.

Let us accompany you on the journey of exploring knowledge and

personal growth!

testbankdeal.com