Engineering Mathematics 2014

SUBJECT NAME : Numerical Methods

SUBJECT CODE : MA2264
MATERIAL NAME : University Questions
MATERIAL CODE : JM08AM1004
REGULATION : R2008
UPDATED ON : November 2014 (Upto N/D 2014 Q.P)

(Scan the above Q.R code for the direct download of this material)

Name of the Student: Branch:

Unit – I (Solution of Equations and Eigenvalue Problems)

 Fixed point iteration

1) Solve y  3 e x  3 x  0 by the method of fixed point iteration. (M/J 2012)

2) Find a real root of the equation x 3  x 2  1  0 by iteration method. (M/J 2012)

3) Find a positive root of the equation cos x  3 x  1  0 by using iteration method.

(M/J 2013)

 Newton’s method (or) Newton Raphson method

1) Solve for a positive root of the equation x 4  x  10  0 using Newton – Raphson
method. (A/M 2010)

2) Find by Newton-Raphson method a positive root of the equation 3 x  cos x  1  0 .

(N/D 2014)

3) Solve x sin x  cos x  0 using Newton-Raphson method. (M/J 2012)

4) Find an iterative formula to find the reciprocal of a given number N and hence find the
1
value of . (N/D 2011)
19

5) Find theNewton’s iterative formula to calculate the reciprocal N and hence find the
1
value of . (N/D 2012)
23

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 1

Engineering Mathematics 2014

6) Find an iterative formula to find N , where N is a positive number and hence find
5. (N/D 2012)

 Solution of linear system by Gaussian elimination method

1) Apply Gauss elimination method to find the solution of the following system:
2 x  3 y  z  5, 4 x  4 y  3z  3, 2 x  3 y  2z  2 . (N/D 2012)

 Solution of linear system by Gaussian-Jordon method

1) Apply Gauss-Jordan method to find the solution of the following system:

10 x  y  z  12; 2 x  10 y  z  13; x  y  5z  7 . (N/D 2011)

2) Using Gauss-Jordan method to solve 2 x  y  3z  8;  x  2 y  z  4;

3 x  y  4z  0 . (N/D 2014)

3) Solve the system of equation by Gauss-Jordan method: 5 x1  x2  9 ,

 x1  5 x2  x3  4 ,  x2  5 x3  6 . (M/J 2014)

 Solution of linear system by Gaussian-Seidel method

1) Using Gauss-Seidel method, solve the following system of linear equations
4 x  2 y  z  14 , x  5 y  z  10 , x  y  8z  20 . (M/J 2014)

2) Apply Gauss-Seidal method to solve the equations

20 x  y  2z  17; 3 x  20 y  z  18; 2 x  3 y  20z  25 . (M/J 2012),(N/D 2014)

3) Solve, by Gauss-Seidel method, the following

system: 28 x  4 y  z  32; x  3 y  10z  24; 2 x  17 y  4z  35 . (N/D 2011)

4) Use Gauss – Seidal iterative method to obtain the solution of the equations:
9 x  y  2z  9; x  10 y  2z  15; 2 x  2 y  13z  17 . (A/M 2010)

5) Solve the following system of equations using Gauss-Seidel method:

10 x  2 y  z  9; x  10 y  z  22;  2 x  3 y  10z  22 . (N/D 2012)

6) Solve, by Gauss-Seidel method, the following

system: 8 x  3 y  2z  20; 4 x  11 y  z  33; 6 x  3 y  12z  35 . (N/D 2012)

7) Solve, by Gauss-Seidel method, the equations 27 x  6 y  z  85; 6 x  15 y  2z  72;

x  y  54z  110 . (M/J 2013)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 2

Engineering Mathematics 2014

 Inverse of a matrix by Gauss Jordon method

2 2 6 

1) Using Gauss-Jordan method, find the inverse of 2 6 6 .
 (M/J 2014)
 
 4 8 8 

1 1 3

