CE5610 Finite Element Analysis
January – May 2014
Tutorial – 4: Numerical Integration and Ordinary Differential
equations
Class work problems: 1a,2a,3
1. Calculate the following integrals using (a) Simpson rule (b) Composite Simpson
rule with 4 subintervals (c) Gaussian Quadrature with two evaluation points and
comment on the results
3
1
a. dx
0
0.8 x 0.5 x 2
2
cos
2
b. (0.5 x) sin 4 (0.5 x)dx
0
x
c. y exp x * (t n t m ) * exp t dt , with m = 2, n = 1, = 2 for x = 20
0
2. Evaluate the following integral using (a) Simpson rule (b) two-point Gauss
Quadrature for both the variables.
1 1
r
4 4
a. s drds
1 1
1 1
r s drds
2 4
b.
1 1
3. Find the constants, (ai ) in the polynomial,
u(x) = 1+a1x(x-3)+a2(x(x-3))2+ a3(x(x-3))3
that best approximates the solution to the following second order ordinary
differential equation:
d 2u
2
3x 4 , over 0 x 3 , given u(0) 1, u(3) 1 , using:
dx
a. Collocation method for equally spaced points
b. Least squares method
c. Sub-domain method using equally spaced sub-domain
and compare the results with the exact solution and draw conclusions.