Numerical Methods Two Marks
Numerical Methods Two Marks
x rp - N
= xr - --------
Pxp-1
(p-1)xrp + N
= -----------------
Pxp-1
9. Say true or false: Newton’s method is useful when the graph of the
function when it crosses the x – axis is nearly vertical. (True).
10. Say true or false: Newton’s method is useful in cases where the
graph of the function when it crosses the x – axis is nearly
horizontal. ( False ) .
12. What is condition for applying the fixed point iteration (successive
approximation method) method to find the real root of the equation
x=f(x) (or)
If g(x) is continuous in [a, b], then under what conditions the
iterative method x=g(x) has a unique solution in [a, b]?
Let x=r be a root of x=g(x).Let I be an interval
combining the point x=r. If |g’(x)| < 1 for all x in I, the sequence of
approximation x1, x2, …. X n will converge to the root r, provided that the
initial approximation x 0 is chosen in I.
16. What is the condition for the convergence of the iteration method for
solving x=Φ Φ(x)?
Soln: Φ (x) < 1 in the range.
’
7. (a)Solve for positive root of the equation x-cos x=0, by Newton’s method.
(b).Solve by Gauss elimination method, the system of equations
3 x – y + 2 z = 13 , X + 2 y + z = 3, 2 x + 3 y + 3 z = 10.
8. (a) Solve by Gauss seidel method,
X + y + 54z = 110, 27 x + 6 y – z = 85, 6 x + 15 y + 3 z = 72.
2 2 3
(b).Using Gauss Jordan method,find inverse of the matrix 2 1 1
1 3 5
9. (a) Using Gauss Jacobin method, solve
X + y + 5 z = 16, 2 x + 3 y + z = 4, 4 x + y – z = 4.
(b).Using Gauss elimination method solves the following equations.
4 x + 4 y – 3 z = 4; 10 x + 8 y – 6 z = 5; 20 x – 4 y + 22 z = 7.
10.(a) Solve the equation x 6 – 5 x 2 + 136 = 0 and y 4 – 3 x 4 y + 80 = 0
assuming that ( 2 , 3 ) is an approximate solution. Obtain the solution
correct to 2 decimals.
(b). Solve X + 3 y + 3 z = 16, X + 4 y + 34 z = 18, X + 3 y + 4 z = 19,
by Gauss Jordan method.
M.KUMARASAMY COLLEGE OF ENGINEERING
THALAVAPALAYAM,KARUR.
UNIT-2 TWO MARKS
(x-x0) (x-x1)
+ -------------- y2
(x2-x1) (x2-x1)
( y – y 0 ) ( y – y 2 )…… ( y – y n )
+ x 1
( y 1 – y 0 ) ( y 1 – y 2 )…. ( y 1 – y n )
+……………………………………….
( y – y 0 ) ( y – y 2 )………( y – y n – 1 )
+ x n .
( y n – y 0) ( y n – y 2 )……( y n – y n – 1 )
2! 3!
Where ν =(x-xn)/n
2! 3!
Where u=(x-x0)/n
2! 3!
Where u=(x-x0)/n
2! 3!
Where ν =(x-xn)/n
2! 3!
Where ν =(x-xn)/n
Here x = x n ⇒ν = 0
Y( x ) = y n
Y ( x + h ) – y (x)
Y ′( x ) = yn′=lim
h→0 h
y′′(x)=y n′′
14. Newton’s forward interpolation formula used only for ----- intervals.
Equidistant or equal intervals.
15.Say true or false:
Newton’s interpolation formula are not suited to estimate the
value of a function near the middle of a table.
Soln: true
2! 3!
Where u=(x-x0)/n
Here X 0 = 0 , h = 1 – 0 = 1 , u = x
x x(x-1)
Y(x)=-2+ (4)+ (2)
1! 2!
=-2+4x+x(x-1)
=-2+4x+x2-x
=x2+3x-2
∴The roots of the equation f( x ) = 0 i.e x 2+ 3x - 2 = 0 are
- 3±√9+8 -3±√17
X = =
2 2
18. Obtain the interpolation polynomial for the given data by using
Newton’s forward difference formula.
X: 0 2 4 6
Y: -3 5 21 45
x y ∆y ∆2y ∆3y
0 -3
8
2 5 8
16 0
4 21 8
24
6 45
X–x0 x-0 x
u = = =
h 2 2
u (u-1)
y(x)=y0+n∆y0+ ∆2y0
2
(x/2) (x/2) ((x/2)-1)
Y(x)=-3+ (8)+ (8)
1! 2!
2
(x /4)- (x/2)
=-3+4x+ (8)
2!
2
=-3+4x+x -2x
=x2+2x-3
1. From the table below, Find the number of students who secured marks
between 40 and 45
Marks; 30 - 40 40 - 50 50 - 60 60 - 70 70 – 80
No.of students; 35 48 70 40 22
2. Define quadrature.
The process of evaluating of definite integral from a set of
tabulated values of a function of single variable is known as quadrature
(d 3 y / d 3 x ) x = x n = 1 / h 3 [∇ 3 y n + 3 / 2 ( ∇ 4 y n ) + …………. ]
1
6. Evaluate ∫ (1/x) dx by trapezoidal rule dividing the range in to 4 equal
parts. 0.5
Here h = [ 1 - ( 1 / 2 ) ] / 4 = 1 / 8 ,
Y = 1 / x . x 0 = 4 / 8 = 1 / 2 . X n =1,
X 1 = 5 / 8, x 2 = 6 / 8, x 3 = 7 / 8, x 4 = 8 / 8.
Y 1 = 1 / 4, y 2 = 1 / 5, y 3 = 8 / 6. y 4 = 8 / 7, y 5 = 8 / 8
A= sum of the first and last ordinates
= 8/4 +8/8 = 3.
B = sum of the remaining ordinates.
= 8/5 + 8/6 +8/7
= 856/210
1
∫ (1/x) dx = h/2(A+2B)
0 .5
=1/16( 3+( 856 *2/210)).
= 0.6971.
π
∫ sinx dx = h/2[ ( y0 + yn)+ 2( remaining terms)]
0
= ( π / 12 ){ 7.464}
=0.622π.
15. A curve passes through ( 2,8) (3,27) (4,64) and (5, 125) . Find the
area
of the curve between X axis and the lines x= 2 and x=5 by
trapezoidal
rule.
5
∫ y d x = 1/2[(8+125) +2(27+64)]
2 = 157.5 sq units.
12. Evaluate ∫ ∫dxdy /( 1+x+ y) , using Simpson’s method with step length
of
h=k= 0.5. 0 0
13. Find the value of sec ( 31 ). From the following table.
θ ( deg) : 31• 32• 33 • 34 •
tanθ :0.6008 0.6249 0.6494 0.6746
6
14. Using Simpson’s 1/3 rd rule ,Evaluate ∫ dx/(1+x2), by dividing the range
in
to 6 equal parts. 0
15. Find the maximum and minimum value of y .
X : -2 -1 0 1 2 3 4
Y: 2 -0.25 0 -0.25 2 15.75 56
16. The following data gives the velocity of a particle for 20 sec at an
interval
of 5 sec, Find the initial acceleration using the data.
Time( sec) : 0 5 10 15 20
Velocity(m/sec) : 0 3 14 69 228
1
17. Using Trapezoidal rule Evaluate ∫ dx( 1+x2) , using 8 intervals and also
by
-1
3 point Gaussian formula..
18. Obtain the value of f ′(0.04) using an approximate formula for the
given
data :
X: 0.01 0.02 0.03 0.04 0.05 0.06
F(x) : 0.1023 0.1047 0.1071 0.1096 0.1122 0.1148
2 2
19. Using Trapezoidal rule , Evaluate ∫ ∫ 1/(X +Y) dx dy. Taking 4 sub
intervals. 1 1
20. Find first and second derivative at x=1.2 for the data.
X: 1 1.1 1.2 1.3 1.4 1.5
1.6
Y : 7.989 8.403 8.781 9.12 9.451 9.750
10.031
5 2
21. Using three point formula Evaluate ∫(1/x)dx & ∫(x2+ 2x +1)dx/(
(1+(x+1)4)).
5 1 0
22. Evalute ∫dx/(4x+5) , by Simpson’s rule and hence find the value of log
e5 (n=10) 0
2 2
23. Evaluate ∫ ∫ f( x , y ) dx dy by Trapezoidal rule for the following data.
0 0
y x 0 0.5 1 1.5 2
0 2 3 4 5 5
1 3 4 6 9 11
2 4 6 8 11 14
NUMERICAL METHODS
UNIT-IV
INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL
EQUATIONS
8. State the third order R.K method algorithm to find the numerical
solution of the First order differential equation.
To solve the differential equation dy / dx= f ( x , y ) by the third
order R.K method , we use the following algorithm
K 1= h f ( x , y )
K 2= h f [ x + h / 2 , y + k 1/ 2 ]
K 3= h f [ x + h , y + 2 k 2– k 1]
and ∆y = 1 / 6 [ k 1 + 4 k 2 + k 3 ] .
9. Write down the Runge - Kutta formula for the fourth order to solve
dy / dx = f ( x , y ) with y ( x 0 ) = y 0.
Let h denotes the interval between equidistant values of x. If the
initial values are (x0, y0), the first increment in y is computed from the
formulas. K 1 = h f ( x 0 , y 0 )
K 2= h f [ x 0+ h / 2 , y 0+ k 1/ 2 ]
K 3= h f [ x 0+ h / 2 , y 0+ k 2/ 2 ]
K 4 = h f [ x 0+ h , y 0+ k 3]
∆ y =1 / 6 [ k 1 + 2 k 2 + 2 k 3 + k 4 ]
x 1 = x 0 + h ,y 1 = y 0 + ∆ y .
19. Given two examples for single-step and multi step method.
Single step methods:
(i) Taylor series method.
(ii) Euler method.
(iii) Runge-Kutta method.
Multi step Method:
(i) Milne’s method.
(ii) Adams-Bash forth method.
9. For the following mesh in solving ∇2u = 0, find one set of rough values of u at
interior mesh points
1 2
U1 U2
U3 U4
4 5
By symmetry, u2 = u3.
Assume u2 = 3
The rough values are u1 = (1/4)(1 +1 + 2u2) = 2
U2 = 3
U3 = (1/4)(5 +5+ 2u2) = 4
U4 = (1/4)(u1 + u4 + 2 + 4) = 3
10. Write the Laplace equations uxx + uyy = 0 in difference quotients.
uxx + uyy = 0
⇒ ui-1,j – 2ui,j + ui+1,j + ui,j-1 – 2ui,j + ui,j +
h2 k2
11. Define a difference quotient
A difference quotient is the quotient obtained by dividing the difference between
two values of a function by the difference between two corresponding values of
the independent variable
12. State Liebmann’s iteration process formula
U(n+1) i,j = (1/4)[u(n+1)i-1,j +u(n) i+1,j+ u(n)i,j-1 + u(n+1)i,j+1].
13. State the conditions for the equation, Auxx + Buxy + Cuyy + Dux + Euy +fu = G
where A, B, C, D, E, F, G are functions of x and y, to be (i) elliptic, (ii) parabolic,
(iii) hyperbolic
The given equation is said to be (i) elliptic at a point (x,y) in the plane if B2 – 4ac
< 0(ii) parabolic if B2 – 4ac = 0 (iii) hyperbolic if B2 – 4ac>0
14. At which region, the equation yuxx + uyy = 0 is hyperbolic?
Here A = y, B = 0, C = 1
B2 – 4ac = 0 – 4y = -4y
The equation is hyperbolic in the region (x, y) where B2 – 4ac > 0 (ie) –4y >0
⇒y<0
⇒ it is hyperbolic in the region y<0
15. What is the classification of fx – fyy = 0?
Here A = 0, B = 0, C = -1. B2 – 4ac = 0 –4x0x-1 = 0. So the equation is parabolic
16. Give an example of a parabolic equation
The one dimensional heat equation ∂ u/∂t = α2∂2u/∂x2, is parabolic where α2 =
(1/a)
17. Write down the finite difference form of the equation ∇2u =f(x,y). [Poisson
equation]
(or) Write down the five point formula t6o solve the Poisson equation uxx + uyy =
f(x,y)
ui-1,j +ui,j-1 + ui,j+1 – 4ui,j = h2 f(ih,jh)
The solution value at any point (I,j+1) on the (j+1)th level is expressed in terms of
the `solution values at the points (I-1,j) on the previous jth level, (not on the same
level). Hence, this formula is called explicit formulas.
24. What is the value of k to solve ðu/ðt = (½ ) uxx by Bender – Schmidt method
with h= 1 if h & k are the increments of x and t respectively?
Given ðu/ðt = (1/2) uxx
Ie) uxx = 2ut
Here a = 2, h = 1
We know that,λ = k/(ah2) = ½ [in Bender Schmidt]
2
Ie) k/(ah ) =1/2
Ie)K/2x1 = ½
Ie)K = 1
25. What is the classification of one dimensional heat flow equation?
One dimensional heat flow equation is ∂ 2u/∂x2 = a ∂u/∂t
Here A = 1, B = 0, C = 0. B2 – 4ac = 0 ⇒ It is parabolic
The solution value at any point (I,j+1) on the (j+1)th level is dependent on the
solution values at the neighboring points on the same level and on three values on
the jth level. Hence it is an implicit formula.
29. What type of equations can be solved by using Crank – Nicholson’s difference
formula?
It is used to solve parabolic equation of the form, uxx = aut [one dimensional
heat equation]
30. Express the implicit form of Crank – Nicholson formula to solve uxx = aut
Ui,j+1 = (1/4)[ui-1,j+1 + ui+1,j+1 + ui-1,j + ui+1,j] where λ = 1 ie)k = ah2
31. Write a note on the stability and of the solution of the difference equation
corresponding to the hyperbolic equation Utt = a2Uxx
a. For λ = 1/a, the solution of the differential equation is stable and coincides
with the solution of the differential equation
b. For λ>1/a, the solution is unstable.
c. forλ<1/a, the solution is stable but not convergent
32. Write down one-dimensional wave equation and its boundary conditions. Utt =a2
uxx
With boundary conditions
U (0,t) =0
U (l,t) =0
And the initial conditions
U(x,0) =f(x)
Ut(x,0) =0.
33. State the explicit scheme formula for the solution of the wave equation , utt =a2
uxx.
Ui,j+1 =2(1-λ2a2)ui,j+λ2a2(ui-1,j +ui+1,j)-ui,j-1, where λ=(k/h).
34.State the explicit formula for the one dimensional wave equation with 1-λ2α2 =0
where
λ=(k/h) & α2=(T/m)
The solution value at any point (I,j+1) on the(j+1)th level is expressed in terms of
solution values on the previous j and (j-1) levels (and not in terms of values on the
same level). Hence this is an explicit difference formula.
37. In one dimensional heat equation,∂u/∂t = α2∂2u/∂x2
Where α= k/ρc, what do the values kρc stand for?
K = thermal conductivity
ρ = Density
c = specific heat