2) Using Gauss Jordan method, find the inverse of the matrix 1 3 3  .(M/J 2012)

 2 4 4 

 0 1 2
 
3) Find the inverse of the matrix by Gauss – Jordan method:  1 2 3  . (A/M 2010)
 3 1 1
 

1 2 6 
 
4) Using Gauss Jordan method, find the inverse of the matrix 2 5 15 .(N/D 2012)
 
 6 15 46 

4 1 2
 
5) Find, by Gauss-Jordan method, the inverse of the matrix A  2 3 1 .(M/J 2013)
 
 1 2 2 

 Eigen value of matrix by power and Jacobi method

 25 1 2 
 
1) Find the numerically largest eigenvalue of A  1 3 0 and its corresponding
 
 2 0 4 
eigenvector by power method, taking the initial eigenvector as  1 0 0  .
T

(M/J 2014),(N/D 2014)

2) Determine the largest eigenvalue and the corresponding eigenvector of the matrix

 2 1 0 
  1 2  1 . (M/J 2012)
 
 0 1 2 
1 6 1
 
3) Find the largest eigenvalue of  1 2 0  , by using Power method.
 0 0 3
 
(A/M 2010),(N/D 2011),(M/J 2012),(N/D 2012)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 3

Engineering Mathematics 2014

 1 2 2 
 
4) Find all the eigenvalues and eigenvectors of the matrix  2 3 2  using Jacobi
 
 2 2 1 
method. (N/D 2012)

5) Using Jacobi method find the all eigenvalues and their corresponding eigenvectors of the
 2 3
matrix A   . (M/J 2013)
 3 2

Unit – II (Interpolation and Approximation)

 Lagrange Polynomials and Divided differences

1) Using Lagrange’s interpolation, calculate the profit in the year 2000 from the following
data: (M/J 2012)

Year: 1997 1999 2001 2002

Profit Lakhs Rs. 43 65 159 248

2) Using Lagrange’s interpolation formula find the value of y when x  10 , if the values
of x and y are given below: (M/J 2012)

x : 5 6 9 11
y : 12 13 14 16

3) Find the polynomial f ( x ) by using Lagrange’s formula and hence find f (3) for the
following values of x and y: (N/D 2014)
x: 0 1 2 5
y: 2 3 12 147

4) Apply Lagrange’s formula, to find y(27) to the data given below. (M/J 2013)
x : 14 17 31 35
y : 68.8 64 44 39.1
5) Use Lagrange’s formula to find a polynomial which takes the values
f (0)  12, f (1)  0, f (3)  6 and f (4)  12 . Hence find f (2) . (A/M 2010)

6) Using Lagrange’s interpolation formula, find y(2) from the following data:
y(0)  0; y(1)  1; y(3)  81; y(4)  256; y(5)  625 . (M/J 2014)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 4

Engineering Mathematics 2014

7) Use Lagrange’s method to find log10 656 , given that

log10 654  2.8156, log10 658  2.8182, log10 659  2.8189 and
log10 661  2.8202 . (N/D 2012)

8) Determine f ( x ) as a polynomial in x for the following data, using Newton’s divided

difference formulae. Also find f (2) . (N/D 2011)

x: –4 –1 0 2 5
f ( x ) : 1245 33 5 9 1335

9) Find f (3) by Newton’s divided difference formula for the following data: (M/J 2014)

x: –4 –1 0 2 5
f ( x ) : 1245 33 5 9 1335

10) Find the function f ( x ) from the following table using Newton’s divided difference
formula: (A/M 2010)
x: 0 1 2 4 5 7
f ( x ) : 0 0 –12 0 600 7308

11) By using Newton’s divided difference formula find f (8) , given (N/D 2014)

x: 4 5 7 10 11 13
f ( x ) : 48 100 294 900 1210 2028

12) Using Newton’s divided differences formula determine f (3) from the data:

