DEPARTMENT OF MATHEMATICS

SCHOOL OF ADVANCED SCIENCES

Fall Semester – 2017 ~ 2018
Continuous Assessment Test – I, Aug. 2017
Course Code : MAT 3005 Applied Numerical Methods Slot: C1 + TC1+TCC1
Max. Marks : 50 Duration: 90 Minutes
---------------------------------------------------------------------------------------------------------------------
Answer ALL Questions
(5 x 10 marks = 50 marks)

1 N
1. Show that the iteration formula for finding square root of N is xn 1   xn   and
2 xn 

hence find the value of 15 using Newton’s method correct to three decimal places.

2. Solve the following system of nonlinear equations by Newton’s method correct to three
decimal places: x 2 + y = 11 and y 2 + x = 7, starting with (x0,y0) = (2.5, 1.5) up to three
iterations, correct to three decimal places.

3. Solve the following system of equations using Gauss Seidel method, correct to three
decimal places:
8 x  3 y  2 z  20; 4 x  11y  z  33;6 x  3 y  12 z  35.

4. Solve the following tri - diagonal system of equations by Thomas algorithm.

2x1 – 2x2 = 1, x1 – 2x2 + 3x3 = 2, 2x2 – 2x3 + 4x4 = 1, x3 – x4 = 3.

5. From the following table values of x and f(x), determine using Newton’s interpolation

(i) f(0.23) and (ii) f(0.29).

x 0.20 0.22 0.24 0.26 0.28 0.30

f(x) 1.6596 1.6698 1.6804 1.6912 1.7024 1.7139

***********
 Department of Mathematics
School of Advanced Sciences
Winter Semester 2018-19
Continuous Assessment Test - II, March - 2019
Course Code: MAT3005 Slot: A1+TA1+TAA1
Course Title: Applied Numerical Methods xxxxx
Max Marks: 50 Duration: 90 Minutes

Answer ALL the questions

1. A plane is being tracked by radar, and data are taken every second in polar coordi- (15 M)
nates θ and r.

r
t, s 200 202 204 206 208 210

pe
θ, rad 0.75 0.72 0.70 0.68 0.67 0.66
r, m 5120 5370 5560 5800 6030 6240

Pa
If the velocity and acceleration expressions in polar coordinates are given by:
   
2
~v = ṙ~er + rθ̇~eθ and ~a = r̈ − rθ̇ ~er + rθ̈ + 2ṙθ̇ ~eθ ,
n
find the vector expressions for velocity ~v , and acceleration ~a at t = 202 s.
io
2. The outflow concentration from a reactor is measured at a number of times over a (10 M)
st

24−hr period:
ue

t, hr 0 1 5 10 12 14 16 18 20 24
c, mg/L 1 1.5 2.3 2.1 4 5 5.5 5 3 1.2
Q

The flow rate for the outflow in m3 /s can be computed with the following equation:
 
2π
el

Q(t) = 20 + 10 sin (t − 10) ,

24
od

where the arguments of the sine are expressed in degrees. Use the best numerical
integration technique to determine the weighted concentration
M

Z t
wc = Q(t)c(t) dt
0

leaving the reactor over the 24−hr period.

3. The following differential equation can be used to compute the position and velocity (13 M)
of a soccer ball that is kicked straight up in the air with an initial velocity of 40 m/s:
d2 y
 
Cd dy  dy 
+g+ =0
dt2 m dt  dt 
where y = upward distance (m), t = time (s), v = dy dt
upward velocity (m/s), g =
gravitational constant (= 9.81 m/s2 ), Cd = drag coefficient (kg/m) and m = mass
(kg). Find the displacement and velocity of the ball at t = 0.5 seconds by assuming
m = 0.4 kg, Cd = 0.52.

Page 1 of 2
4. The temperature distribution in a tapered conical cooling fin is described by the (12 M)
following differential equation, which has been nondimensionalised

d2 u
  
2 du
+ − pu = 0
dx2 x dx

where u = temperature (0 ≤ u ≤ 1), x = axial distance (0 ≤ x ≤ 1), and p is a nondi-

mensional parameter that describes the heat transfer and geometry. Considering the
boundary conditions u(x=0) = 0, u(x=1) = 1, solve the equation for the temperature
distribution for p = 10 and by dividing the interval [0, 1] into 4 sub-intervals.

r
pe
Pa
n
io
st
ue
Q
el
od
M

Page 2 of 2
 SCHOOL OF ADVANCED SCIENCES

DEPARTMENT OF MATHEMATICS

Final Assessment Test

Course Code : MAT 3005 Slot :

Course Name : Applied Numerical Methods

Max. Marks : 100 Answer Any FIVE Questions Duration: 3hrs.

--------------------------------------------------------------------------------------------------------------------

1 a. Find the root of the equation 𝑐𝑜𝑠𝑥 − 𝑥𝑒 𝑥 = 0 by

