Numerical Solution of Linear
Algebraic Systems Tutorial
Engineering and Numerical Analysis
Department of Electrical Power
Dr. Hayder Yaqoob
2019-2020
1
Difference
• For systems containing only a few equations Cramer’s rule is a useful method of
solution.
• The Gaussian elimination method is the most widely used.
• A method called partial pivoting helps stop rounding error from growing.
• In LU decomposition, the whole point of Gaussian elimination is to replace
the coefficient matrix with one that is triangular.
• There are occasions when direct are not the best way to solve a system of
equations.
2
Cramer’s Rule
Use Cramer’s rule to solve the system
x1 = 1, x2 = 2, x3 = 4
Mike Renfro
3
Gaussian Elimination
Use Gaussian elimination to solve the system
Subtract the first row to the second one and we add the third row with the first.
Multiply the first row by 1/2 and the second one by 1/3
4
Gaussian Elimination
Add the first row with the second one multiplied by 1/2 and the third row with the
second one multiplied by -2
https://www.matesfacil.com/english/high/solving-systems-by-Gaussian-Elimination.html
5
LU Decomposition
Find an LU decomposition of
6
LU Decomposition
7
Jacobi Iteration
Carry out two iterations of the Jacobi method to approximate the solution of
Solution: First Iteration
Second Iteration
8
Gauss-Seidel Iteration
Carry out two iterations of the Gauss-Seidel method to approximate the solution
of
Solution: First Iteration
Second Iteration
9
Exercises
Solve the following using Cramer’s rule
Carry out row operations (with partial pivoting if necessary) to reduce these matrices
to upper triangular form.
10
Exercises
11