x: 0 1 2 4 5 (M/J 2012)
f ( x ) : 1 14 15 5 6

13) Use Newton’s divided difference formula to find f ( x ) form the following data.

x: 1 2 7 8 (M/J 2013)
y: 1 5 5 4

14) Using Newton’s divided difference formula, find f ( x ) from the following data and
hence find f (4) . (N/D 2012)

x 0 1 2 5

f ( x) 2 3 12 147

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 5

Engineering Mathematics 2014

 Interpolating with a cubic spline

1) Fit the cubic splines for the following data: (M/J 2014)

x: 1 2 3 4 5
y: 1 0 1 0 1

2) The following values of x and y are given:

x: 1 2 3 4
y : 1 2 5 11
Find the cubic splines and evaluate y(1.5) and y(3) . (M/J 2012)

3) From the following table:

x: 1 2 3
y : - 8 -1 18
Compute y(1.5) and y(1) using cubic sphere. (M/J 2012),(M/J 2013)

4) Obtain the cubic spline for the following data to find y(0.5) . (N/D 2012),(N/D 2014)

x : -1 0 1 2
y : -1 1 3 35

5) If f (0)  1 , f (1)  2, f (2)  33 and f (3)  244 . find a cubic spline approximation,
assuming M(0) = M(3) = 0. Also, find f (2.5) . (A/M 2010)

 Newton’s forward and backward difference formulas

1) Find the cubic polynomial which takes the following values:
(M/J 2012),(M/J 2013),(N/D 2014)

x: 0 1 2 3
f ( x ) : 1 2 1 10

2) Find the value of y when x using Newton’s interpolation formula from the following
table: (N/D 2012)
x 4 6 8 10
y 1 3 8 16

3) Find the value of y at x  21 and x  28 from the data given below (N/D2012)

x: 20 23 26 29
y : 0.3420 0.3907 0.4384 0.4848

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 6

Engineering Mathematics 2014

4) The population of a town is as follows:

x Year: 1941 1951 1961 1971 1981 1991
y Population in
thousands: 20 24 29 36 46 51
Estimate the population increase during the period 1946 to 1976. (N/D 2011)

5) Given the following table, find the number of students whose weight is between 60 and
70 lbs: (A/M 2010),(M/J 2012)

Weight (in lbs) x: 0 – 40 40 – 60 60 – 80 80 – 100 100 – 120

No. of students: 250 120 100 70 50

Unit – III (Numerical Differentiation and Integration)

 Differentiation using interpolation formulae

1) A slider in a machine moves along a fixed straight rod. Its distance x cm along the rod is
given below for various values of the time ' t ' seconds. Find the velocity of the slider
when t  1.1 second. (M/J 2012)

t : 1.0 1.1 1.2 1.3 1.4 1.5 1.6

x : 7.989 8.403 8.781 9.129 9.451 9.750 10.031

2) For the given data, find the first two derivatives at x  1.1 . (M/J 2014)

x : 1.0 1.1 1.2 1.3 1.4 1.5 1.6

y : 7.989 8.403 8.781 9.129 9.451 9.750 10.031

3) Find f ( x ) at x  1.5 and x  4.0 from the following data using Newton’s formulae
for differentiation. (N/D2012)

x: 1.5 2.0 2.5 3.0 3.5 4.0

y  f ( x ) : 3.375 7.0 13.625 24.0 38.875 59.0
4) Find the first three derivatives of f ( x ) at x  1.5 by Newton’s forward interpolation
formula to the data given below. (M/J 2013)

x: 1.5 2.0 2.5 3.0 3.5 4.0

y  f ( x ) : 3.375 7.0 13.625 24.0 38.875 59.0

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 7

Engineering Mathematics 2014

5) Given the following data, find y(6) and the maximum value of y (if it exists)

x: 0 2 3 4 7 9 (A/M 2010)
y : 4 26 58 112 466 922

