Mathematical Principles for
Scientific Computing and Visualization

by Gerald Farin & Dianne Hansford

published by A K Peters, Ltd., May 2008

ISBN: 978-1-56881-321-9

Errata
If you find errors in the text or figures, please dianne@farinhansford.com!

 

 

Chapter 2 Computational Basics Posted
Page 7,
Section 2.2

There is an inconsistency in the description of the IEEE standard for 64-bit number representation. The range of exponents for 64-bit numbers is -308 to 308. The range given (-38 to 38) is that of 32-bit numbers.

Found by Mark McLaughlin, Walt Disney Animation Studios

8/26/08
Page 8,
Section 2.2,
Example 2

Strictly speaking the truncated polynomial would be of degree 2n-1 although it would only have n terms.

Found by Mark McLaughlin, Walt Disney Animation Studios

8/26/08
Page 8,
Section 2.3,
Example 3

The equation shown, 1/10 = .0001100110011..., might be unclear.
This expression for 1/10 is
1/2^4 + 1/2^5 + 1/2^8 + 1/2^9 + 1/2^{12} + 1/2^{13} ...
or 1/16 + 1/32 + 1/256 + 1/512 + 1/4096 + 1/8192 ...

Found by Mark McLaughlin, Walt Disney Animation Studios

8/26/08

 

Chapter 3 Coordinate Systems Posted
Page 14

Figure 3.1 clarification: each grid element represents 1/5th in x and y.

Found by Mark McLaughlin, Walt Disney Animation Studios

9/15/08
Page 20 Should be "20 latitudinal bands" rather than "longitudinal pieces" 9/1/16

 

Chapter 4 Background: Numerical Linear Algebra Posted
Page 31

Last paragraph of section 4.2.
Note that m <= n.

Found by Mark McLaughlin, Walt Disney Animation Studios

9/15/08

 

Chapter 5 Solving Linear Systems Posted
Page 47

The clause "this the norm of choice" should be "this is the norm of choice" .

Found by Mark McLaughlin, Walt Disney Animation Studios

9/15/08
Page 51

First paragraph: "This is, will the sequence..." should be "That is, will the sequence..." .

Found by Mark McLaughlin, Walt Disney Animation Studios

9/15/08

 

Chapter 6 Eigen-Problems Posted
Page 59

For the given eigenvalues, the eigenvectors are reversed. They should be:
u_1 = [1,1]^T and u_2 = [1, -1]^T.
Modify the column vectors of U to reflect this change, thus making the second row +0.707 and -0.707.

Found by Mark McLaughlin, Walt Disney Animation Studios
and Eric Schafer, ILM

9/15/08
Page 60

The vector sequence equation should be: v^(i) = Av^(i-1)

Found by Mark McLaughlin, Walt Disney Animation Studios

9/15/08
Page 61

In Step 1 we pick v^(0). To be consistent with the development of method, this should be v^(1).

Found by Mark McLaughlin, Walt Disney Animation Studios

9/15/08
Page 64

The (2,2) element of the matrix U is incorrect. The columns of the matrix should be [.707,.707]^T and [-.707, .707]^T

On page 59 the matrix U holds a rotation and reflection. This is possible because of the degree of freedom in the solution to the homogeneous system (A-lambda I)u_2 = 0.

Found by Matthias Schweinoch

11/12/14

 

Chapter 9 Computing Dynamic Processes Posted
Page 121 Figure 9.9 was created with 2r(t) -fr rather than the expression in (9.12) 8/1/08