Numerical Methods Bank
Numerical Methods Bank
com
www.rejinpaul.com
NUMERICAL METHIODS
www.rejinpaul.com www.rejinpaul.com
2 1 0 4. Find all the eigen values of the matrix 1 2 1 by power method (Apply 0 1 2 only 3 iterations).
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
Unit II : Interpolation and Approxiamtion Part A
www.rejinpaul.com
1. Using Lagranges interpolation formula find y(10) given that y(5) = 12, y(6) = 13, y(9) = 14 and y(11) = 16. 2. Find the missing term in the following table x:0 1 2 3 4 y:1 3 9 - 81
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
3. From the data given below find the number of students whose weight is between 60 to 70. Wt (x) : 0-40 40-60 60-80 80-100 100-120
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
No of students : 250 120 100 70 50
www.rejinpaul.com www.rejinpaul.com
6. Given log 10 654 = 2.8156, log 10 658 = 2.8182, log 10 659 = 2.8189, log 10 661 = 2.8202, find using Lagranges formula the value of log 10 656. 7. Fit a Lagrangian interpolating polynomial y = f(x) and find f(5) x:1 3 4 6 y : -3 0 30 132 8. Find y(12) using Newton forward interpolation formula given x : 10 20 30 40 50 y : 46 66 81 93 101 5. Given sin 450 = 0.7071, sin 500 = 0.7660, sin 550 = 0.8192, sin 600 = 0.8660, find sin 520 using Newtons forward interpolating formula.
4. From the following table find y(1.5) and y(1) using cubic spline. X : 1 2 3 Y : -8 -1 18
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
Unit III Differentiation and Integration
www.rejinpaul.com
Part A
www.rejinpaul.com www.rejinpaul.com
6. State two point Gaussian quadratue formula to evaluate
1
1. What the errors in Trapezoidal and Simpsons rule. 2. Write Simpsons 3/8 rule assuming 3n intervals. 1 dx 3. Evaluate using Gaussian quadrature with two points. 1 x 4 1 4. In Numerical integration what should be the number of intervals to apply Trapezoidal, Simpsons 1/3 and Simpsons 3/8. 1 x 2 dx 5. Evaluate using Gaussian three point quadrature formula. 1 x 4 1
f ( x)dx .
www.rejinpaul.com www.rejinpaul.com
10. Using trapezoidal rule evaluate sin xdx by dividing into six equal parts.
0
7. Using Newton backward difference write the formula for first and second order derivatives at the end value x = x0 upto fourth order. d2y dy 8. Write down the expression for and at x = x0 using Newtons forward dx dx 2 difference formula. 9. State Simpsons 1/3 and Simpsons 3/8 formula.
www.rejinpaul.com
Part B
1. Using Newtons backward difference formula construct an interpolating polynomial of degree three and hence find f(-1/3) given f(-0.75) = - 0.07181250, f(-0.5) = - 0.024750, f(-0.25) = 0.33493750, f(0) = 1.10100. 2. Evaluate
1 x y
dxdy
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
2 2
5. Evaluate 6.
1 1
www.rejinpaul.com
7.
www.rejinpaul.com
8.
www.rejinpaul.com
9.
direction. Find the value of sec (31) from the following data 32 33 34 (deg ree) : 31 Tan : 0.6008 0.6249 0.6494 0.6745 Find the value of x for which f(x) is maxima in the range of x given the following table, find also maximum value of f(x). X: 9 10 11 12 13 14 Y : 1330 1340 1320 1250 1120 930 The following data gives the velocity of a particle for 20 seconds at an interval of five seconds. Find initial acceleration using the data given below Time(secs) : 0 5 10 15 20 Velocity(m/sec): 0 3 14 69 228 7 dx Evaluate using Gaussian quadrature with 3 points. 2 1 x 3
www.rejinpaul.com
d2y dy 10. For a given data find and at x = 1.1 dx dx 2 X : 1.0 1.1 1.2 1.3 1.4 Y: 7.989 8.403 8.781 9.129 9.451
1.5 9.750 1.6 10.031
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
UNIT IV : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS PART A
dy y x2 dx
is
satisfied
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
18. Given dy x 3 y, y(0) 2, y(0.2) 2.443214, y(0.4) 2.990578, y(0.6) 3.823516 dx find y(0.8) by Milnes predictor-corrector method taking h = 0.2. 19. Using R.K.Method of order 4, find y for x = 0.1, 0.2, 0.3 given that dy xy y 2 , y (0) 1 also find the solution at x = 0.4 using Milnes method. dx dy 20. Solve y x 2 , y(0) = 1. dx Find y(0.1) and y(0.2) by R.K.Method of order 4. Find y(0.3) by Eulers method. Find y(0.4) by Milnes predictor-corrector method. 21. Solve y 0.1(1 y 2 ) y y 0 subject to y(0) 0, y(0) 1 using fourth order Runge-Kutta Method. Find y(0.2) and y (0.2) . Using step size x 0.2 . xy 22. Using 4th order RK Method compute y for x = 0.1 given y given y(0) = 1 x2 1 taking h=0.1. dy 23. Determine the value of y(0.4) using Milnes method given xy y 2 , y(0) 1 , dx use Taylors series to get the value of y at x = 0.1, Eulers method for y at x = 0.2 and RK 4th order method for y at x=0.3. dy 24. Consider the IVP y x 2 1, y(0) 0.5 dx (i) Using the modified Euler method, find y(0.2). (ii) Using R.K.Method of order 4, find y(0.4) and y(0.6). (iii) Using Adam- Bashforth predictor corrector method, find y(0.8). 25. Consider the second order IVP y 2 y 2 y e 2t S int, with y(0) = -0.4 and y(0)=-0.6. (i) Using Taylor series approximation, find y(0.1). (ii) Using R.K.Method of order 4, find y(0.2).
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
UNIT-5 PART-A
www.rejinpaul.com
2. Write down the standard five point formula used in solving laplace equation U xx + U yy = 0 at the point ( ix, jy ). 3. Derive Crank-Niclson scheme.
www.rejinpaul.com
4. State Bender Schmidts explicit formula for solving heat flow equations 5. Classify x 2 f xx + (1-y 2 ) f yy = 0
www.rejinpaul.com
6. What is the truncation error of the central difference approximation of y ' (x)?
www.rejinpaul.com
7. What is the error for solving Laplace and Poisssons equation by finite difference method. 8. Obtain the finite difference scheme fore the difference equations 2
www.rejinpaul.com
9. Write dowm the implicit formula to solve the one dimensional heat equation. 10. Define the diagonal five point formula .
d2y + y = 5. dx 2
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
PART B
www.rejinpaul.com
u(1,t) =0 using Crank- Nicholson method taking h = 1/3 k = 1/36(do on time step) 2. Solve U xx + U yy = 0 for the following square mesh with boundary values 1 2
www.rejinpaul.com www.rejinpaul.com
u1 u2 u4 1 2 4 u3 5
www.rejinpaul.com
4 5
www.rejinpaul.com
3. Solve U xx = U tt with boundary condition u(0,t) = u(4,t) and the initial condition u t (x,0) = 0 , u(x,0)=x(4-x) taking h =1, k = (solve one period)
www.rejinpaul.com
5. Solve the boundary value problem xy II -2y + x = 0, subject to y(2) = 0 =y(3).Find y(2.25),y(2.5),y(2.75).
www.rejinpaul.com
6 . Solve the vibration problem
www.rejinpaul.com
y(0,t)=0,y(8,0)=0 and y(x,0)=
1 up 2
1 x(8-x).Find y at x=0,2,4,6.Choosing x = 2, t = 2
www.rejinpaul.com
compute to 4 time steps. 7. Solve 2 u = -4(x + y) in the region given 0 x 4, 0 y 4. With all boundaries kept
www.rejinpaul.com
www.rejinpaul.com
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com
at 0 0 and choosing x = y = 1.Start with zero vector and do 4 Gauss- Seidal iteration.
www.rejinpaul.com
00 00 00 00 00 00 00
www.rejinpaul.com
00
www.rejinpaul.com
00 00 00 00 00
www.rejinpaul.com
8. Solve u xx + u yy = 0 over the square mesh of sid following conditions . e 4 units, satisfying the
www.rejinpaul.com
u(x,0) =3x for 0 x 4 u(x, 4) = x 2 u(0,y) = 0, for 0 x 4
www.rejinpaul.com
u(4,y) = 12+y for 0 y 4
for 0 y 4
www.rejinpaul.com
9. Solve h=1.Find the values of u upto t =5.
www.rejinpaul.com
10. Solve y tt = 4y xx subject to the condition y(0,t) =0, y(2,t)=o, y(x,o) = x(2-x),
y ( x,0) 0 . Do 4steps and find the values upto 2 decimal accuracy. t
www.rejinpaul.com www.rejinpaul.com
www.rejinpaul.com