6) Find the first two derivatives of x1/ 3 at x  50 and x  56 , for the given table:
(N/D 2011)

x: 50 51 52 53 54 55 56
y  x1/ 3 3.6840 3.7084 3.7325 3.7563 3.7798 3.8030 3.8259

7) Find the first and second derivative of the function tabulated below at x  0.6 .
(N/D 2012)

x: 0.4 0.5 0.6 0.7 0.8

y : 1.5836 1.7974 2.0442 2.3275 2.6511

dy d2 y
8) Find and at x  51 , from the following data: (M/J 2012)
dx dx 2

x : 50 60 70 80 90
y : 19.96 36.65 58.81 77.21 94.61

 Numerical integration by trapezoidal, Simpson’s 1/3, 3/8

rules & Romberg’s method
1
dx
1) Evaluate  1 x
0
and correct to 3 decimal places using Romberg’s method and hence

find the value of log e 2 . (N/D 2014)

1
dx
2) Use Romberg’s method to compute  1 x
0
2
correct to 4 decimal places. Also evaluate

the same integral using tree-point Gaussian quadrature formula. Comment on the
obtained values by comparing with the exact value of the integral which is equal to  / 4 .
(M/J 2012)

1
dx
3) Evaluate  1 x
0
2
by using Romberg’s method correct to 4 decimal places. Hence

deduce an approximate value for  . (M/J 2012)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 8

Engineering Mathematics 2014

1
dx
4) Using Romberg’s integration to evaluate  1 x
0
2
. (A/M 2010)

1
2
x
5) Evaluate  sin x dx correct to three decimal places using Romberg’s method. (M/J 2014)
0

1
dx
6) Using Trapezoidal rule, evaluate  1 x
1
2
by taking eight equal intervals. (M/J 2013)

 /2
7) Compute  sin x dx using Simpson’s 3/8 rule.
0
(N/D 2012)

2
8) Evaluate I   sin x dx by dividing the range into ten equal parts, using
0

(i) Trapezoidal rule (N/D2012)

(ii) Simpson’s one-third rule, Verify your answer with actual integration.

6
1
9) Evaluate I   1  x dx by using (i) Direct integration (ii) Trapezoidal rule (iii) Simpson’s
0

one-third rule (iv) Simpson’s three-eighth rule. (N/D 2011)

1.3
10) Taking h  0.05 , evaluate 
1
x dx using Trapezoidal rule and Simpson’s three-eighth

rule. (M/J 2014)

11) The velocity  of a particle at a distance S from a point on its path is given by the table
below:
S (meter) 0 10 20 30 40 50 60
 (m / sec) 47 58 64 65 61 52 38
Estimate the time taken to travel 60 meters by Simpson’s 1/3rd rule and Simpson’s 3/8th
rule. (A/M 2010),(N/D 2014)

 Two and Three point GaussianQuadrature formulae

5
dx
1) Use Gaussian three-point formula and evaluate 
1
x
. (M/J 2012)

2
dx
2) Evaluate  1 x
1
3
using 3 point Gaussian formula. (N/D 2014)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 9

Engineering Mathematics 2014

x2  2 x  1
2
3) Evaluate  1   x  1
0
2
dx by Gaussian three point formula. (M/J 2013)

 Double integrals using Trapezoidal and Simpsons’s rules

2 2
dxdy
1) Using Trapezoidal rule, evaluate  x
1 1
2
 y2
numerically with h  0.2 along x -

direction and k  0.25 along y -direction. (M/J 2012)

1 1
1
2) Evaluate   1  x  y dxdy by trapezoidal rule.
0 0
(N/D 2014)

1.2 1.4
dxdy
3) Evaluate 
1 1
x y
by trapezoidal formula by taking h  k  0.1 . (A/M 2010)

2 1
1 1
4) Evaluate   4 xy dxdy using Simpson’s rule by taking h  4 and k  2 .(N/D 2012)
0 0

