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 root of a nonlinear equation f(x)=0f(x)=0, and state its order of convergence. Discuss two situations in which the method fails. [7]

(b) Using the Bisection method, find a root of the equation x34x9=0x^3 - 4x - 9 = 0 correct to three decimal places. Show the tabulation of iterations. [5]

roots-of-nonlinear-equationsbisection-methodnewton-raphson
2long12 marks

(a) Derive Newton's divided difference interpolation formula and show that the divided differences are symmetric with respect to their arguments. [6]

(b) The following table gives the values of a function f(x)f(x):

x56911
f(x)12131416

Using Lagrange's interpolation formula, estimate the value of f(x)f(x) at x=10x = 10. [6]

interpolationnewton-divided-differencelagrange-interpolation
3long12 marks

(a) Explain the Gauss-Seidel iterative method for solving a system of linear equations. State the condition of convergence (diagonal dominance) for the method. [5]

(b) Solve the following system of equations using the Gauss-Seidel iteration method, performing three iterations and starting with the initial approximation (0,0,0)(0,0,0):

20x+y2z=17,3x+20yz=18,2x3y+20z=2520x + y - 2z = 17,\quad 3x + 20y - z = -18,\quad 2x - 3y + 20z = 25

[7]

solution-of-linear-systemsgauss-eliminationiterative-methods
4long12 marks

(a) Derive the standard five-point finite-difference formula for the Laplace equation 2ux2+2uy2=0\dfrac{\partial^2 u}{\partial x^2} + \dfrac{\partial^2 u}{\partial y^2} = 0 over a square mesh. [6]

(b) Classify the partial differential equation Auxx+Buxy+Cuyy+f(x,y,u,ux,uy)=0A u_{xx} + B u_{xy} + C u_{yy} + f(x,y,u,u_x,u_y) = 0 as elliptic, parabolic or hyperbolic, and give one physical example of each type. [6]

numerical-solution-of-pdeslaplace-equationfinite-difference
B

Section B: Short Answer Questions

Attempt all / any as specified.

7 questions
5short8 marks

Find a real root of the equation cosx=3x1\cos x = 3x - 1 correct to four significant figures using the Secant method. Take the initial approximations x0=0.5x_0 = 0.5 and x1=0.6x_1 = 0.6 and tabulate the successive iterations.

roots-of-nonlinear-equationssecant-method
6short8 marks

Evaluate the integral 06dx1+x2\displaystyle\int_{0}^{6} \frac{dx}{1+x^2} by dividing the interval into six equal subintervals using (a) the Trapezoidal rule and (b) Simpson's 1/3 rule. Compare your results with the exact value.

numerical-integrationtrapezoidal-rulesimpsons-rule
7short8 marks

From the following table of values, find the first and second derivatives of f(x)f(x) at x=1.1x = 1.1 using Newton's forward difference formula:

x1.01.21.41.61.82.0
f(x)0.0000.1280.5441.2962.4324.000
numerical-differentiationfinite-difference
8short8 marks

Fit a straight line y=a+bxy = a + bx to the following data using the method of least squares and estimate the value of yy when x=2.5x = 2.5:

x12345
y1.21.83.34.56.3
curve-fittingleast-squares
9short8 marks

Using the fourth-order Runge-Kutta method, find an approximate value of yy at x=0.2x = 0.2 for the initial value problem dydx=x+y2\dfrac{dy}{dx} = x + y^2, y(0)=1y(0) = 1, taking a step size h=0.1h = 0.1.

numerical-solution-of-odesrunge-kutta-method
10short8 marks

Using the LU (Doolittle) decomposition method, factorize the coefficient matrix and solve the system:

2x+3y+z=9,x+2y+3z=6,3x+y+2z=8.2x + 3y + z = 9,\quad x + 2y + 3z = 6,\quad 3x + y + 2z = 8.
solution-of-linear-systemslu-decompositionmatrix-inverse
11short8 marks

Using the explicit (Bender-Schmidt) finite-difference scheme, solve the one-dimensional heat conduction equation ut=2ux2\dfrac{\partial u}{\partial t} = \dfrac{\partial^2 u}{\partial x^2} subject to u(0,t)=u(1,t)=0u(0,t)=u(1,t)=0, u(x,0)=sinπxu(x,0)=\sin\pi x for 0x10\le x\le 1. Take h=0.25h=0.25 and choose an appropriate time step satisfying the stability criterion; compute uu for the first two time levels.

numerical-solution-of-odesheat-equationnumerical-solution-of-pdes