Computational Physics

Errata

Errata to versions 2.0+1.0

• Section 2.1, p 77, $\vec r(t)= x(t) \,\hat x+y(y)\, \hat y \changeto \vec r(t)= x(t) \,\hat x+y(t)\, \hat y$
• Equations (14.77) - (14.80): Lower limit of integration is $-\infty$: $\int_{\infty}^{+\infty} \changeto \int_{-\infty}^{+\infty}$
• Problem 4.4: Case (c) is to be worked out only numerically and not to be compared with an analytic solution.

Errata to versions 1.0

(Fortran edition, found in the printed versions, vol. I and vol. II. Corrected in version 1.1 or higher.)

• 12 lines after Equation (3.15): $0\le |f'(x_2)|=|2-r|<1 \changeto 0\le |f'(x_2^*)|=|2-r|<1$
• 6 lines after Equation (3.18): $f^{(2)}{}'(x_{4}^2)= 1 \changeto f^{(2)}{}'(x_{4}^*)= 1$
• Page 150, line 4: $r_c^{(2)} < r< r_c^{(3)}\approx 3.5441 \changeto r_c^{(3)} < r < r_c^{(4)} \approx 3.5441$
• Before Equation (3.36): $\lambda=\lim\limits_{n\to\infty}\lim\limits_{\epsilon\to 0}|\Delta x_n|/\epsilon \changeto \lambda=\lim\limits_{n\to\infty}\lim\limits_{\epsilon\to 0}\frac{1}{n}\ln |\Delta x_n|/\epsilon$
• Equation (3.36): $\tilde x_2 = f(\tilde x_1) = f(x_0 + \epsilon f'(x_0)) \changeto \tilde x_2 = f(\tilde x_1) = f(x_1 + \epsilon f'(x_0))$
• Equation (3.36): $\tilde x_3 = f(\tilde x_3) \changeto \tilde x_3 = f(\tilde x_2)$
• Chapter 4, footnote 1:We have $N$ discrete time points $t_i\equiv t_1, \ldots, t_{N-2}, t_{N-1}\equiv t_f$ $\changeto$ We have $N$ discrete time points $t_i\equiv t_1, \ldots, t_{N-1}, t_{N}\equiv t_f$
• Chapter 5, section 5.1, program rk2.f90, subroutine RKSTEP: Correct bug in code:
  k14=f1(t+h ,x1+h *k13,x2+h *k23,x3+h *k33,x4+h2*k43) $\changeto$ k14=f1(t+h ,x1+h *k13,x2+h *k23,x3+h *k33,x4+h *k43)
  and similarly correct k24, k34, k44. The same corrections must be applied to the program scatter.f90.
• Equation (11.2): $z \changeto x_{i-1}$
• Chapter 9, p. 375, Equation (9.14), $p_{n m} \changeto -p_{n m}$
• Chapter 9, p. 395, code of program doublewell.f90, 5 lines from bottom: (add minus sign)
  if(j.ge.1) iP(i,j) = -isqrt2*sqrt(DBLE(m+1))
• Chapter 9, p. 396, code of program doublewell.f90, first line: (remove minus sign)
  if(j.le.DIM) iP(i,j) = isqrt2*sqrt(DBLE(m))
• Chapter 13, p. 518, file init.f90: Replace the lines:
  if( beta < 0.0D0 ) beta = obeta $\changeto$if(beta < 0.0D0)call locerr('beta has not been set.')
read(f_in,'(A4,I5,A4,I5,A6,G28.17,A6,25I16)')buf,OL,buf,OL,buf,beta,buf,seeds $\changeto$read(f_in,'(A4,I5,A4,I5,A6,G28.17,A6,25I16)')buf,OL,buf,OL,buf,obeta,buf,seeds
• Chapter 13, p. 520: Remove the sentence: This way, if the user has not set the temperature ... temperature will be used.