1/ 2 1/ 2
sin( xy ) 1
5) Evaluate 
0 0
1  xy
dxdy using Simpson’s rule with h  k  .
4
(M/J 2012),(M/J 2014)

1.2 2.4
1
6) Evaluate 
1 2
xy
dxdy using Simpon’s one-third rule. (M/J 2013)

Unit – IV (Initial Value Problems for Ordinary Differential Equations)

 Taylor series method

1) Use Tailor series method to find y(0.1) and y(0.2) given that
dy
 3e x  2 y, y(0)  0 , correct to 4 decimal accuracy. (A/M 2010)
dx

dy
2) Using Taylor series method solve  x 2  y , y(0)  1 at x  0.1, 0.2, 0.3 . Also
dx
compare the values with exact solution. (M/J 2012)

3) Using Taylor series method to find y(0.1) if y  x 2  y 2 , y(0)  1 . (M/J 20113)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 10

Engineering Mathematics 2014

4) Using Taylor’s series method, find y at x  1.1 by solving the equation

dy
 x 2  y 2 ; y(1)  2 . Carryout the computations upto fourth order derivative.
dx
(M/J 2014)
dy
5) Using Taylor’s series method, find y at x  0 if  x 2 y  1 , y(0)  1 . (N/D 2014)
dx

 Euler method for first order equation

y x
1) Solve y  , y(0)  1 at x  0.1 by taking h  0.02 by using Euler’s method.
y x
(M/J 2013)

2) Apply modified Euler’s method to find y(0.2) and y(0.4) given y  x 2  y 2 , y(0)  1
by taking h  0.2 . (N/D 2014)

dy
3) Using Modified Euler’s method, find y(4.1) and y(4.2) if 5 x  y 2  2  0;
dx
y(4)  1 . (N/D 2012)

 Fourth order Runge – Kutta method

dy y 2  x 2
1) Use Runge – Kutta method of fourth order to find y(0.2) , given  ,
dx y 2  x 2
y(0)  1 , taking h  0.2 . (A/M 2010)

2) Consider the second order initial value problem y  2 y  2 y  e 2 t sin t with
y(0)  0.4 and y(0)  0.6 using Fourth order RungeKutta algorithm, find y(0.2) .
(M/J 2012)

3) Solve y  xy  y given y(0)  3, y(0)  0 to find the value of y(0.1) , using RK-
method of order 4. (M/J 2012)

4) Given y  xy  y  0, y(0)  1, y(0)  0 . Find the value of y(0.1) by using Runge-
Kutta method of fourth order. (N/D 2011),(N/D 2014)

5) Given y  x  y   y 2  0 , y(0)  1 , y(0)  0 . Find the value of y(0.2) by using

2

Runge-Kutta method of fourth order, by taking h  0.2 . (N/D 2012)

6) Using Runge-Kutta method find y(0.2) if y  xy 2  y 2 , y(0)  1 , y(0)  0, h  0.2 .
(M/J 2013)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 11

Engineering Mathematics 2014

7) Solve for y(0.1) and z(0.1) from the simultaneous differential equations
dy dz
 2 y  z;  y  3 z; y(0)  0, z(0)  0.5 using Runge-Kutta method of the
dx dx
fourth order. (N/D 2012)

 Milne’s and Adam’s predictor and corrector methods

1) Using RungeKutta method of fourth order, find the value of y at x  0.2, 0.4, 0.6 given
dy
 x 3  y , y(0)  2 . Also find the value of y at x  0.8 using Milne’s predictor and
dx
corrector method. (M/J 2014)

2) Given that y  xy  y  0, y(0)  1, y(0)  0 obtain y for x  0.1,0.2 and 0.3 by
Taylor’s series method and find the solution for y(0.4) by Milne’s method.(M/J 2012)

  1  x 2  y 2 ; y(0)  1; y(0.1)  1.06; y(0.2)  1.12 and