Secant method correct to 4 decimal places. (8)

b. Solve the system of equations 𝑙𝑜𝑔𝑒 (𝑥 2 + 𝑦) − 1 + 𝑦 = 0 𝑎𝑛𝑑 𝑥𝑦 + √𝑥 = 0 has an

starting with (𝑥0 , 𝑦0 ) = (2.4, −0.6) (12)

2. a For the steady-state case, a differential equation based on heat conservation can be
written for such a system as 𝑇 ′′ +ℎ′ (𝑇𝛼 − 𝑇) = 0. where T = temperature (𝐶 𝑜 , ) 𝑥 =
distance along the rod (m), h’= a heat transfer coefficient between the rod and the
surrounding air (𝑚−2), and 𝑇𝛼 = the air temperature (𝐶 𝑜 ). If h’= 0.01, 𝑇𝛼 = 20, 𝑇0 =40, and 𝑇5
= 200, the governing differential equation yields the system of equations −𝑇0 +
2.04 𝑥𝑇1 − 𝑇2 = 0.8; −𝑇1 + 2.04𝑇2 − 𝑇3 = 0.8; −𝑇2 + 2.04𝑇3 − 𝑇4 = 0.8 and −𝑇3 +
2.04𝑇4 − 𝑇5 = 0.8. Hence find the heat at interior nodes of the rod of length 10cm the
nodes are taken at equal distance ∆𝑥 = 2 by Thomas algorithm. (12)

−30 10 20
b. The stress matrix describing the state of stress at a point is ( 10 40 −50)
20 −50 −10
.Determine the largest principal stress by the power method. (8)

3. a. From the following data, estimate the number of students who obtained marks
between 40 and 45 (8)

Marks 30 − 40 40 − 50 50 − 60 60 − 70 70 − 80
No. of students 31 42 51 35 31
b. Fit the following 4 points by the cubic splines

𝑖 0 1 2 3
𝑥𝑖 1 2 3 4
𝑦𝑖 1 5 11 8
′′ ′′
Use the end conditions 𝑦0 = 𝑦3 = 0 . Hence compute 𝑦(1.5) 𝑎𝑛𝑑 𝑦 ′ (2). (12)

4. a. A slider in a machine moves along a fixed straight rod. Its distance 𝑥 𝑐𝑚 along the rod
is given below for various values of the time 𝑡 𝑠𝑒𝑐𝑜𝑛𝑑𝑠. Find the velocity of the slider and
its acceleration when 𝑡 = 0.3. (6)

t 0 0.1 0.2 0.3 0.4 0.5 0.6

x 30.13 31.62 32.87 33.64 33.95 33.81 33.24
1
1 𝑥
b. Using Romberg integration with step size ℎ = 16 calculate ∫02 𝑑𝑥. (8)
𝑠𝑖𝑛𝑥

1 log(1+𝑥)
c. Compute by Gaussian quadrature 𝐼 = ∫0 𝑑𝑥. (6)
√𝑥(1−𝑥)

5 a. Apply the Adams- predictor corrector method with ℎ = 0.2 and starting values from
the R.K.method to initial value problem 𝑦 ′ = 𝑦 − 𝑡 2 + 1, 0 ≤ 𝑡 ≤ 2 𝑤𝑖𝑡ℎ 𝑦(0) = 0.5.
Given that 𝑦(0.2) = 0.82929; 𝑦(0.4) = 1.2141 𝑎𝑛𝑑 𝑦(0.6) = 1.64892 and hence find
𝑦(1.0). (12)

b. Using the finite difference method, find 𝑦(0.25), 𝑦(0.5) 𝑎𝑛𝑑 𝑦(0.75) satisfying the
differential equation 𝑦 ′′ + 𝑦 = 𝑥 subject to the boundary conditions 𝑦(0) = 0, 𝑦(1) = 2.

(8)

𝜕𝑢 𝜕2 𝑢
6. Classify the nature of the partial differential equation = 𝜕𝑥 2 𝑖𝑛 0 < 𝑥 < 5; 𝑡 ≥ 0 given
𝜕𝑡
that 𝑢(𝑥, 0) = 20, 𝑢(0, 𝑡) = 0; 𝑢(5, 𝑡) = 100. Compute 𝑢 for the time step with ℎ = 1, by
deriving Crank-Nicholson difference Scheme. (20)

𝜋
7. a. Find the extermal of the function𝐼[𝑦(𝑥)] = ∫0 (𝑦 ′2 − 𝑦 2 )𝑑𝑥 with boundary condition
𝜋
𝑦(0) = 0, 𝑦(𝜋) = 1 and subject to the constraint ∫0 𝑦 𝑑𝑥 = 1. (10)

b.Using Galerkins method solve 𝑦 ′′ + 𝑦 + 𝑥 = 0 ,0 < 𝑥 < 1 with the boundary condition
𝑦(0) = 0 = 𝑦(1). (10)