Practice Questions of Numerical Integration

Q1). Evaluate the following integral

4

 (1  x  4 x  2 x 5 ) dx
3

2

(a) Analytically (b) single application of the trapezoidal rule (c) composite trapezoidal rule
(d) Simpson one third and three-eighth rules (e) Compute Percentage Error.

Q2). Suppose that the upward force of air resistance on a falling object is proportional to the square
of the velocity. For this case, the velocity can be computed as
gm  g cd 
v (t )  tanh  t
cd  m 

Where cd = a second-order drag coefficient, g = 9.81 ms-2, m = 68.1 kg, and cd = 0.25 kgm-1
(a) Apply analytical integration to determine how far the object falls in 10 sec.
(b) Make the same evaluation, but compute the integral by using the Simpson rule when no.
of sub-intervals equal to 6.

Q3). An 11-m beam is subjected to a load, and the shear force follows the equation
V ( x)  5  0.25x 2
Where V is the shear force and x is length in distance along the beam. We know that
V  dM dx , and M is the bending moment. Integration gives the relationship
x
M  M 0   V dx
0

If Mo is zero and x = 11, calculate M using (a) analytical integration, (b) trapezoidal rule and
Simpson rules (one-third, and three-eighth) when no. of intervals equal to 7.

Q4). The work produced by a constant temperature, pressure volume thermodynamic process can
be computed as
W   p dV
Where W is work, p is pressure, and V is volume. Using a combination of the trapezoidal
rule, Simpson rules (one-third, and three-eighth), use the following data to compute the work
in kJ (kJ = kN - m)

Pressure (kPa) 336 294.4 266.4 260.8 260.5 249.6 193.6 165.6
Volume (m3) 0.5 2 3 4 6 8 10 11

Q5). The total mass of a variable density rod is given by

m    ( x) Ac ( x) dx
Where m = mass,  (x) =density, Ac(x) = cross-sectional area, x = distance along the rod, and
L = the total length of the rod. The following data have been measured for a 10 m rod length.
Determine the mass in kilograms to the best possible accuracy.

1
 x, m 0 2 3 4 6 8 10
 , g/cm 4.00
3
3.95 3.89 3.80 3.60 3.41 3.30
Ac, cm2 100 103 106 110 120 133 150

Q6). A projectile is fired from the ground level, the acceleration of the projectile is described the
1
integral  f ( x) dx where f ( x)  2 x  1. Find the velocity of the projectile by using the
0
Simpson’s rule when interval is divided into seven sub-intervals and also find the total
distance covered by the projectile after 25 seconds and also compute the percentage error?

Q7). The cross-sectional area of a channel can be computed as

B


Ac  H ( y ) dy
0

Where B is the total channel width (m), H is the depth (m), and y is the distance from the
bank (m). In a similar fashion, the average flow Q (m3 s-1) can be computed as
B


Q  U ( y ) H ( y ) dy
0

Where U is the water velocity (m/s). Use these relationships and a numerical method to
determine Ac and Q for the following data

x, m 0 2 4 5 6 9
H, m 0.5 1.3 1.25 1.7 1 0.25
U, m/s 0.03 0.06 0.05 0.12 0.11 0.02

Q8). If a capacitor initially holds no charge, the voltage across it as a function of time can be
computed as
t
1
V (t ) 
C0 
i (t ) dt

If C = 10 - 5 farad, use the following current data to develop a plot of voltage versus time

t, sec 0 0.2 0.4 0.6 0.8 1.0 1.2

i, 10-3 A 0.2 0.3683 0.3819 0.2282 0.0486 0.0082 0.1441

Q9). Velocity data for air are collected at different radii from the centerline of a circular 16-cm-
diameter pipe as tabulated below

r, cm 0 1.60 3.20 4.80 6.40 7.47 7.87 7.95 8.0

v, ms-1 10 9.69 9.30 8.77 7.95 6.79 5.57 4.89 0

Use numerical integration to determine the mass flow rate, which can be computed as
R

 ρ v2π r dr
0
where ρ  density  1.2 kg m3 .

2
Q10). A river is 80 meter wide. The depth d in feet at a distance x ft. from one bank is given from
the following table.

x 0 10 20 30 40 50 60 70 80
y 0 4 7 9 12 16 14 7 4

Find approximately the area of cross-section. The area of cross-section can be computed as
x

A  f ( x) dx
0

Q11). A rocket is launched from the ground level, its acceleration is registered during the first 45
45

seconds is described by the integral  f ( x) dx where

0
f ( x)  x e  x , by applying the

Simpson compute the velocity of the rocket at 45 seconds when no. of interval is divided in
to seven sub-intervals..

Q12). A solid of revolution is formed by rotating about x axis, the area between the x axis, the
lines x = 0 and x = 1, and the curve passing through the points with the following coordinates

x 0 2.5 5.0 7.5 10.0 12.5 15.0

y 5 5.5 6.0 6.75 6.25 5.5 4.0

Estimate the volume of the solid is generated. The volume of the solid of revolutions is
computed as
15

V   y 2 dx 
0