dy 1
3) Given that
dx 2
y(0.3)  1.21 , evaluate y(0.4) and y(0.5) by Milne’s predictor corrector method.
(N/D 2011)
dy
4) Given that  1  y 2 ; y(0.6)  0.6841, y(0.4)  0.4228, y(0.2)  0.2027,
dx
y(0)  0 , find y( 0.2) using Milne’s method. (N/D 2012)

dy  1  x  y
2 2

5) Use Milne’s predictor – corrector formula to find y(0.4) , given  ,

dx 2
y(0)  1, y(0.1)  1.06, y(0.2)  1.12 and y(0.3)  1.21 . (A/M 2010)

1
6) Given y  , y(0)  2, y(0.2)  2.0933, y(0.4)  2.1755 , y(0.6)  2.2493
x y
find y(0.8) using Milne’s method. (M/J 2012)
7) Given 5 xy  y 2  2, y(4)  1, y(4.1)  1.0049, y(4.2)  1.0097, y(4.3)  1.0143 .
Compute y(4.4) using Milne’s method. (N/D 2014)
dy
8) Given  xy  y 2 , y(0)  1, y(0.1)  1.1169, y(0.2)  1.2774, y(0.3)  1.5041 .
dx
Use Adam’s method to estimate y(0.4) . (A/M 2010)

9) Using Adams method find y(1.4) given y  x 2 (1  y ), y(1)  1, y(1.1)  1.233,

y(1.2)  1.548 and y(1.3)  1.979 . (M/J 2012)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 12

Engineering Mathematics 2014

10) Given that y  y  x 2 ; y(0)  1; y(0.2)  1.1218; y(0.4)  1.4682 and

y(0.6)  1.7379 , evaluate y(0.8) by Adam’s predictor-corrector method. (N/D 2012)

11) Using Adam’s method to find y(2) if y  ( x  y ) / 2, y(0)  2, y(0.5)  2.636,

y(1)  3.595, y(1.5)  4.968 . (M/J 2013)
12) Using Adam’s Bashforth method, find y(4.4) given that 5 xy  y  2, 2

y(4)  1, y(4.1)  1.0049, y(4.2)  1.0097 and y(4.3)  1.0143 . (M/J 2014)

Unit – V (Boundary Value Problems in ODE and PDE)

 Finite Difference Solution of Second Order ODE

1) Solve the equation y  x  y with the boundary conditions y(0)  y(1)  0 using
finite differences by dividing the interval into four equal parts. (M/J 2012),(M/J 2014)

2) Solve y  y  0 with the boundary conditions y(0)  0 and y(1)  1 . (N/D 2012)

3) Solve, by finite difference method, the boundary value problem y( x )  y( x )  0 ,

where y(0)  0 and y(1)  1 , taking h  0.25 . (M/J 2012)

4) Using the finite difference method, compute y(0.5) , given

y  64 y  10  0, x  (0,1), y(0)  y(1)  0 , subdividing the interval into (i) 4 equal
parts (ii) 2 equal parts. (N/D 2011),(N/D 2012)

 One Dimensional Heat equation by explicitmethod

 2 u u
1) Solve  , subject to u(0, t )  u(1, t )  0 and u( x,0)  sin  x, 0  x  1 ,
x 2  t
using Bender-Schmidt method. (M/J 2012)

 2 u u
2) Using Bender-Schmidt’s method Solve  , given
x 2  t
u(0, t )  0, u(1, t )  0 u( x,0)  sin  x, 0  x  1 and h  0.2 . Find the value of
u upto t  0.1 . (M/J 2014)

3) Solve by Bender-Schmidt formula upto t  5 for the equation uxx  ut , subject to

 
u(0, t )  0, u(5, t )  0, and u( x,0)  x 2 25  x 2 , taking h  1 . (N/D 2012)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 13

Engineering Mathematics 2014

