Department of Mathematical and Computational Sciences
National Institute of Technology Karnataka, Surathkal
Numerical Methods- MA207
Problem Sheet on System of Equations
Default number of decimal places is 4 in your calculations.
1. Solve the following system of equations by Jacobi iterative method and Gauss-Seidel
mehtod:
i) x + 5y − z = 10; x + y + 8z = 20; 4x + 2y + z = 14.
ii)
4 2 3 x 8
3 5 2 y = −14
2 3 8 z 27
2. The upward velocity of a rocket is given at three different times in the following
table:
Time, t (s) Velocity, v (m/s)
5 106.8
8 177.2
12 279.2
The velocity data is approximated by a polynomial as v(t) = a1 t2 + a2 t + a3 , 5 ≤
t ≤ 12. Find the values of a1 , a2 , and a3 using the Gauss-Seidel method. Assume
an initial guess of the solution as [1 2 5]T and conduct two iterations.
3. Solve the following system by Gauss -Seidel method:
2.38x + 1.95y − 3.27z + 1.58t = 2.16
3.21x − 0.86y + 2.42z − 3.20t = 3.28
1.44x + 2.95y − 2.14z + 1.86t = 1.42
4.17x + 3.62y − 1.68z − 2.26t = 5.21
4. i) Write a program to apply Gauss- seidel method to solve AX = B after
checkimg the given matrix A is diagonaly dominant or not.
ii) Solve the following system by your own program in Matlab wtih tolerence10−4
with initial value X(0) = 0
4x1 − x2 − x4 = 0
−x1 + 4x2 − x3 − x5 = 5
−x2 + 4x3 − x6 = 0
−x1 + 4x4 − x5 = 6
−x2 − x4 + 4x5 − x6 = −2
−x3 − x5 + 4x6 = 6.
5. Use SOR with ω = 1.25 to solve
4x1 + 3x2 = 24
3x1 + 4x2 − x3 = 30
−x2 + 4x3 = −24
with X(0) = [1, 1, 1]T
6. Use the power method to find the dominant eigenvalue and the coressponding eigen-
vector for the following matrices
0 11 −5
i) A = −2 17 −7 with X(0) = [1, 1, 1]T .
4 26 −10
1 3 −1
ii) A = 3 2 4 with X(0) = [0.5, 0.5, 0.5]T .
1 4 10
7. Use the power method to find the smallest eigenvalue and the coressponding eigen-
vector for the following matrices [Use the fact that eigenvalues of A−1 are reciprocals
of those of A, so smallest eigenvalue of A is reciprocal of largest eigenvalue of A−1 .]
3 6
i) A = with X(0) = [1, 1]T .
6 14
1 3 −1
ii) A = 3 2 4 with X(0) = [0.5, 0.5, 0.5]T .
1 4 10