Mth603 Finalterm Subjective: Question No: 50 (Marks: 5)
Mth603 Finalterm Subjective: Question No: 50 (Marks: 5)
com
Answer:
f ( x ) =sin x =1 + x 3
f ( x ) =sin x - 1 - x 3
f (1) =sin1 - 1 - 13 =0.0174 - 1 - 1 =- 1.9826
f (2) =sin 2 - 1 - 23 =0.0348 - 1 - 8 =- 8.9652
f ¢( 1) =- cos x - 3 x 2 =- cos1 - 3(1) 2 =- 0.9998 - 3 =- 3.9998
f ¢( 2 ) =- cos x - 3 x 2 =- cos 2 - 3(2) 2 =- 0.9993 - 12 =- 12.9993
f ( x0 )
x1 = x0 -
f ¢( x0 )
- 8.9652
x1 =- 3.9998 - =- 3.9998 - 0.6896 =- 4.6894
- 12.9993
Answer:
X 0 0.2 0.4 0.6 0.8 1
www.vustudy.com
www.vustudy.com
h
I=
2
( y0 + y5 + 2 ( y1 + y2 + y3 + y4 ) )
0.2
I= ( 1.0 + 2 + 2(1.2 + 1.4 + 1.6 + 1.8) )
2
0.2
I= ( 3 + 12 )
2
3
I = =1.5
2
5
Q: Evaluate the integral sin x cos x dx By Simpson`s 3/8 rule. Take the
3
interval h=1.
Answer:
X 3 4 5
Y=(sin x+ cos x) 1.0509 1.0673 1.0833
5
3h
sin x cos x dx 8 y
3
0 y2 3 y1
5
3 1
sin x cos x dx
3
8
1.0509 1.0833 3(1.0673)
5
3
sin x cos x dx 8 2.1342 3.2019
3
5
3 16.0083
sin x cos x dx 8 5.3361
3
8
2.0010375
Answer:
dy x y
, y (0) 1,
dx x y
find y (0.1)
x y 0 1 1
y0 1
x y 0 1 1
now, u sin g Pr edictor formula
4h
y1 y0 2 y0
3
Q: yʹ =2t+3y then find next two derivatives w.r.t t or y.
www.vustudy.com
www.vustudy.com
Answer:
X 0 3 4
Y 4 7 9
0 -1 1/3 y = 0
0 0 -1/3 z 1
y1 y0 hf x0 , y0
1 1
y1 1 1 1 1 1 0.5 1.5
x0 y0 11
4
Q: Evaluate the integral by Simpson’s 1/3 rule ,
sin x 1 dx taking interval
0
p
of 4 .
Answer:
X 0 p
4
Y=(sin x+ 1) 1 1.7071
www.vustudy.com
www.vustudy.com
5
h
sin x 1 dx 3 y
3
0 y1
5
3 sin x 1 dx 43 1 1.7071
5
3.14 8.500294
sin x 1 dx
3
12
2.7071
12
0.7083
2012 paper1:
Q.1. State the sufficient condition of convergence of the iterative solution to the
exact solution? 2m
Answer: page 72
A sufficient condition for convergence of the iterative solution to the exact solution is
n
aii aii , i=1, 2...n When this condition (diagonal dominance) is true, Jacobi’s
j 1
j 1
method converges
5
Take h=1…..2m
Answer:
X 3 4 5
www.vustudy.com
www.vustudy.com
Q5. Obtain numerically the solution of y` = x2 +2x+y2, y (0) =1 Using Euler’s method
to find y at x=1, h=1 …. 3m
Answer:
Y1 = y0 + hf(x0, y0)
Y1 = 1+ (1) (1+1)
Y1 = 1+2=3
Q6. If f(2)=-2.6146 and f(3)=4.7610, then find the first approximation using the
Regula-Falsi method? 3m
Answer:
x2 - x1
x3 = x2 - f ( x2 )
f ( x2 ) - f ( x1 )
3- 2
x3 =3 - ( 4.7610 )
4.7610 - ( - 2.6146 )
4.7610
x3 =3 -
7.3756
4.7610
x3 =3 - =3 - 0.6455 =2.3545
7.3756
Q7. If in a solving a giving differential equations
y0 1, y1/ 1.1, y2 / 1.2, y3/ 1.3, h 1
……3m
then find y4 by Mi ln e Pr edictor formula
Answer:
www.vustudy.com
www.vustudy.com
4h
P : yn +1 = yn - 3 + ( 2 yn¢- 2 - yn¢- 1 + 2 yn¢)
3
4h
P : y4 = y0 + ( 2 y1¢- y2¢ + 2 y3¢)
3
4(1)
P : y4 =1 + ( 2(1.1) - 1.2 + 2(1.3) )
3
4
P : y4 =1 + ( 2.2 - 1.2 + 2.6 )
3
4
P : y4 =1 + ( 3.6 )
3
14.4
P : y4 =1 +
3
3 + 14.4 17.4
P : y4 = = =5.8
3 3
Q8. Find F(h), using Richardson’s extrapolation limit, while finding
1
y / (0.01) to the function y with h 0.005 …………..3m
x
Answer:
y ( x + h) - y ( x - h)
F ( h) =
2h
1 1
-
F (h) = 0.01 + 0.005 0.01 - 0.005
2(0.005)
1 1
-
F (h) = 0.015 0.005
0.01
66.66 - 200
F ( h) =
0.01
- 133.34
F ( h) = =- 13334
0.01
Q.9
Evaluate the int egral
3
2
ò( x
0
+ x)dx
3
u sin g Simpson ' s rule
8 take h=1………5m
Answer:
X 0 1 2 3
2
Y=(x +x) 0 2 6 12
www.vustudy.com
www.vustudy.com
3
2 3
ò( x
0
+ x )dx = h [ y0 + y3 + 3( y1 + y2 ) ]
8
3
3(1)
ò( x
2
+ x )dx = [ 0 + 12 + 3(2 + 6) ]
0
8
3
3 3 108
ò( x
2
+ x )dx = [ 12 + 24 ] = ( 36 ) = =13.5
0
8 8 8
Q10. Construct a backward difference table from the following values of x and
y…..5m
x -1 0 1 2 3
Y=f(x) 10 2 10 62 80
Answer:
x Y=f(x) Ñ Ñ 2y Ñ3y Ñ4y
y
-1 10
0 2 -8
1 10 8 16
2 62 52 44 28
3 80 18 -34 -78 -106
Q11. From the following table of values, construct forward difference table…..5m
X 1.00 1.05 1.10 1.15 1.20
y 1.000 1.0247 1.0488 1.0724 1.0954
Answer:
x Y=f(x) D D2y D3y D4 y
y
1.00 1.000
1.05 1.0247 0.0247
1.10 1.0488 0.0241 0.0006
1.15 1.0724 0.0236 -0.0005 -0.0011
1.20 1.0954 0.023 -0.0006 -0.0001 0.0000
Q12. Use Runge-Kutta Method of order four to find the values of k1k2k3k4 for the
initial value problem
y/= ½(2x3 +y),y(1)=2 taking h=0.1…………….5m
Answer:
www.vustudy.com
www.vustudy.com
1
here f x, y
2
2 x 3 y , h 0.1, x1 1, y1 2
1
K1 hf x1 , y1 0.1 2 13 2
2
6 0.6
K1 hf x1 , y1 0.1 0.3
2 2
h k 0.1 0.3
K 2 hf x1 , y1 1 0.1 1 ,2
2 2 2 2
h k
K 2 hf x1 , y1 1 0.1 1.05, 2 0.15
2 2
h k
K 2 hf x1 , y1 1 0.1 1.05,2.15 0.1 1.05 2.15 0.32
2 2
h k 0.1 0.32
K 3 hf x1 , y1 2 0.1 1 ,2
2 2 2 2
h k
K 3 hf x1 , y1 2 0.1 1.05,2 0.16
2 2
h k
K 3 hf x1 , y1 2 0.1 1.05,2.16 0.1 1.05 2.16 0.321
2 2
K 4 hf x1 h, y1 k3 0.1 1 0.1, 2 0.321
K 4 hf x1 h, y1 k3 0.1 1.1, 2.321 0.1 0.1 2.21575 0.3421
Now,
1
y2 y1 k1 2k2 2k3 k4
6
1
y2 2 0.3 2 0.32 2 0.321 0.3421
6
1 1.9241
y2 2 0.3 0.64 0.642 0.3421 2 2 0.320683 2.320683
6 6
2012 Paper2
Q1. Obtain numerically the solution of y/=x+y2,y(0)=2 Using Euler`s method to find
y at x=1, h=1….2m
Answer:
h =1, x =1 and y (0) =2
y = y0 + hf ( x0 , y0 )
y =2 + 1 ( 1 + 2 )
y =2 + 3 =5
www.vustudy.com
www.vustudy.com
Q2. If f(0)=3 and f(1)=9, then the next approximate value of the function using
secant method….2m
Answer:
x f ( x1 ) x1 f ( x0 )
x2 0
f ( x1 ) f ( x0 )
(0)9 (1)3
x2
93
0 3 3
x2 0.5
93 6
Q3. Write a backward difference formula of D2f(x)….2m
Answer:
Q5.
Evaluate the int egral
5
(log x x )dx
3
…..3m
3
u sin g Simpson ' s rule
8
take h 1
Answer:
X 3 4 5
Y=(log x+ x) 3.477 4.602 5.6989
1 0
5
3
3 (log x x)dx 8 h y0 y2 3 y1
5
3
(log x x)dx 8 (1) 3.4771 5.6989 3(4.6020)
3
5
3
(log x x)dx 8 9.176 13.806
3
5
3 68.946
(log x x)dx 8 (22.982)
3
8
8.61825
www.vustudy.com
www.vustudy.com
Q6. Write the two steps of solving the linear equations using Gaussian Elimination
method….3m
Answer:
In this method, the solution to the system of equations is obtained in two stages.
I) the given system of equations is reduced to an equivalent upper triangular form using
elementary transformations
ii) The upper triangular system is solved using back substitution procedure
Q.7
Use Runge Kutta Method of order four to find the values of
k1 and k2 for the int ial value problem ….3m
y yx , y (0) 1, h 0.1
/ 2
Answer:
f ( x, y ) yx 2 , y (0) 1 and h 0.1, x0 0
k1 hf ( x0 , y0 ) 0.1 1 02 0
h k 0.1 0
k2 hf x0 , y0 1 0.1 0 ,1 0.1 0.05,1 0.1 0.05 1 0.105
2 2 2 2
1
y1 y0 k1 2k2
6
1 0.21
y1 1 0 2 0.105 1 1.035
6 6
Q.8
Evaluate the int egral
3
(x x )dx
2
0 …..3m
U sin g T apezoidal rule
take h 1
Asnwer:
x 0 1 2 3
Y=(x2 + x) 0 2 6 12
www.vustudy.com
www.vustudy.com
3
h
(x
2
x )dx y0 y3 2( y1 y2 )
0
2
3
1
(x
2
x )dx 0 12 2(2 6)
0
2
3
1
(x
2
x )dx 28
0
2
3
28
(x x )dx 14
2
0
2
Q.9
Evaluate the int egral
3
( x 1)dx
0
…..5m
3
U si ng Simpson ' s rule
8
Answer:
x 0 1 2 3
Y=(x + 1) 1 2 3 4
3
3
0 ( x 1)dx 8 h y0 y3 3( y1 y2 )
b a 3 0 3
h 0.3
3n 3(3) 9
3
3
(x
2
x )dx (0.3) 1 4 3(2 3)
0
8
3
0.9
(x
2
x )dx 5 15
0
8
3
0.9 18
(x x )dx (20) 2.25
2
0
8 8
Q10. Find the dominant Eigen-value and the corresponding eigenvector of the
matrix
é8 1 2ù
ê0 10 - 1úby Power Method with vector with vector v 0 =(1,1,1)T as theinitial vector
ê ú
ê
ë6 2 15 ûú ….5
m
Answer:
www.vustudy.com
www.vustudy.com
vector v 0 (1,1,1)T as theinitial vector. Then, compute first iteration
8 1 2 1 8 1 2 1 11
u A v 0 10 1 1 0 10 1 1 9
1 0
6 2 15 1 6 2 15 1 23
11
23
9
u1 23 q1v1
23
1
sec ond iteration
11 88 9 143
23 23 23 2 23
8 1 2
1 9 90 67
u A v 0 10 1
2
0 1
23 23 23
6 2 15
1 66 18 429
15
23 23 23
143
u 23 67 q2v 2
2
429
third iteration
8 1 2 143 8 1 2 143 1573
u A v 0 10 1 67 0 10 1 67 603
3 2
6 2 15 429 6 2 15 429 9867
1573
9867 0.0159
0.0159
603
u 9867
3 0.0611 9867 0.0611 q3v 2
9867
1.0 1.0
The l arg est eign value and the corresponding eignvector as accurate to two decimals are
0.01
9867 and X 0.06
1.0
Q11. From the following table of values, construct backward difference table….5m
X 0.1 0.2 0.3 0.4
F(x) 1.10517 1.22140 1.34986 1.49182
www.vustudy.com
www.vustudy.com
Answer:
x F(x) Ñ f (x) Ñ 2 f ( x) Ñ 3 f ( x)
0.1 1.10517
0.2 1.22140 0.11623
0.3 1.34986 0.12846 0.01223
0.4 1.49182 0.14196 0.0135 0.00127
2 Marks Questions:
Q. Evaluate the integral
/2
Answer:
x 0 1 2 3 4
Y= cos 2x 0 0.999 0.9975 0.9945 0.9902
3
www.vustudy.com
www.vustudy.com
2
h
cos 2 xdx 3 y
0
0 y4 4( y1 y3 ) 2 y2
2
3 Marks Questions:
Q. Evaluate the integral
Answer:
www.vustudy.com
www.vustudy.com
X p 3p p
2 4
Y=sinx 1.0 0.7071 0.0
3h
sin xdx 8 y 0 y2 3 y1
2
3
4
sin xdx 8 1.0 0.0 3(0.7071)
2
3
Answer:
here f ( x, y ) =1 + xy , h =0.2, x0 =0, y0 =2
K1 =hf ( x0 , y0 ) =0.2 ( 1 + ( 0 ) ( 2 ) )
K1 =hf ( x0 , y0 ) =0.2 ( 1 + 0 ) =0.2 ( 1) =0.2
æ h k ö æ 0.2 0.2 ö
K 2 =hf çx0 + , y0 + 1 ÷=0.2 ç0 + ,2 + ÷
è 2 2ø è 2 2 ø
æ h k ö
K 2 =hf çx0 + , y0 + 1 ÷=0.2 ( 0.1, 2 + 0.1)
è 2 2ø
æ h k ö
K 2 =hf çx0 + , y0 + 1 ÷=0.2 ( 0.1, 2.1) =0.2 ( 0.1 + 2.1) =0.44
è 2 2ø
Now,
1
y1 = y0 + ( k1 + 2k2 )
6
1
y1 =2 + ( 0.2 + 2 ( 0.44 ) )
6
1 1.08
y1 =2 + ( 0.2 + 0.88 ) =2 + =2 + 0.18 =2.18
6 6
www.vustudy.com
www.vustudy.com
Q. Find the residuals by relation method
6 x1 - 3 x2 + x3 =11
2 x1 + x2 - 8 x3 =- 15
x1 - 7 x2 + x3 =10
With starting vector (0,0,0)
Answer:
Re-arrange the equation:
6x1-3x2+x3=11
X1-7x2+x3=10
2x1+x2-8x3=-15
Thus, we get
R1=11-6x1+3x2-x3=11-6(0) +3(0)-(0) =11
R2=10-x1+7x2-x3=10-(0) +7(0)-(0) =10
R3=-15-2x1-3x2+8x3=-15-2(0)-3(0)+8(0)=-15
The largest residual in magnitude is R3
R3 - 15
- =- =1.875
a33 8
Dx3=
Similarly, we find the new residuals of large magnitude and relax it to zero,
and so on.
5 Marks Questions:
Q. Evaluate the integral
3
2
ò( x + 1) dx
0 h =1
Using Simpson’s 3/8 rule. Where
Answer:
x 0 1 2 3
Y=(x2 + 1) 1 2 5 10
www.vustudy.com
www.vustudy.com
3
3
(x
2
1)dx h y0 y3 3( y1 y2 )
0
8
3
3 1
( x 1)dx
2
1 10 3(2 5)
0
8
3
3
(x
2
1)dx 11 21
0
8
3
3 96
(x 1)dx (32) 12
2
0
8 8
www.vustudy.com
www.vustudy.com
1
x = [ 8 + y + 2z]
8
1
y =- [ - 4 - x - 2 z ]
7
1
z = [ 12 - 2 x - y ]
9
initial approximantion is (0, 0, 0)
first iteration
1 8
x1 = [ 8 + 0 + 2(0) ] = =1
8 8
1 -5
y1 =- [ - 4 - 1 - 2(0)] =- =0.71
7 7
1 9.29
z1 = [ 12 - 2(1) - 0.71] = =1.032
9 9
sec ond iteration
1 1 10.774
x2 = [ 8 + 0.71 + 2(1.032) ] = [ 8 + 0.71 + 2.064 ] = =1.34675
8 8 8
1 1 - 7.41075
y2 =- [ - 4 - 1.34675 - 2(1.032) ] =- [ - 4 - 1.34675 - 2.064 ] =- =1.0586
7 7 7
1 1 8.2479
z2 = [ 12 - 2(1.34675) - 1.0586 ] = [ 12 - 2.6935 - 1.0586 ] = =0.91643
9 9 9
Answer:
x 0 1 2 3 4
Y=(x2 + x+2) 2 4 8 14 22
4
h
(x
2
x 2)dx
3
y0 y4 4( y1 y3 ) 2 y2
0
4
1
(x
2
x 2)dx
3
2 22 4(4 14) 2 8
0
4
1
(x
2
x 2)dx 24 72 16
0
3
4
1 112
(x
2
x 2)dx 112 37.33
0
3 3
www.vustudy.com
www.vustudy.com
Q. Find the 2nd derivative of f(x) at x=0.3 using three point equation
x 0.1 0.2 0.3 0.4 0.5 0.6
f(x) 0.125 0.352 0.652 0.756 0.812 0.924
Answer:
f ( x3 ) - 2 f ( x2 ) + f ( x1 )
f ¢¢( 0.3 ) =
h 2 =0.01
f (0.4) - 2 f ( 0.3) + f ( 0.2 )
f ¢¢( 0.3 ) =
h 2 =0.01
0.756 - 2 ( 0.652 ) + 0.352
f ¢¢( 0.3 ) =
h 2 =0.01
0.756 - 1.304 + 0.352 0.196
f ¢¢( 0.3 ) = 2
=- =- 19.6
h =0.01 0.01
Final 2013:
Q. Find Newton`s forward difference table from the following data
…5m
Answer:
X Y=f(x) Dy D2 y D3y D4 y
0.0 1
www.vustudy.com
www.vustudy.com
x - x0 x - 0.0 x 0.0 x x
p= = = - = - 0=
h 0.1 0.1 0.1 0.1 0.1
p ( p - 1) 2 p ( p - 1) ( p - 2) 3 p ( p - 1) ( p - 2)( p - 3) 4
y = y0 + pD y0 + D y0 + D y1 + D y2
2! 3! 4!
Dy0 = -0.0952,
x0 =0.0, y0 =1, =- = D 2 y =-= 0.091, D= 3
y = -0.009, D 4 y = 0.001
x æ x ö x æ x öæ x ö x æ x öæ x öæ x ö
ç - 1÷ ç - 1÷ç - 2÷ ç - 1÷ç - 2 ÷ç - 3÷
æx ö 0.1 è 0.1 ø
y =1 + ç ÷( - 0.0952 ) + ( 0.091) + 0.1 è 0.1 øè 0.1 ø - 0.009 + 0.1 è 0.1 øè 0.1
( ) øè 0.1 ø 0.001
( )
è0.1 ø 2 6 24
this is required Newton`s int erpolating polynomial.
Q5.
Evaluate the int egral
5
3
ò(log x + 2)dx u sin g Simpson ' s 8 rule. take h =1
3
…..3m
Answer:
X 3 4 5
Y=(log x+ 2) 2.477 2.602 2.6989
1 0
5
3
ò3 (log x + 2)dx =8 h ( y0 + y2 + 3 y1 )
5
3
ò(log x + 2)dx =8 (1) ( 2.4771 + 2.6989 + 3(2.6020) )
3
5
3
ò(log x + 2)dx =8 ( 5.176 + 7.806 )
3
5
3 38.946
ò(log x + 2)dx =8 (12.982) =
3
8
=4.86825
Q5.
Evaluate the int egral
2
3 …..2m
0 (cos x 2)dx u sin g Simpson ' s
8
rule. take h
4
www.vustudy.com
www.vustudy.com
Answer:
X 0 p p
4 2
Y=(cos x+ 2) 3 2.707 2
1
2
3h
(cos x 2)dx 8 y
0
0 y2 3 y1
2 3
4
0 (cos x 2)dx 8 3 2 3(2.7071)
3
2
æh ö h
Fç ÷ F1
è2 ø
Q. if F (h) =256.2354 and =257.1379, then find 2 using Richardson`s
extrapolation limit….2m
Answer:
h h
4f ( )- f ( )
h 2 2
F1 ( ) =
2 41 - 1
h 4(257.1379) - (256.2354)
F1 ( ) =
2 3
h 1028.5516 - 256.2354
F1 ( ) =
2 3
h 772.3162
F1 ( ) = =257.43873
2 3
7. if in solving given differential equation, 0.1, ,,, h=1, y3=2 then find the y4 by Adam–
moulton’s predictor formula .
Answer:
www.vustudy.com
www.vustudy.com
h
P : yn+1 = yn + ( 25 yn¢ - 59 yn¢- 1 + 37 yn¢- 2 - 9 yn¢- 3 )
24
1
P : y3+1=4 = y3 + ( 25 y3¢ - 59 y3¢- 1=2 + 37 y3¢- 2=1 - 9 y3¢- 3=0 )
24
1
P : y4 = y3 + ( 25 y3¢ - 59 y2¢ + 37 y1¢- 9 y0¢)
24
1
P : y4 =2 + ( 25 ( 1.3) - 59 ( 1.2 ) + 37 ( 1.1) - 9 ( 0.1) )
24
1
P : y4 =2 + ( 32.5 - 70.8 + 40.7 - 0.9 )
24
1 1.5
P : y4 =2 + ( 1.5 ) =2 + =2 + 0.0625 =2.0625
24 24
12. using Adam-moulton’s predictor formula find the predicted value of f(0.4) from
ordinary differential equation y’=1+2xy; y(0)=0; h=0.1
x 0 0.1 0.2 0.3
Answer:
www.vustudy.com
www.vustudy.com
www.vustudy.com
www.vustudy.com
Q. using Adam-moulton’s predictor formula find the predicted value of f(0.4) from
ordinary differential equation y’=1+xy; y(0)=0; h=0.1
x 0 0.1 0.2 0.3
Answer:
www.vustudy.com
www.vustudy.com
www.vustudy.com
www.vustudy.com
2
Q.2. Evaluate the integral cos xdx Using trapezoidal rule
0
p
Take h= 4 …..2m
Answer:
X 0 p
4
Y= cosx 1 0.7071
p
Using trapezoidal rule, taking h= 4
www.vustudy.com
www.vustudy.com
h
I= ( y0 + y1 )
2
p
I = 4 ( 1 + 0.7071)
2
3.14
I= ( 1.7071)
2
5.360294
I= =2.680147
2
Q. obtain numerically the solution of
y x y 2 , y (0) 2 u sin g Euler `s method to find y at x 1, h 1.
Answer:
Here: y(0)=2, x=1, h=1
y y0 hf x0 , y0
y 2 1 1 4 2 5 7
Q5.
p
3 p
Evaluate the int egral òsin xdx u sin g Simpson ' s rule. take h =
p 8 4
2
…..2m
Answer:
X p 3p p
2 4
Y=sinx 1 0.707 0
1
3h
sin xdx 8 y 0 y2 3 y1
2
3
4
sin xdx 8 1 0 3(0.7071)
2
3
www.vustudy.com
www.vustudy.com
1
Q. Prove that E 2
2
Answer:
1
E
2
2
1 1
1
1 1
1
E2 E 2
E E2 E 2
2
2 2
1 1 1 1
1
E2 E 2
E2 E 2
E
2
2 2
1 1 1 1 1
1
E E
2 2
E E2 2
2E 1 2
E
2
E2
2 2
2
Take h=1…..2m
Answer:
X 0 1 2
Y= (x2-3x+4) 4 2 2
take h 1
Answer:
x 0 1 2 3 4
Y=(x2 + 2) 2 3 6 11 18
www.vustudy.com
www.vustudy.com
4
h
(x
2
2)dx
3
y0 y4 4( y1 y3 ) 2 y2
0
4
1
(x
2
2)dx
3
2 18 4(3 11) 2 6
0
4
1
(x
2
2)dx 20 56 12
0
3
4
1 88
(x
2
2)dx 88 29.33
0
3 3
3
Answer:
x 0 1 2 3
Y=x 0 1 2 3
3
3h
0 xdx 8 y0 y3 3( y1 y2 )
3
3
xdx 8 0 3 3 1 2
0
3
3
xdx 8 3 9
0
3
3 36
xdx 8 12
0
8
4.5
Final 2014
1- Write three equations for the following system of linear equation for Gauss Seidal
Approximation Method (2)
83x+11y-4z=95
7x+52y++13z=104
3x+8y+29z=71
ò (log x + 2) dx
3
www.vustudy.com
www.vustudy.com
X 3 4 5
Y=(log x+2) 2.4771 2.6020 2.6989
5
h
3 (log x 2)dx 2 y0 y2 2 y1
5
1
(log x 2)dx 2 2.4771 2.6989 2 2.6020
3
5
1
(log x 2)dx 2 2.4771 2.6989 5.204
3
5
1 10.38
(log x 2)dx 2 10.38
3
2
5.19
3- Write a formula for finding the value of k2 in fourth order R-K mrthod?(2)
Answer:
4- If y’=t2-y, then find the next two derivatives in terms of t and or/y. (2)
Answer:
y0 1, y '1 1.1, y ' 1.2, y '3 1.3,h 1 Then find y4 by Milne’s Predictor formula. (3)
Answer:
www.vustudy.com
www.vustudy.com
4h
P : yn +1 = yn - 3 + ( 2 yn¢- 2 - yn¢- 1 + 2 yn¢)
3
4h
P : y4 = y0 + ( 2 y1¢- y2¢ + 2 y3¢)
3
4(1)
P : y4 =1 + ( 2(1.1) - 1.2 + 2(1.3) )
3
4
P : y4 =1 + ( 2.2 - 1.2 + 2.6 )
3
4
P : y4 =1 + ( 3.6 )
3
14.4
P : y4 =1 +
3
3 + 14.4 17.4
P : y4 = = =5.8
3 3
x dx
2
Y’=x2+2x+y2, y(6)=1
Using Euler’s method to find y at x=1, h=1. (3)
Answer:
www.vustudy.com
www.vustudy.com
y (6) =1, x6 =1, h =1
y7 = y6 + hf ( x6 , y6 )
y7 =1 +1( x 2 + 2 x + y 2 )
y7 =1 +1( 1 + 2 +1) =1 + 4 =5
9- Using predictor formula find the predicted value of y(0,4) from ordinary differential
equation
www.vustudy.com
www.vustudy.com
F(x)=ex-3x=0. (Note only 1 iteration is required upto 4 decimal places is required) (5)
Answer:
f x e x 3x 0
f 1 e1 3(1) 2.7182 3 0.2818
f 2 e 2 3(2) 7.3890 6 1.389
x2 x1
x3 x2 f ( x2 )
f ( x2 ) f ( x1 )
2 1
x3 2 (1.389)
1.389 ( 0.2818)
1.389
x3 2 2 0.8313 1.1687
1.6708
Answer:
x 0 0.5 1 1.5 2 2.5 3 3.5 4
Y=x 0 0.5 1 1.5 2 2.5 3 3.5 4
www.vustudy.com
www.vustudy.com
4
h
òxdx = 3 ( y
0
0 + y8 + 4( y1 + y3 + y5 + y7 ) + 2 ( y2 + y4 + y6 ) )
b- a 4- 0 4
h= = = =0.5
2n 2 ( 4) 8
4
0.5
òxdx = 3 ( 0 + 4 + 4(0.5 +1.5 + 2.5 + 3.5) + 2 ( 1 + 2 + 3 ) )
0
4
0.5
òxdx = 3 ( 4 + 32 + 12 )
0
4
0.5 24
òxdx = 3 ( 48 ) = 3 =8
0
12- Find y(1) for the following initial value problem y y’=x, y(0)=2
www.vustudy.com
www.vustudy.com
decimal. Suppose you have the root from last iteration 0.8324516 and
0.8324525 if we subtract both and consider absolute value of the difference
0.0000009 now it can be written as 0.09*10-5 so the TOL in this case is 10-
5.similarly if we have been provided that you have to for the TOL 10-2 you
will check in the same way. In the given equation you will solve the
equation by any method and will consider some specific TOL and try to go
to that TOL. Some time no TOL is provided and you are asked to perform to
some specific no of iterations.
Question: When the forward and backward interpolation formulae are used?
Answer: In interpolation if we have at the start then we use the forward difference
formula and the formula to calculate p is x-x0/h. If the value of x lies at the
www.vustudy.com
www.vustudy.com
end then we use Newton’s backward formula and formula to calculate the
value of p is x-xn/h. Now I come to your question as in this case the value
lies at the end so 6 will be used as the xn. This procedure has been followed
by the teacher in the lectures. But some authors also use another technique
that is if you calculate the value of p and that is negative then the origin is
shifted to that value for which the value of p becomes positive. And then
according to that origin the values of differences are used and you need not
follow this procedure.
Question: What is forward and backward difference operator and the construction of
their table.
Answer: For forward Dfr =fr+1 –fr Df0 = f1-f0 In terms of y Dyr+1=yr+1-yr D
stands for the forward difference operator For backward Dfr =fr –fr-1 Df1 =
f1-f0 In terms of y Dy1=y1-y0 Here D stands for backwards operator Now
the construction of the difference table is based on X Y 1st forward 2nd
forward 3rd forward x1 Y1 Y2-Y1=Dy0 x2 Y2 Y3-Y2=Dy1 x3 Y3 Y4-
Y3=Dy2 x4 Y4 Now consider the construction of table for the backward
table X Y 1st forward 2nd forward 3rd forward x1 Y1 Y2-Y1=Dy1 x2 Y2
Y3-Y2=Dy2 x3 Y3 Y4-Y3=Dy3 x4 Y4 Dear student this is the main
difference in the construction of the forward and backwards difference table
when you proceed for forward difference table you get in first difference the
value Dy0 but in the construction of backwards difference table in the first
difference you get Dy1 and in the second difference in the forward
difference table you get D2 y0 and in the backward difference table the first
value in the second difference is D2 y1. I think so you have made it clear.
www.vustudy.com
www.vustudy.com
0.985 y2 y3-y2/x3-x2=sy2 sy3- sy2/x4-x1=s2y2 4 1.2365 y3 y4-y3/x4-
x3=sy3 5 2.3651 y4 This is the complete table showing the central
differences . Here s stands for difference
Question: State the sufficient condition for the convergence of the system of equation
by iterative methods.
Answer: A sufficient condition for convergence of iterative solution to exact solution
www.vustudy.com
www.vustudy.com
is |a1| => (greater or equal) |b1|+|c1| |b2| => (greater or equal) |a2|+|c2| |c3|
=> (greater or equal) |a3|+|b3| For the system a1x+b1y+c1z=d1,
a2x+b2y+c2z=d2, a3x+b3y+c3z=d3 Similarly for the system with more
variables we can also construct the same condition
Question: The calculation for numerical analysis should be done in degree or radians.
Answer: All the calculation for numerical analysis should be done in radians not in
degrees set your calculator in radians mode and suppose the value of
pi=3.14.
www.vustudy.com
www.vustudy.com
more iteration and the example of open method is Newton raphson method.
In bracketing method we bracket the root and find that interval in which root
lies means we need two initial approximations for the root finding.
Bisection method is an example of the bracketing method.
Question: Should the system be diagonally dominant for gauss elimination method?
Answer: The system of equation need not to be diagonally dominant for Gauss
elimination method and gauss Jordan method for both the direct method it is
not necessary for the system to be diagonally dominant .it should be
diagonally dominant for iterative methods like jacobie and gauss seidel
method.
www.vustudy.com
www.vustudy.com
approximation are known as direct methods like Graffee’s root squaring
method.
www.vustudy.com
www.vustudy.com
computer, which is obviously in binary form, should be converted into
decimal form-a form understandable to the user-before their print out.
Therefore, an additional error is committed at this stage too. This error is
called local round-off error.
Question: Which equation models the rate of change of any quantity with respect to
another?
Answer: An ordinary differential equation models the rate of change of any quantity
with respect to another.
Question: If a system has the base M, then how many different symbols are needed to
represent an arbitrary number? Also name those symbols.
Answer: M 0, 1, 2, 3, … , M-1
Is the order of global error in Simpson’s 1/3 rule equal to the order of
global error in Simpson’s 3/8 rule?
www.vustudy.com
www.vustudy.com
Answer: Yes. The order of global error in Simpson’s 1/3 rule equal to the order of
global error in Simpson’s 3/8 rule.
Question: What is the formula for finding the width of the interval?
Answer: Width of the interval, h, is found by the formula h=(b-a)/n
Important Topics
Newton's Forward Diffrence & Backward Difference (Interpolation)
Newton's Forward Diffrence & Backward Difference (Differentiation)
Trapezoidal Rule (***Very Important***)
Simpson's 1/3 Rule
Simpson's 3/8 Rule
Euler's Method (**Important**)
Two Point Formula for Differentiation (**Important**)
Three Point Formula for Integration (**Important**)
Runge-Kutta Second Order Method (**Important**)
Runge-Kutta Fourth Order Method (***Very Important***)
Adam-Moulton Method {Predictor, Corrector and Error} (**Important**)
Milne Method {Predictor, Corrector and Error}
mcqs was mostly from forward backward interpolating and runge kutta and lagrange 1/3 simpson and
trapezoidal rule
2 marks question : adam moultan predictor corrector formula
2 marks question : milne predictor corrector formula
3 marks question : formula of k3 of runge kutta
3 marks question : 1/3 simpson question
3 marks : trapezoidal rule question
5 marks : euler related
5 marks: trapezoidal
www.vustudy.com
www.vustudy.com
5 marks: milne predictor corrector formula question in the handouts on page 211(find f(0.4) from
ordinary differencial equation )
www.vustudy.com