Browse papers
A

Section A: Long Answer Questions

Attempt all / any as specified.

4 questions
1long12 marks

(a) Derive the iterative formula for the Newton-Raphson method for finding a real root of the equation f(x)=0f(x)=0, and state its order of convergence. Mention two situations in which the method fails. [6]

(b) Find a real root of the equation x32x5=0x^3 - 2x - 5 = 0 correct to four decimal places using the Newton-Raphson method, taking x0=2x_0 = 2. Compare the number of iterations required with the Bisection method to achieve the same accuracy in the interval [2,3][2, 3]. [6]

roots-of-equationsnewton-raphsonbisection
2long12 marks

(a) Construct Newton's divided difference table for the following data and hence estimate the value of f(3)f(3).

xx01245
f(x)f(x)2312147326

[7]

(b) Explain why Newton's divided difference formula is preferred over Lagrange's interpolation formula when a new data point is added to an existing set. Illustrate your answer briefly. [5]

interpolationnewton-divided-differencelagrange
3long12 marks

(a) Derive the composite Simpson's 1/31/3 rule for numerical integration and write down the expression for its error term. [6]

(b) Evaluate 06dx1+x2\displaystyle\int_{0}^{6}\frac{dx}{1+x^2} by taking 6 sub-intervals using (i) the Trapezoidal rule and (ii) Simpson's 1/31/3 rule. Compare both results with the exact value and comment on the accuracy. [6]

numerical-integrationtrapezoidal-rulesimpsons-rule
4long12 marks

(a) Describe the fourth-order Runge-Kutta method for solving a first-order ordinary differential equation dydx=f(x,y)\dfrac{dy}{dx}=f(x,y) with the initial condition y(x0)=y0y(x_0)=y_0, clearly writing all four slope expressions. [6]

(b) Using the fourth-order Runge-Kutta method, compute y(0.2)y(0.2) for the initial value problem dydx=x+y,  y(0)=1\dfrac{dy}{dx}=x+y,\; y(0)=1, taking a step size h=0.1h=0.1. [6]

numerical-solution-odesrunge-kuttaeuler-method
B

Section B: Short Answer Questions

Attempt all / any as specified.

8 questions
5short7 marks

Solve the following system of linear equations using Gauss elimination with partial pivoting, and explain why partial pivoting is necessary:

2x+y+z=10,3x+2y+3z=18,x+4y+9z=16.2x + y + z = 10,\quad 3x + 2y + 3z = 18,\quad x + 4y + 9z = 16.
linear-systemsgauss-eliminationpartial-pivoting
6short7 marks

State the condition for convergence of the Gauss-Seidel iterative method. Apply the method to solve the following system, performing three iterations starting from (0,0,0)(0,0,0):

20x+y2z=17,3x+20yz=18,2x3y+20z=25.20x + y - 2z = 17,\quad 3x + 20y - z = -18,\quad 2x - 3y + 20z = 25.
linear-systemsgauss-seideliterative-methods
7short7 marks

Fit a straight line of the form y=a+bxy = a + bx to the following data using the method of least squares, and estimate yy when x=6x = 6.

xx12345
yy1427405568
curve-fittingleast-squareslinear-regression
8short6 marks

The following table gives the value of a function f(x)f(x). Using Newton's forward difference formula, find the first and second derivatives f(x)f'(x) and f(x)f''(x) at x=1.0x = 1.0.

xx1.01.11.21.31.4
f(x)f(x)7.9898.4038.7819.1299.451
numerical-differentiationforward-differencebackward-difference
9short6 marks

(a) Distinguish between round-off error and truncation error with a suitable example of each. [3]

(b) The number π=3.14159265\pi = 3.14159265 is rounded to 3.14163.1416. Compute the absolute error, the relative error, and the percentage error. [3]

error-analysisround-off-errortruncation-error
10short6 marks

(a) Write down the iterative formula of the Secant method and explain how it differs from the Newton-Raphson method. [3]

(b) Perform two iterations of the Secant method to find a root of f(x)=cosxxex=0f(x) = \cos x - x e^{x} = 0 taking initial approximations x0=0x_0 = 0 and x1=1x_1 = 1. [3]

roots-of-equationssecant-methodfixed-point-iteration
11short7 marks

Using Newton's forward difference interpolation formula, estimate the population in the year 1925 from the following census data, and state the assumption on which the formula is based.

Year19111921193119411951
Population (in thousands)1215202739
interpolationnewton-forward-differencefinite-differences
12short6 marks

Using the Modified Euler (Heun's) method, compute y(0.1)y(0.1) and y(0.2)y(0.2) for the initial value problem dydx=x2+y,  y(0)=1\dfrac{dy}{dx} = x^2 + y,\; y(0) = 1 with step size h=0.1h = 0.1, iterating the corrector until two successive values agree to three decimal places.

numerical-solution-odesmodified-eulerpredictor-corrector