MA6459 Numerical Methods
MA6459 Numerical Methods
com
10. Can we apply iteration method to find the root of the equation 2 x cos x 5 in 0 , ?
2
11. Find the positive root of x 2 + 5x -3 = 0 using fixed point iteration method starting with 0.6 as
first approximation
12. Find the root of xex 3 = 0 in 1< x < 1.1 by Iteration method.
13. State the condition for Convergence of Iteration method.
1 3
Power method.
2 7
15.Find the first iteration values of x, y, z by Gauss seidel method
28x 4 y z 32 ; x 3 y 10z 24 ; 2 x 17 y 4 z 35
14. Find inverse of A =
16. State the condition for the convergence of Gauss Seidel iteration method for solving a system of
linear equation
17. By Gauss Elimination method solve (i)x + y =2 and 2x + 3y = 5 (ii) 2x y =1 , x 3y +2 = 0
(iii) x 2y = 0 , 2x + y = 5
18. Give two direct method to solve system of linear equations.
19. Compare Gauss Elimination, Gauss Jordan method.
20. Compare Gauss seidel method , Gauss Jordan method.
PART B
1. (i) Find the positive real root of 3x cosx 1 = 0 using Newton Raphson method.
1 6 1
(ii) Find the dominant eigen value and vector of A = 1 2 0 using Power method.
0 0 3
2. (i) Find the positive real root of 2x log10 x - 6 = 0 using Newton Raphson method.
1 2 3
(ii) Find the dominant eigen value and vector of A = 0 4 2 using Power method.
0 0 7
3. (i) Find the positive real root of x3 5x +3= 0 using Iteration method.
1
3
1
4. (i) Find the positive real root of 2x3 -3x 6 = 0 using Iteration method.
2
3 1
(ii) Find the inverse of the matrix 2 3 1 using Gauss Jordan method.
1 2
1
1 1 2
7. (i) Find the inverse of the matrix 1 2 3 using Gauss Jordan method.
2 3 1
2
3 1
8. (i) Find the inverse of the matrix 2 3 1 using Gauss Jordan method.
1 2
1
(ii) Apply Gauss seidel method to solve system of equations 30x 2y +3z = 75 ;
2x + 2y + 18z = 30 ; x + 17 y 2z =48
25 1 2
9. (i) Find the dominant eigen value and vector of A = 1 3 0 using Power method.
2 0 4
(ii) Apply Gauss seidel method to solve system of equations 28x + 4y z =32 ;
x + 3y + 10z =24 ; 2x + 17y + 4z =35
1 2 1
(ii) Apply Gauss seidel method to solve system of equations x + y + 54z = 110 ;
27x + 6y z = 85 ; 6x + 15y 2z =72
********************************
UNIT - II INTERPOLATION AND APPROXIMATIOIN
PART A
1. State Lagranges interpolation formula for unequal intervals.
2. Using Lagranges find a polynomial through (0,0) (1 ,1 ) and ( 2 ,2 ).
3. Give inverse Lagranges interpolation formula.
4.Write the Lagranges formula for y, if three sets of values
are given?
7. Form the divided difference table for the following data ( 0,1) , ( 1,4) ,( 3,40) and ( 4,85)
8. Find the interpolating polynomial for the given data: x :
9.
f(x): 1
8
3
f(x):48
10
11
13
100
294
900
1210
2028
y:
12
x:
y:
24
Temp :
140
150
160
170
180
Pressure: 3.685
4.854
6.302
8.076
10.225
(ii) Fit a polynomial using Inverse Lagranges method x: 1
2
7
y: 1
5
5
8. (i) Find y(0.5) and y(1) given that M0 = M2 = 0 using Cubic Spline
x: 0
1
2
y: -5
-4
3
(ii) Find the polynomial of degree 3 from the following using Newtons formula
X: 0
1
2
3
4
5
6
7
Y: 1
2
4
7
11
16
22
29
9. (i) Find y (43) from the following
X: 40
50
60
70
80
90
Y: 184
204
226
250
276
304
(ii) Using Cubic Spline to the following data, Find y(1.5)
x:
1
2
3
4
y:
1
2
5
11
10.(i) Find the polynomial and hence find f(1) using Lagranges method
x: -1
0
2
3
y: -8
3
1
12
2
3
4
5
(ii) Find the spline interpolation x: 1
y: 1
0
1
0
1
8
4
**************************************
UNIT III NUMERICAL DIFFERENTIATION AND INTEGRATION
d2y
dy
1. State Newtons backward difference formula to find
and 2
dx x xn
dx x xn
2. What do you mean by numerical differentiation.
1
1
1 x
1 x
dx
dx
16. Write down the abscissa and weights of the three point Gaussian integration.
1
rd
rule to find I =
1 x 2 dx taking h = 0.1
19. Evaluate
th
rule x : 1
f(x) : 1
27
64
(ii) By dividing the range into 10 equal parts , evaluate sin xdx using Simpsons 1/3 rule.
0
2.(i) Obtain first and second derivative of y at x = 0.96 from the data
x:
0.96
0.98
1
1.02
1.04
y:
0.7825 0.7739 0.7651
0.7563
0.7473
6
1 x
(ii) Evaluate
3.(i) The table given below reveals the velocity of the body during the time t specified .Find
its acceleration at t =1.1
t : 1.0
1.1
1.2
1.3
1.4
v: 43.1
47.7 52.1 56.4 60.8
2.4 4.4
(ii) Evaluate
rd
dx dy
x y
5. (i) From the following table find the value of x for which f (x) is maximum
x:
60
75
90
105
120
f(x):
28.2 38.2 43.2 40.9 37.7
1.5
r 2
dr
0.2
6.(i) Find the value of f (8) from the table given below
x: 6
7
9
12
f(x) : 1.556
1.690 1.908
2.158
1.4
2.4
(ii) Evaluate
7. (i) Evaluate
0
dx dy
1 x y
dy
d2y
and
at x = 3.5 from the data given below by Newtons backward
dx
dx 2
difference formula x : 1.5
2
2.5
3
3.5
4
y : 3 . 375
7
13.625 24 38.875
59
/
//
8.(i) Obtain f (0.02 ) and f (0.05)
x : 0.01
0.02
0.03
0.04
0.05
0.06
y : 0.1023
0.1047
0.1071
0.1096
0.1122
0.1148
(ii) Calculate
dt
1 t
1.2
8.781
1.3
1.4
1.5
9.129
9.451
9.75 find
dy
d2y
and
at x = 1.1
dx
dx 2
dx
1 x
10.(i) A river is 80 meter wide the depth d in meters at a distance x meters from one bank
is given below. Calculate the area of the cross section of the river using Simpson rule
x:0
10
20
30
40
50
60
70
80
d:0
4
7
9
12
15
14
8
3
3
(ii) Find the first second derivative of the function f ( x) = x 9x -14 at x = 3.0 using the
values given below x : 3
3.2
3.4
3.6
3.8
4
f(x) : -14 -10.03 -5.296
-0.256
-6.672
14
dy
x y when
dx
with
given that
dy
y x 2 given
dx
, for
.Take
dy
1 2 xy y(0) = 0 with h = 0.2.
dx
dy
f ( x, y), y( x0 ) y 0
dx
12. Use modified Eulers method to find y (0.4) given
dy
x y y(0) = 1, h = 0.1
dx
14. Write down the Runge Kutta method of order 4 for solving initial value problems in
ordinary differential equation
15. Use R-K method of second order to find y (0.4) given
13 .Apply fourth order Runge Kutta method to find y(0.1) given
dy 4t
t. y given y(0) = 3 calculate y at x = 0.1 and 0.2.
dx y
dy
2.(i) Using Taylor series method find y at x = 0.1 given
2 y 3e x y(0) = 0.
dx
for
(ii) Solve 2y x y = 0 given y(0) = 2 , y(0.5) = 2.636, y(1) = 3.595, y(1.5) = 4.968 to get y(2)
by Adams method.
3.(i) Using Taylor series method find correct to 4 decimal places the value of y(0.1) given
dy
x 2 y 2 y(0) =1.
dx
dy
(ii) Obtain y(0.6) given
x y ,y(0) = 1 using h = 0.2 by Adams method if
dx
y(-0.2) = .8373, y(0.2) = 1.2427, and y(0.4) = 1.5834.
4.(i) Using Eulers method , solve numerically the equation y = x + y, y(0) = 1
for x = 0.0 (0.2) (1.0). Check your answer with the exact solution.
(ii) Solve and find y(2) by Milnes method
dy 1
( x y ) , given y(0) = 2 , y(0.5) = 2.636,
dx 2
dy
x 2 y 2 y(0) =1and h = 0.1
dx
dy
x 2 (1 y ) y(1) = 1 , y(1.1) = 1.233, y(1.2) = 1.548 , y(1.3) = 1.979, evaluate y(1.4)
dx
dy
xy y 2 y(0) = 1 , using Taylor series method
dx
dy 1
1 x 2 y 2 and y(0) = 1, y(0.1) = 1.06, y(0.2) = 1.12, y( 0.3) = 1.21, evaluate y(0.4)
dx 2
dy
y x 2 , y(0) = 1. Find y(0.1) and y(0.2) by R-K method of fourth order , find y(0.3)
dx
dy
x 2 y 1 y(0) =1.
dx
(ii) Given
dy
xy y 2 y(0) =1,y(0.1) = 1.1169, y(0.2) = 1.2773, y(0.3) = 1.5049, evaluate y(0.4)
dx
by Milnes method
dy y 2 x 2
10(i) Using Runge-kutta method of 4 solve
given that y(0) = 1 at x = 0.2 and 0.4
dx y 2 x 2
th
(ii)Given
dy
xy y 2 y(0) =1,y(0.1) = 1.1169, y(0.2) = 1.2773, y(0.3) = 1.5049,
dx
1.
PART-A
Define difference quotient of a function y (x)
d2y
and state the order of truncation error
dx 2
d2y
3y 2
3. Write down the finite difference scheme for the differential equation
dx 2
d2y
y5
4. Obtain the finite difference scheme for the differential equation 2
dx 2
2. State the finite difference approximation for
5. Express y (x) + a(x) y(x) + b(x) y(x) = f(x) in terms of difference quotients
6. Write standard five point formula and diagonal five point formula used in solving Laplace
equation Uxx + U yy = 0 at the point ( ix , jy)
7. Write down Laplaces equation and its finite difference analogue and the standard five
point formula
8. Write down the diagonal five point formula in Laplace equation
9. Write down the two dimensional Laplaces equation and Poissons equation
10. Write down Poissons equation and its finite difference analogue
11. What is the order of error in solving Laplace and Poissons equation by using finite
difference method?
12. Write down the finite difference scheme for solving the Poissons equation
13. State one dimensional heat equation and its boundary conditions
14. Name at least two numerical methods that are used to solve one dimensional diffusion
equation
15. State the implicit finite difference scheme for one dimensional heat equation
16. What is Bender Schmidt recurrence equation? For what purpose it is used?
17. Write down one dimensional wave equation and its boundary conditions
18. State the explicit formula for the one dimensional wave equation with 1 2 a2 = 0 where =
k / h and a2 = T /m
2
2u
2 u
19. State explicit finite difference scheme for one dimensional wave equation
t 2
x 2
at 0oand choosing x =y =1. Start with zero vector and do 4 Gauss seidel iterations.
2.(i) Solve ytt = 4yxx subject to the condition y(0,t) = 0; y(2,t) = 0; y(x,0) = x(2-x) ; ut(x,0) = 0, Do 4
steps. Find the values upto 2 decimal accuracy.
(ii) Solve u(x, 0) = 3x 0 x 4, u(x,4) = x2 0 x 4, u(0,y) = 0 0 x 4,
u(4,y) = 12 + y 0 x 4,
3.(i) Solve y = x + y with the conditions y(0) = y(1) = 0 by finite difference method, taking h = 0.25
2 x , when 0 x 2.5
2 u 2 u u
(ii) 25 2 2 ,
( x,0) 0 u(0, t ) 0 , u(5, t ) 0 and u( x,0)
t
x
t
10 2 x, when 2.5 x 5
by the method derived above taking h =1 and for one period of vibration , ( i.e. up to t =2)
4.(i) Solve the Laplace equation Uxx+ Uyy = 0 for the square mesh with the boundary values , using
Leibmanns method correct to integers
0 100
200
100 0
0 100 200
100 0
5. (i)Solve the elliptic equation Uxx + Uyy = 0 for the following square mesh with boundary values
as shown , using Liebmans iteration procedure .
0
11.1 17 19.7
18.6
6.(i) Solve Uxx + Uyy = 0 over the square mesh of side 4 units, satisfying the following conditions
u(x, 0) = 3x 0 x 4, u(x,4) = x2 0 x 4, u(0,y) = 0 0 x 4, u(4,y) = 12 + y 0 x 4,
7.(i) Solve the 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 is 1.
u
2u
4 2 and the boundary
8(i) Find the value of u( x, t ) satisfying the parabolic equation
t
x
2
x
conditions u(0,t) = 0 , u(8,t) = 0 and u(x,0) = 4x at x = i, i = 0,1,2.7and t = 1/8j j = 0,1,2,3
2
y (2.5 ) and y (2.75)
(ii) Solve : Uxx + Uyy =0 given
A
1
2u 2u
2 0 x 1 subject to the initial
t 2
x
conditions u ( 0 ,t ) = 0 , u ( 1 , t ) =0 , ut(x,0) = 0 , u( x ,0 ) = sin x , 0 x 1,
(ii) Solve :
u 2 u
2 , given u ( 0 ,t ) = 0 , u ( 4 , t ) =0 , u( x ,0 ) = x ( 4 x ) , assuming h = k =1
t
x
u 2 u
2 given u ( 0 ,t ) = 0 , u ( 5 , t ) =0 ,
10(i) Using Bender Schimidt formula solve :
t
x
u( x ,0 ) = x2 ( 25 x2 ) , assuming x =1 .Find the value of u upto t =5
(ii) Solve by Crank Nicholsons method the equation 16 Ut =Uxx 0 <x <1 and t > 0subject to
u(x , 0) = 0 , u ( 0,t )=0 and u ( 1 ,t ) =100 t .Compute one time step , taking x =