UNIVERSITY OF HONG KONG
DEPARTMENT OF MATHEMATICS
MATH 3601 Numerical Analysis (2020–2021)
Tutorial 10
1. (Richardson extrapolation) By using the repeated Richardson extrapolation, find f 0 (1) from the
following values:
x 0.6 0.8 0.9 1.0 1.1 1.2 1.4
f (x) 0.707178 0.859892 0.925863 0.984007 1.033743 1.074575 1.127986
f (x+h) f (x h)
Apply the approximate formula f 0 (x) = 2h
with h = 0.4, h = 0.2, h = 0.1.
2. (Numerical method for ODE) The f 0 (x) can be approximated by the following finite di↵erence scheme
f 0 (x)
af (x) + bf (x h) + cf (x 2h)
gh (x) =
2h
(a) Determine the values of a, b, and c so that the finite di↵erent scheme gh (x) approximates f 0 (x)?
(b) What is the local truncation error of the above finite di↵erence scheme?
(c) Calculate f 0 (2.0) as accurately as possible by the Richardson extrapolation with the given data
(some data may be useless)
x 1.2 1.4 1.6 1.7 1.8 1.9 2.0
f (x) 0.550630 0.604826 0.665766 0.699730 0.736559 0.776685 0.820576
3. (Numerical method for ODE) Consider the IVP
x0 = f (t, x) := tx2 , x(0) = 2.
(a) Let h = 0.1. Use the second order Taylor series method to compute x(0.1).
(b) The tableau of a second order RK method is given below
1/3 0
2/3 3/4 3/4
Verify that the approximated value of x(0.1) agrees with the one in (a) if the same step size is
used.
dy
4. (Error analysis) Consider the initial value problem, dt
= f (t, y), y(t0 ) = y0 , t0 t T .
(a) The second-order two-stage Runge-Kutta method is of the form
8
< K1 = f (tk , yk ),
K2 = f (tk + ↵h, yk + hK1 )
:
yk+1 = yk + h(aK1 + bK2 ),
where a + b = 1, ↵b = 1/2, b = 1/2, and tk = t0 + kh. Do one step of the second-order Runge-
Kutta method with b = 1 to approximate y(1.2), where y(t) is the solution of dy dt
= t y3,
y(1) = 1.
(b) Assume the solution y(t) is exact at time t = tk , the Taylor expansion of y(t) at t = tk is
y(tk+1 ) = y(tk ) + hy 0 (tk ) + 12 h2 y 00 (tk ) + 16 h3 y 000 (⇠k ), where ⇠k lies between t = tk and t = tk+1 .
What is the local truncation error of the approximation method given in part (a)?
5. (Boundary value problem) Use a second order method for the solution y = y(x) of the boundary
value problem
y 00 = xy + 1, x 2 [0, 1],
with the boundary condition y 0 (0) + y(0) = 1 and y(1) = 1. The mesh size h = 0.25.
1