u 1  2 u
4) Solve  with the condition u(0, t )  0  u(4, t ) , u( x,0)  x(4  x ) taking
t 2 x 2
h  1 employing Bender-Schmidt recurrence equation. Continue the solution through
10 time steps. (M/J 2012)

5) Solve uxx  32ut , h  0.25 for t  0, 0  x  1, u(0, t )  0, u( x,0)  0, u(1, t )  t .

(M/J 2013)

 One Dimensional Heat equation by implicit method

kc 2
1) Obtain the Crank – Nicholson finite difference method by taking    1 . Hence,
h2
 2 u u
find u( x , t ) in the rod for two times steps for the heat equation  , given
x 2  t
u( x,0)  sin( x ), u(0, t )  0, u(1, t )  0 . Take h = 0.2. (A/M 2010)

 2 u u
2) Solve by Crank-Nicolson’s method  for 0  x  1, t  0 given that
x 2  t
u(0, t )  0, u(1,t)  0 and u( x,0)  100 x(1  x ) . Compute u for one time step with
1 1
h  and k  . (N/D 2014)
4 64

 One Dimensional Wave equation

 2u  2u
1) Solve the equation  , 0  x  1, t  0 satisfying the conditions
x 2 t 2
u 1
u( x , 0)  0, ( x, 0)  0, u(0, t )  0 and u(1, t )  sin  t . Compute u( x , t ) for 4
t 2
1
time-steps by taking h  . (N/D 2012)
4

 2u  2u
2) Solve the wave equation  , 0  x  1, t  0 ,
x 2 t 2
1, 0  x  0.5 u
u(0, t )  u(1, t )  0,t  0 u( x , 0)   and ( x , 0)  0, using
 1, 0.5  x  1 t
h  k  0.1 , find u for three time steps. (M/J 2014)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 14

Engineering Mathematics 2014

1  2u  2u
3) Find the pivotal values of the equation  with given conditions
4 t 2 x 2
u
u(0, t )  0, u(4, t )  0 , u( x,0)  x(4  x ) and ( x , 0)  0 by taking h  1 , for 4
t
time steps. (M/J 2012)
4) Solve 4utt  uxx , u(0, t )  0, u(4, t )  0 , u( x,0)  x(4  x ), ut ( x,0)  0, h  1
upto t  4 . (M/J 2013)

 Two Dimensional Laplace and Poisson equations

1) Solve the elliptic equation uxx  uyy  0 for the following square mesh with boundary
values as shown: (M/J 2012)

2) *Deduce the standard five point formula for  2 u  0 . Hence, solve it over the square
region given by the boundary conditions as in the figure below: (A/M 2010)

 2u  2u
3) By iteration method, solve the elliptic equation   0 over a square region of
x 2  y 2
side 4, satisfying the boundary conditions.

(i) u(0, y )  0, 0  y  4 (ii) u(4, y )  12  y, 0  y  4

(ii) u( x,0)  3 x, 0  x  4 (iv) u( x,4)  x 2 , 0  x  4

By dividing the square into 16 square meshes of side 1 and always correcting the
computed values to two places of decimals, obtain the values of u at 9 interior pivotal
points. (N/D 2014)

 
4) Solve the Poisson equation  2 u  10 x 2  y 2  10 over the square mesh with sides
x  0, y  0, x  3 and y  3 with u  0 on the boundary and mesh length 1 unit.
(N/D 2012),(M/J 2014)

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 15

Engineering Mathematics 2014

5) Solve  2 u  8 x 2 y 2 for square mesh given u  0 on the four boundaries dividing the
square into 16 sub-squares of length 1 unit. (N/D 2011)

6) Solve  2 u  8 x 2 y 2 over the square x  2, x  2, y  2, y  2 with u  0 on the

boundary and mesh length = 1. (M/J 2013)

----All the Best----

Prepared by C.Ganesan, M.Sc., M.Phil., (Ph:9841168917) Page 16