Ma2264-Numerical Methods: Unit-I: Solution of Equations and Eigenvalue Problems
Ma2264-Numerical Methods: Unit-I: Solution of Equations and Eigenvalue Problems
Ma2264-Numerical Methods: Unit-I: Solution of Equations and Eigenvalue Problems
1 1 2
1 2 3
3. Find the inverse of the given matrix by Gauss-Jordan method A= 2 3 1
[ ]
Dr. D. Saravanan, Professor of Mathematics
2
25 1 2
1 3 0
4. Find the numerically largest eigenvalue of A= 2 0 −4
[ ] by power method and the
corresponding eigenvector (correct to three decimal points). Start with initial eigenvector
[]
0
0
6. Determine the largest eigen value and the corresponding eigen vector of the matrix
1 3 −1
[ 3 2 4
−1 4 10 ] by power method
7. Using Newton-Raphson’s method, solve xlog10x=12.34. Start with x0=10
2 2 3
6 0 1
matrix
[ ]
0 −1 0
1 0 6
1 1 3
1 3 −3
14. Using Gauss-Jordan method, find the inverse of A= −2 −4 −4
[ ]
15. Solve the following equations by Gauss-Seidel method 4x+2y+z=14; x+5y-z=10; x+y+8z=20
16. Solve x2+y2=1; xy=x+y given that x0=0.5, y0=-1 by Newton-Raphson method (2iterations)
17. Solve x+3y+3z=16; x+4y+3z=18; x+3y+4z=19 by Gauss-Jordan method.
1 −3 2
4 4 −1
18. Find the numerically largest eigen value of A= 6 3 5 [ ] by power method.
19. Solve the system ex-y=0 and xy-ex=0 by Newton-Raphson method to obtain the solution near
(0.95, 2.7). Perform only 3 iterations.
2 0 1
PART-A
1. Obtain the interpolation quadratic polynomial for the given data by using Newton’s forward
difference formula.
x: 0 2 4 6
y: -3 5 21 45
2. Obtain the divided difference table for the following data
x: 2 3 5
y: 0 14 102
3. Find a polynomial for the following data by Newton’s backward difference formula
4. Show that
Δ
abc
( 1a )=− abcd
1
PART-B
9. Find f(x) as a polynomial in x for the following data by Newton’s divided difference formula.
x: -4 -1 0 2 5
f(x): 1245 33 5 9 1335
10. Find f(8) by Newton’s divided difference formula for the data:
x: 4 5 7 10 11 13
f(x): 48 100 294 900 1210 2028
11. Find the polynomial f(x) by using Lagrange’s formula and hence find f(3) for
x: 0 1 2 5
f(x): 2 3 12 147
12. Using Newton’s divided differences method find f(1.5) using the data f(1.0)=0.7651977,
f(1.3)=0.6200860, f(1.6)=0.4554022, f(1.9)=0.2818186 and f(2.2)=0.1103623.
13. Construct an approximating polynomial for the following data:
x f(x) f’(x)
0.8 0.22363 2.16918
1.0 0.65810 2.04670
14. Compute cosh0.56 from the data cosh0.5=1.127626, cosh0.6=1.185465,
Cosh0.7=1.255169 and cosh0.8=1.337435.
14. The following values of x and y are given:
X 1 2 3 4
Y 1 2 5 11
Find the cubic spline and evaluate y(1.5).
PART-A
dy d2 y
and 2
1. Write down the expressions for dx dx at x=x0 by Newton’s forward difference formula.
4
∫ f ( x)dx
2. Evaluate 1 from the table by Simpson’s 3/8 rule
x: 1 2 3 4
f(x): 1 8 27 64
3. What is the condition for Simpson’s 3/8 rule and state the formula.
4. Using Newton’s backward difference formula, write the formulae for the first and second order
derivatives at the end value x=xn up to the fourth order difference term.
xn
∫ ydx
x0
5. In order to evaluate by Simpson’s 1/3 rule as well as by Simpson’s 3/8 rule, what is the
restriction on the number of intervals?
6. State three point Gaussian quadrature formula.
π
∫ sin xdx
7. Using Trapezoidal rule evaluate 0 by dividing the range into 6 equal parts.
1
dx
∫ 1+x 2
8. Evaluate −1 by Gaussian two point formula.
6
∫ f ( x)dx
9. Write down the Trapezoidal rule to evaluate 1 with h=0.5
10. By differentiating Newton’s backward difference formula, find the first derivative of the function
f(x).
∫ f ( x)dx
of the form a ?
PART-B
1
∫ 1+dxx2
1. Compute 0 by using Trapezoidal rule, taking h=0.5 and h=0.25. Hence find the value of
the above integration by Romberg’s method.
1
dt
∫ 1+t
2. Evaluate I= 0 by using three point Gaussian quadrature formula.
2 2
∫ ∫ dxdy
x+ y
3. Evaluate 1 1 by using Simpson’s 1/3 rule, taking Δx=Δy =0.25.
4. Obtain the value of f1(0.04) using an appropriate 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
∫ ∫ x+1 y
5. Using Trapezoidal rule evaluate 1 1 dx dy taking 4 sub intervals
6. For the given data
X: 1.0 1.1 1.2 1.3 1.4 1.5 1.6
F(x): 7.989 8.403 8.781 9.129 9.451 9.750 10.031
dy d2 y
and 2
Find dx dx at x=1.1.
5
∫ ( 1x ) dx.
7. Use three point Gauss formula to evaluate 1
8. The following data gives the velocity of a particle for 20 seconds at an interval of 5 seconds. Find
the initial acceleration using the entire data.
1 2
2 xydxdy
∫∫
12. Evaluate 0 1 ( 1+x 2 ) ( 1+ y 2 ) by Trapezoidal rule with h=k=0.25
13. Find the value of sec310 from the following data:
(deg): 310 320 330 340
tan: 0.6008 0.6249 0.6494 0.6745
6
∫ 1+dxx2
14. Using Simpson’s 3/8th rule evaluate 0 , by dividing the range into 6 equal parts.
15. Find the maximum and minimum value of y tabulated below:
x: -2 -1 0 1 2 3 4
y: 2 -0.25 0 -0.25 2 15.75 56
16. Find f!(3) and f”(3) for the following data:
x: 3.0 3.2 3.4 3.6 3.8 4.0
f(x): -14 -10.032 -5.296 -0.256 6.672 14
∫ 4 dx
x+5
17. Evaluate 0 by Simpson’s one-third rule and hence find the value of loge5 (n=10)
2 2
∫ ∫ f ( x, y )dxdy
18. Evaluate 0 0 by Trapezoidal rule for the following data:
2
x 2 +2 x+1
∫ 4
dx
19. Evaluate 0 1+( x+1) by Gaussian three point formulae.
20. Find the first and second derivatives at x=1.6 for the function represented by the
following tabular data:
x: 1 1.5 2.0 3.0
y: 0 0.40547 0.69315 1.09861
1
dx
∫
22. Using three point Gaussian quadrature, evaluate 0 √1+x 4
1 2
23. Evaluate 0
{
∫ ∫ ( 2 ) ( 2)
1 1+ x 1+ y
2 xydy
}dx
, using Simpson’s 1/3 rule with step length h=k=0.25
Single step methods: Taylor series method - Euler methods for first
order – Runge – Kutta method for solving first and second order
equations – Multistep methods: Milne’s and Adam’s predictor and
corrector methods.
PART-A
1. State Adams-Bashforth predictor and corrector formula.
2. Find y(1.1), given dy/dx=x+y, y(1)=2 by Euler’s method.
3. What is the condition to apply Adams-Bashforth method?
4. In the derivation of fourth order Runge-Kutta formula, why it is called fourth order.
5. Using modified Euler’s method, find y(0.1) if dy/dx=x2+y2, y(0)=1
6. Write down the formula to solve 2nd order differential equation using Runge-Kutta method.
7. What are the values of k1 and l1 to solve y”+xy’+y=0; y(0), y’(0)=0 by Runge-kutta method of
fourth order.
8. Solve the differential equation dy/dx=x+y+xy, y(0)=1 by Taylor series method to get the value of y
at x=h.
13. Apply Milne’s method, to find a solution of the differential equation dy/dx=x-y2 at x=0.8, given
the values.
x: 0 0.2 0.4 0.6
PART-A
1. Obtain the finite difference scheme for the differential equation 2d2y/dx2+y=5.
2. State standard five points finite difference formula for solving uxx+uyy=0.
3. Write down the finite difference scheme for the differential equation d2y/dx2-3dy/dx=2
4. State the implicit finite difference scheme for one dimensional heat equation.
5. What is the error for solving Lap lace and Poisson’s equations by finite difference method?
6. Write down the Bendre-Schmidt recurrence relation for one dimensional heat equation
7. Write down the Crank-Nicloson formula to solve ur=uxx.
8. Write the diagonal five-point formula to solve the Laplace equation uxx+uyy=0
9. Write the down the Implicit formula to solve one dimensional heat flow equation uxx= (1/c2)ut.
∂u 1
= u xx
10. What is the value of k to solve ∂t 2 by Bender-Schmidt method with h=1 if h and k are
the increments of x and t respectively?
11. Derive the difference scheme for the one dimensional heat equation using backward difference for
the time derivative.
12. Set up a finite difference scheme for the boundary value problem u”=u, u’ (1) =a and u’ (3) =b
with h=0.5 using central differences.
PART-B
2 1
2 1
6. Solve 2u=-10(x2+y2+10) over the square mesh with sides x=0, y=0, x=3, y=3 with u=0 on the
boundary and mesh length of 1 unit.
∂u ∂2 u
2 =
7. Solve ∂ t ∂ x2 given u(0,t)=0, u(4,t)=0, u(x,0)=x(4-x), taking x=1 and t=1. Find the
value of u up to t=3 using Bender-Schmidt’s explicit finite difference scheme.
∂u ∂ 2 u
=
8. By Crank-Nicholson method solve the equation ∂t ∂ x 2 subject to u(x,0)=0, u(0,t)=0 and
u(1,t)=t for two time steps.
∂2 u ∂2 u ∂u
2
= 2
9. Solve ∂t ∂ x , 0<x<1, t>0 given u(x,0)=0, ∂t (x,0)=0, u(0,t)=0 and u(1,t)=100sint.
Compute u(x, t) for 4 times steps with h=0.25
10. Using finite difference method, find y(0.25), y(0.5) and y(0.75) satisfying the differential
d2 y
+ y =x
equation dx 2 , subject to the boundary conditions y(0)=0, y(1)=2.
∂u ∂ 2 u
=
11. Solve ∂t ∂ x 2 , 0x1, t0 with u(x,0)=x(1-x), 0<x<1 and u(0,t)=u(1,t)=0 t>0 using
1 , 0≤x≤1 /2
t>0, u(x,0)=
{−1,1/2<x ≤1
∂u
and ∂t
( x , 0 )=x
using h=k=0.1 for 3 time steps.
18. Approximate the solution to the following elliptic partial differential equation
2 2
∂ u ∂ u xy 2 2
+ =e (x + y )
∂ x2 ∂ y2 , 0<x<1, 0<y<1, u(0,y)=1, u(1,y)=ey, 0y1 and u(x,0)=1, u(x,1)=ex,
0x1, using h=k=1/3.
19. Solve 2U = 0 in the region 0 x 4, 0 y 4 under the condition ( h = k = 1 ).
U(0,y) = 0, U(4,y) = 12 + y for 0 y 4
U(x, 0) = 3x, U(x, 4) = x2 for 0 x 4
20. Solve the Poisson equation Uxx + Uyy = -81xy ; 0 < x < 1 and 0 < y < 1 given that U (x, 0) = 0,
U(0, y) = 0, U(1, y) = 100, U(x, 1) = 100 and h = 1/3.
∂u ∂ 2 u
=
21. Solve ∂t ∂ x 2 in 0 < x < 5, t > 0 given that U (x, 0) = 20, U(0, t) = 0, U(5, t) = 100.
Compute U for one time step with h = 1 by Cronk-Nicholson method.