Example 1
Uploaded by
s140917Example 1
Uploaded by
s140917Example 1: Solve the following system of liner equations using Naïve Gauss Method.
3 x 1−0.1 x2−0.2 x 3=7.85
0.1 x 1+7 x 2−0.3 x 3=−19.3
0.3 x 1−0.2 x 2+10 x 3=71.4
Solution:
Step (1): Arrange the Equations in matrix [A] and matrix [B].
[ ][ ] [ ]
3 −0.1 −0.2 x 1 7.85
0.1 7 =
−0.3 x 2 −19.3
0.3 −0.2 10 x3 71.4
Step (2): Forward Elimination (First Elimination).
[ ] []
a11 a12 a13 1 ⇒ b1 1
a21 a22 a23 R2−(a21 /a 11) R1 b2 r 2−(a21 /a 11)r 1
a31 a32 a33 R3−(a 31/a 11) R1 b3 r 3−(a 31/a 11)r 1
⇒
A [ i , : ]−( A [ i , j ] / a [ j , j ] ) A [ j, : ] ⇒ B [ i ] −( A [ i , j ] / a [ j, j ] ) B [ j ]
A [ i , : ]−( A [ i , j ] / a [ j , j ] ) A [ j, : ] B [i]−( A [i , j]/a [ j , j])B [ j]
[ ] [ ]
3 −0.1 −0.2 1 ⇒ 7.85 1
0.1 7 −0.3 R 2−(0.1/3) R1 −19.3 r 2−(0.1/3)r 1
0.3 −0.2 10 R 3−(0.3/3)R1 71.4 r 3−(0.3 /3)r 1
[ ][ ] [ ]
3 −0.1 −0.2 x1 7.85
=
0 7.003333 −0.293333 x2 −19.5617
0 −0.19 10 x3 70.615
Step (2*): Forward Elimination (Second Elimination).
[ ]
3 −0.1 −0.2 1
0 7.003333 −0.293333 1
0 −0.19 10 R 3−(−0.19/7.003333) R2
[ ]
⇒ 7.85 1
−19.5617 1
70.615 r 3−(−0.19/7.003333)r 2
⇒
A [ i , : ]−( A [ i , j ] / a [ j , j ] ) A [ j, : ] ⇒ B [ i ] −( A [ i , j ] / a [ j, j ]) B [ j ]
A [ i , : ]−( A [ i , j ] / a [ j , j ] ) A [ j, : ] B [i]−( A [i , j]/a [ j , j])B [ j]
1
[ ][ ] [ ]
3 −0.1 −0.2 x1 7.85
0 7.003333 −0.293333 x2 = −19.5617
0 0 10.012 x3 70.0843
Step (3): Backward Substitution.
d3
1. x 3=
u33
d 2−u23 x 3
2. x 2=
u22
d 1−u12 x 2−u13 x 3
3. x 1=
u 11
[ ][ ] [ ] [ ][ ] [ ]
u11 u 12 u13 x 1 d1 ⇒ 3 −0.1 −0.2 x1 7.85
0 u 22 u23 x 2 = d2 0 7.003333 −0.293333 x 2 = −19.5617
0 0 u33 x3 d3 0 0 10.012 x3 70.0843
d 3 70.0843
x 3= = =7
u33 10.012
d 2−u23 x 3 −19.5617−(−0.293333)(7)
x 2= = =−2.50
u22 7.003333
d 1−u12 x 2−u13 x 3 7.85−(−0.1)(−2.50)−(−0.2)(7)
x 1= = =3
u 11 3
2
Code Solution: