0% found this document useful (0 votes)
40 views29 pages

Newton Secany

The document contains information about solving systems of equations using various methods including graphical method, LU decomposition, and Gaussian elimination with pivoting. It includes examples of applying each method to solve systems of 4 equations with 4 unknowns and displays the steps and results.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as XLSX, PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
40 views29 pages

Newton Secany

The document contains information about solving systems of equations using various methods including graphical method, LU decomposition, and Gaussian elimination with pivoting. It includes examples of applying each method to solve systems of 4 equations with 4 unknowns and displays the steps and results.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as XLSX, PDF, TXT or read online on Scribd
You are on page 1/ 29

Chart Title

1.2

0.8

0.6

0.4

0.2

0
1
1. f(x) = 〖 4𝑥 〗 ^3−15𝑥2+17x−6
2. f(x) = e-xx
3. f(x) = 3ex-4cos(x)

Newton-Raphson Method

x 0 1 2 3
f(x) -6 0 0 18
f'(x)=12x2−30𝑥+17

n x0 f(x0) f'(x0) x1
0 1.50000 -0.75000 -1.00000 0.75000
1 0.75000 0.00000 23.75000 0.75000
2 0.75000 0.00000 23.75000 0.75000
3 0.75000 0.00000 23.75000 0.75000
4 0.75000 0.00000 23.75000 0.75000
5 0.75000 0.00000 23.75000 0.75000

f(x) = e-xx x -1 0 1 2 3
f'(x)= -xe +e
-x -x
f(x) -2.718282 0 0.367879 0.270671 0.149361

n x0 f(x0) f'(x0) x1
0 -0.50000 -0.82436 2.47308 -0.16667
1 -0.16667 -0.19689 1.37825 -0.02381
2 -0.02381 -0.02438 1.04848 -0.00055
3 -0.00055 -0.00055 1.00111 0.00000
4 0.00000 0.00000 1.00000 0.00000
5 0.00000 0.00000 1.00000 0.00000

f(x) = 3ex-4cos(x) x 0 1 2 3
f'(x)= 3ex+4sinx f(x) -1 6 24 64

n x0 f(x0) f'(x0) x1
0 0.50000 1.43583 6.86387 0.29081
1 0.29081 0.18050 5.15947 0.25583
2 0.25583 0.00478 4.88679 0.25485
3 0.25485 0.00000 4.87921 0.25485
4 0.25485 0.00000 4.87921 0.25485
5 0.25485 0.00000 4.87921 0.25485

a. Graphical Method

0.5x1-x2=-9.5
1.02x1-2x2=-18.8

0.5x1-x2=-9.5
x1 0 2 4 6 8 10 12
x2 9.5 10.5 11.5 12.5 13.5 14.5 15.5

1.02x1-2x2=-18.8
x1 0 2 4 6 8 10 12
x2 9.4 10.42 11.44 12.46 13.48 14.5 15.52

20

18

16
6; 14.5
6; 14.5
14

12

10

4
0.5x1-x2=-9.5
2 1.02x1-2x2=-

0
0 2 4 6 8 10 12 14 16 18

LU DECOMPOSITION
Column1
6x1-2x2+2x3+4x4= 16
12x1-8x2+6x3+10x4= 26
3x1-13x2+9x3+3x4= -19
-6x1+ 4x2+x3-18x4= -34

x1 x2 x3 x4 B
6 -2 2 4 16
m21=12/6= 2 12 -8 6 10 26
m31= 3/6=0.5 3 -13 9 3 -19
m41= -6/6= -1 -6 4 1 -18 -34

x1 x2 x3 x4 B
6 -2 2 4 16
0 -4 2 2 -6
m32=-12/- 3 0 -12 8 1 -27
m42= 2/-4 -0.5 0 2 3 -14 -18

x1 x2 x3 x4 B
6 -2 2 4 16
0 -4 2 2 -6
0 0 2 -5 -9
m43= 4/2 2 0 0 4 -13 -21

x1 x2 x3 x4 B
6 -2 2 4 16
0 -4 2 2 -6
0 0 2 -5 -9
0 0 0 -3 -3

A=LU
6 -2 2 4 1 0 0 0
12 -8 6 10 = 2 1 0 0
3 -13 9 3 0.5 3 1 0
-6 4 1 -18 -1 -0.5 2 1

LY=B 1 0 0 0 y1 16 y1= 16
2 1 0 0 y2 26
=
0.5 3 1 0 y3 -19
=

-1 -0.5 2 1 y4 -34

UX=Y 6 -2 2 4 x1 16 -3x4=-3
0 -4 2 2 x2 -6 x4= 1
=
0 0 2 -5 x3 -9
0 0 0 -3 x4 -3

Therefore:
x1= 3
x2 = -1
x3 = -2
x4= 1

2x1+4x2-2x3-2x4=-4
x1+2x2+4x3-3x4=5
-3x1-3x2+8x3-2x4=7
-x1+x2+6x3-3x4=7

AX=B
x1 x2 x3 x4 x B
2 4 -2 -2 x1 -4
m21=1/2=0 0.5 1 2 4 -3 x2 5
m31= -3/2= -1.5 -3 -3 8 -2 x3 7
m41= -1/2= -0.5 -1 1 6 -3 x4 7

x1 x2 x3 x4 x B
2 4 -2 -2 x1 -4
0 0 5 -2 x2 7
m32= 3/3= 1 0 3 5 -5 x3 1
m42= 0/3= 0 0 3 5 -4 x4 5

x1 x2 x3 x4 x B
2 4 -2 -2 x1 -4
0 3 5 -4 x2 5
0 0 0 -1 x3 -4
m43= 0/5 0 0 0 5 -2 x4 7

x1 x2 x3 x4 x B
2 4 -2 -2 x1 -4
0 3 5 -4 x2 5
0 0 5 -2 x3 7
0 0 0 -1 x4 -4

A=LU
2 4 -2 -2 1 0
1 2 4 -3 = -0.5 1
-3 -3 8 -2 0.5 0
-1 1 6 -3 -1.5 1

LY=B
1 0 0 0 Y1 -4
-0.5 1 0 0 Y2 = 7
0.5 0 1 0 Y3 5
-1.5 1 0 1 Y4 7

UX=Y
2 4 -2 -2 x1 -4
0 3 5 -4 x2 = 5
0 0 5 -2 x3 7
0 0 0 -1 x4 -4

Therefore:
x1= 1
x2= 2
Gaussian Elimination with Pivoting x3=3
x4=4
6x1-2x2+2x3+4x4= 16
12x1-8x2+6x3+10x4= 26
3x1-13x2+9x3+3x4= -19
-6x1+ 4x2+x3-18x4= -34

6 -2 2 4 16
m21=12/6= 2 12 -8 6 10 26 m21=1/2=0
m31= 3/6= 0.5 3 -13 9 3 -19 m31= -3/2=
m41= -6/6= -1 -6 4 1 -18 -34 m41= -1/2=

6 -2 2 4 16 Swap second and fourth row of


0 -4 2 2 -6
m32=-12/- 3 0 -12 8 1 -27
m42= 2/-4 -0.5 0 2 3 -14 -18
6 -2 2 4 16
0 -4 2 2 -6
0 0 2 -5 -9 m32= 3/3=
m43= 4/2 2 0 0 4 -13 -21 m42= 0/3=

6 -2 2 4 16 Swap third and fourth row of th


0 -4 2 2 -6
0 0 2 -5 -9
0 0 0 -3 -3

-3x4=-3 2x3-5x4= -9 -4x2+2x3+2x4=-6 6x1-2x2+2x3+4x4= 16


x4= 1 2x3-5(1)=-9 -4x2+2(-2)+2(1)=-6 6x1-2(1)+2(-2)+4(1)= 16
x3= (-9+5)/2 x2= (-6+4-2)/-4 x1= (16+2+4-4)/6
x3= -2 x2=1 x1=3 m43= 0/5

Therefore:
x1= 3
x2 = 1
x3 = -2
x4= 1

Naive Gaussian Elimination Method

6x1-2x2+2x3+4x4= 16
12x1-8x2+6x3+10x4= 26
3x1-13x2+9x3+3x4= -19
-6x1+ 4x2+x3-18x4= -34

12x1-8x2+6x3+10x4= 26
A21/A11= 12/6 2 - 12x1-4x2+8x3+8x4= 32
-4x2+2x3+2x4 = -6

3x1-13x2+9x3+3x4= -19
A31/A11= 3/6 1/2 - 3x1-x2+x3+2x4= 8
-12x2+8x3+x4= -27

-6x1+ 4x2+x3-18x4= -34


A41/A11= -6/7 -1 - -6x1+2x2-2x3-4x4= -16
2x2+3x3-14x4= -18

6x1-2x2+2x3+4x4= 16
-4x2+2x3+2x4 = -6
-12x2+8x3+x4= -27
2x2+3x3-14x4= -18

-12x2+8x3+x4= -27
A32/A22=-12/-4 3 - -12x2+6x3+6x4 = -18
2x3-5x4= -9

2x2+3x3-14x4= -18
A42/A22=2/-4 -1/2 - 2x2-x3-x4 = 3
4x3-13x4= -21

6x1-2x2+2x3+4x4= 16
-4x2+2x3+2x4 = -6
2x3-5x4= -9
4x3-13x4= -21

4x3-13x4= -21
A43/A33= 4/2 2 - 4x3-10x4= -19
-3x4= -3

Therefore: x4= 1 2x3-5x4= -9 -4x2+2x3+2x4 = -6 6x1-2x2+2x3+4x4= 16


x1= 3 2x3-5(1)= -9 -4x2+2(-2)+2(1)= -6 6x1-2(1)+2(-2)+4(1)=16
x2 = 1 x3= (-9+5)/2 x2= (-6-2+4)/-4 x1= (16-4+4+2)/6
x3 = -2 x3= -2 x2= 1 x1= 3
x4= 1
Secant Method

x 0 1 2 3
f(x) -6 0 0 18

n x0 f(x0) x1 f(x1) x2 f(x2)


0 0 -6 3 18 0.75 0

x -1 0 1 2 3
f(x) -2.718281828 0 0.367879441 0.2706706 0.149361205

n x0 f(x0) x1 f(x1) x2 f(x2)


0 -1.00000 -2.71828 1.00000 0.36788 0.76159 0.35560
1 1.00000 0.36788 0.76159 0.35560 -6.14512 -2866.28682
2 0.76159 0.35560 -6.14512 -2866.28682 0.76074 0.35551
3 -6.14512 -2866.28682 0.76074 0.35551 0.75988 0.35541
4 0.76074 0.35551 0.75988 0.35541 -2.41173 -26.89862
5 0.75988 0.35541 -2.41173 -26.89862 0.71852 0.35026
6 -2.41173 -26.89862 0.71852 0.35026 0.67828 0.34422
7 0.71852 0.35026 0.67828 0.34422 -1.61521 -8.12280
8 0.67828 0.34422 -1.61521 -8.12280 0.58504 0.32592
9 -1.61521 -8.12280 0.58504 0.32592 0.50017 0.30332
10 0.58504 0.32592 0.50017 0.30332 -0.63892 -1.21038
11 0.50017 0.30332 -0.63892 -1.21038 0.27192 0.20718
12 -0.63892 -1.21038 0.27192 0.20718 0.13880 0.12081 n
13 0.27192 0.20718 0.13880 0.12081 -0.04741 -0.04971 0
14 0.13880 0.12081 -0.04741 -0.04971 0.00687 0.00683 1
15 -0.04741 -0.04971 0.00687 0.00683 0.00032 0.00032 2
16 0.00687 0.00683 0.00032 0.00032 0.00000 0.00000 3
17 0.00032 0.00032 0.00000 0.00000 0.00000 0.00000 4
18 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 5
19 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
20 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000

b. Cramer's Rule

x1+x2-x3= -3
6x1+2x2+2x3= 2
-3x1+4x2+x3= 1

x1+x2-x3= -3 1 1 -1
14 16 18 6x1+2x2+2x3= 2 6 2 2
16.5 17.5 18.5 -3x1+4x2+x3= 1 -3 4 1

14 16 18 1 1 -1 = 1[2(1)-4(2)] - 1[6(1)-(-
D= = 1(-6) - 1(12) - 1(30)
16.54 17.56 18.58 6 2 2 = - 48
-3 4 1

-3 1 -1 = -3[2(1)-4(2)] - 1[2(1)-
Dx1=
= -3(-6) - 1(0) + (-1)(6)
2 2 2 = 12
1 4 1

1 -3 -1 = 1[2(1)-1(2)] - (-3)[6(1)
Dx2= 6 2 2 = 1(0) + 3(12) + (-1)(12)
= 24
-3 1 1

1 1 -3 = 1[2(1)-4(2)] - 1[6(1)-(-
Dx3= 6 2 2 = 1(-6) - 1(12) + (-3)(30)
0.5x1-x2=-9.5 = -108
1.02x1-2x2=-18.8
-3 4 1

14 16 18 20
c. Elimination of Unknown

x1+x2-x3= -3
6x1+2x2+2x3= 2
-3x1+4x2+x3= 1
x1+x2-x3= -3 Equation 1
6x1+2x2+2x3= 2 Equation 2
-3x1+4x2+x3= 1 Equation 3

6(x1+x2-x3= -3) Equation 1 6x1+2x2+2x3= 2


-
6x1+2x2+2x3= 2 Equation 2 2(-3x1+4x2+x3= 1)

6x1+6x2-6x3= -18 Equation 1 6x1+2x2+2x3= 2


- +
6x1+2x2+2x3= 2 Equation 2 -6x1+8x2+2x3= 2
4x2-8x3 = -20 Equation 4 10x2+4x3 = 4

10(4x2-8x3 = -20) Equation 4


-
4(10x2+4x3= 4) Equation 5

40x2-80x3 = -200) Equation 4


-
40x2+16x3= 16 Equation 5
-96x3= -216

x3 = 9/4 4x2 - 8x3 = -20 ; x3= 9/4


4x2 - 8(9/4) = -20
4x2 - 18 = -20
Therefore: x2 = (-20+18)/4
x1= -1/4 x2 = -1/2
x2 = -1/2
x3 = 9/4
x1+x2-x3= -3 ; x2= -1/2, x3=
x1+(-1/2)-(9/4)= -3
x1= -3 + 1/2 +9/4
6 -2 2 4 x1= -1/4
0 -4 2 2
0 0 2 -5
0 0 0 -3

2y1+y2=-6 0.5y1+3y2+y3=-19 -y1-0.5y2+2y3-y4=-34


2(16)+y2=26 0.5(16)+3(-6)+y3=-19 -(16)-0.5*(-6)+2*(-9)-Y4=-34
y2=26-32 y3= -19-8+18 y4= -34+16-3+18
y2=-6 y3= -9 y4= -3

2x3-5x4= -9 -4x2+2x3+2x4=-6 6x1-2x2+2x3+4x4= 16


2x3-5(1)=-9 -4x2+2(-2)+2(1)=-6 6x1-2(1)+2(-2)+4(1)= 16
x3= (-9+5)/2 x2= (-6+4-2)/-4 x1= (16+2+4-4)/6
x3= -2 x2=1 x1=3

x1 x2 x3 x4 x B
2 4 -2 -2 x1 -4
swap 0 3 5 -4 x2 5
0 3 5 -5 x3 1
0 0 5 -2 x4 7

x1 x2 x3 x4 x B
2 4 -2 -2 x1 -4
swap 0 3 5 -4 x2 5
0 0 5 -2 x3 7
0 0 0 -1 x4 -4
0 0 2 4 -2 -2
0 0 0 3 5 -4
1 0 0 0 5 -2
0 1 0 0 0 -1

y1= -4 -0.5y1+y2=7 0.5y1+y3=5-1.5y1+y2+y4=7


-0.5(-4)+y2=7 0.5(-4)+y3=-1.5(-4)+5+y4=7
y2= 7-2 y3= 5+2 y4= 7-6-5
y2= 5 y3=7 y4= -4

x4=4 5x3-2x4= 7 3x2+5x3-4x4=5 2x1+4x2-2x3-2x4= -4


5x3-2(4)=7 3x2+5(3)-4(4)=5 2x1+4(2)-2(3)-2(4)= -4
x3= (7+8)/5 x2= (5+16-15)/3 x1= (-4+8+6-8)/2
x3= 3 x2= 2 x1= 1

2x1+4x2-2x3-2x4=-4
x1+2x2+4x3-3x4=5
-3x1-3x2+8x3-2x4=7
-x1+x2+6x3-3x4=7

2 4 -2 -2 x1 -4
0.5 1 2 4 -3 x2 5
1.5 -3 -3 8 -2 x3 7
-0.5 -1 1 6 -3 x4 7

ond and fourth row of 2 4 -2 -2 x1 -4


0 0 5 -2 x2 7
0 3 5 -5 x3 1
0 3 5 -4 x4 5
2 4 -2 -2 x1 -4
0 3 5 -4 x2 5
1 0 3 5 -5 x3 1
0 0 0 5 -2 x4 7

d and fourth row of th 2 4 -2 -2 x1 -4


0 3 5 -4 x2 5
0 0 0 -1 x3 -4
0 0 5 -2 x4 7

2 4 -2 -2 x1 -4
0 3 5 -4 x2 5
0 0 5 -2 x3 7
0 0 0 0 -1 x4 -4

2 4 -2 -2 x1 -4
0 3 5 -4 x2 5
0 0 5 -2 x3 7
0 0 0 -1 x4 -4

x4=4 5x3-2x4= 7 3x2+5x3-4x4=5 2x1+4x2-2x3-2x4= -4


5x3-2(4)=7 3x2+5(3)-4(4)=5 2x1+4(2)-2(3)-2(4)= -4
x3= (7+8)/5 x2= (5+16-15)/3 x1= (-4+8+6-8)/2
x3= 3 x2= 2 x1= 1

Therefore:
x1= 1
x2= 2
x3=3
x4=4

Naive Gaussian Elimination Method

2x1+4x2-2x3-2x4=-4
x1+2x2+4x3-3x4=5
-3x1-3x2+8x3-2x4=7
-x1+x2+6x3-3x4=7

x1+2x2+4x3-3x4=5
A21/A11=1/2 1/2 - x1+2x2-x3-x4=-2
5x3-2x4= 7

-3x1-3x2+8x3-2x4=7
A31/A11= -3/2 -1.5 - -3x1-6x2+3x3+3x4=6
3x2+5x3-5x4=1

-x1+x2+6x3-3x4=7
A41/A11= -1/2 -1/2 - -x1-2x2+x3+x4=2
3x2+5x3-4x4=5

2x1+4x2-2x3-2x4=-4 2x1+4x2-2x3-2x4=-4
5x3-2x4= 7 Swap 3x2+5x3-4x4=5
3x2+5x3-5x4=1 3x2+5x3-5x4=1
3x2+5x3-4x4=5 5x3-2x4= 7

3x2+5x3-5x4=1
A32/A22=3/3 1 - 3x2+5x3-4x4=5
-x4=-4

x4= 4 5x3-2(4)= 7
x3= (7+8)/5
x3= 3

Therefore:
x1= 1
x2= 2
x3=3
x4=4

6x1-2x2+2x3+4x4= 16
6x1-2(1)+2(-2)+4(1)=16
x1= (16-4+4+2)/6
Jacobi Method

6x1-2x2+2x3+4x4= 16 x1= (16+2x2-2x3-4x4)/6


12x1-8x2+6x3+10x4= 26 x2= (-26+12x1+6x3+10x4)/8
3x1-13x2+9x3+3x4= -19 x3= (-19-3x1+13x2-3x4)/9
-6x1+ 4x2+x3-18x4= -34 x4= (34-6x1+4x2+x3)/18

Assumption:
x1=0 x3=0 n x1
x2=0 x4=0 0 0
1 2.66667
2 1.02778
3 5.84362
4 -0.76569
5 11.46482
6 -8.36002
7 25.60053
8 -31.24422
9 64.79759
10 -96.88330
11 175.65036
12 -283.55786
13 490.24472
14 -813.71206
15 1383.51757
16 -2319.07337
17 3920.04050
18 -6593.49680
x 0 1 2 19 11122.64750
f(x) -1 6 24 -1 20 -18730.69956

x0 f(x0) x1 f(x1) x2 f(x2)


0 -1 1 5.993636 0.142987 -0.49803411
1.00000 5.99364 0.14299 -0.49803 0.20874 -0.21681
0.14299 -0.49803 0.20874 -0.21681 0.25943 0.02241
0.20874 -0.21681 0.25943 0.02241 0.25468 -0.00084
0.25943 0.02241 0.25468 -0.00084 0.25485 0.00000
0.25468 -0.00084 0.25485 0.00000 0.25485 0.00000

-3
2
1

= 1[2(1)-4(2)] - 1[6(1)-(-3)(2)] + (-1)[6(4)-(-3)(2)]


= 1(-6) - 1(12) - 1(30)
= - 48

= -3[2(1)-4(2)] - 1[2(1)-(1)(2)] + (-1)[2(4)-(1)(2)] 𝑥1=𝐷𝑥1/D


= -3(-6) - 1(0) + (-1)(6) x1= 12/-48
= 12 x1= -1/4

= 1[2(1)-1(2)] - (-3)[6(1)-(3)(2)] + (-1)[6(1)-(-3)(2)] 𝑥2=𝐷𝑥2/D


= 1(0) + 3(12) + (-1)(12) x2= 24/-48
= 24 x2= -1/2

= 1[2(1)-4(2)] - 1[6(1)-(-3)(2)] + (-3)[6(4)-(-3)(2)] x3=𝐷𝑥3/D


= 1(-6) - 1(12) + (-3)(30) x2= -108/-48
= -108 x3= 9/4
x1+2x2+2x3= 2 Equation 2
3x1+4x2+x3= 1) Equation 3

6x1+2x2+2x3= 2 Equation 2
-6x1+8x2+2x3= 2 Equation 3
10x2+4x3 = 4 Equation 5
-96

Equation 4
Equation 5

Equation 4
Equation 5

4x2 - 8x3 = -20 ; x3= 9/4


4x2 - 8(9/4) = -20
4x2 - 18 = -20
x2 = (-20+18)/4

x1+x2-x3= -3 ; x2= -1/2, x3= 9/4


x1+(-1/2)-(9/4)= -3
x1= -3 + 1/2 +9/4
+3x3+3x4=6

3x2+5x3-5x4=1 2x1+4x2-2x3-2x4=-4
3x2+5(3)-5(4)=2 2x1+4(2)-2(3)-2(4)=-4
x2= (2+20-15)/3 x1= (-4+8+6-8)/2
x2= 1 x1= 1
Gauss Seidel Method

6x1-2x2+2x3+4x4= 16
12x1-8x2+6x3+10x4= 26
3x1-13x2+9x3+3x4= -19
-6x1+ 4x2+x3-18x4= -34

Assumption: n x1
x2 x3 x4 x1=0 x3=0 0 0
0 0 0 x2=0 x4=0 1 2.66667
-3.25000 -2.11111 1.88889 2 2.84877
1.52778 -8.32407 0.16049 3 2.98633
-7.75077 -0.30041 1.42335 4 3.02037
7.06932 -15.72900 -1.79807 5 3.02578
-18.44287 8.95471 2.84124 6 3.02589
24.21481 -33.51949 -5.53365 7 3.02529
-47.84671 37.49706 8.19444 8 3.02465
73.51664 -82.48802 -15.19405 9 3.02403
-130.97490 119.55902 24.05799 10 3.02343
213.68813 -220.91561 -42.17367 11 3.02284
-366.97875 352.90185 69.39649 12 3.02227
611.64754 -613.87381 -118.60641 13 3.02171
-1037.25015 1015.43454 198.22575 14 3.02117
1741.47517 -1729.85149 -335.61302 15 3.02064
-2940.72299 2896.46139 564.01787 16 3.02012
4949.39472 -4899.00057 -951.86311 17 3.01962
-8345.93937 8237.32676 1602.61214 18 3.01912
14058.07100 -13898.24108 -2701.81521 19 3.01865
### 23402.42879 4551.61243 20 3.01818
39922.05838 -39452.17390 -7670.95757 21 3.01772

values does not converge


Fixed-Point Iteration Method

f(x) = x^3-x-1

x3=x+1
x= 3√x+1
g(x)= 3√x+1 a= 1
b= 2
X0= (a+b)/2 X0= 1.5

x 0 1 2
f(x) -1 -1 5

n x0 x1= g(X0)
0 1.5 1.357208808
1 1.357208808 1.330860959
2 1.330860959 1.325883774
3 1.325883774 1.324939363
4 1.324939363 1.324760011
5 1.324760011 1.324725945
6 1.324725945 1.324719475
7 1.324719475 1.324718245
8 1.324718245 1.324718012
9 1.324718012 1.324717968
10 1.324717968 1.324717959
Jacobi Method

x1= (16+2x2-2x3-4x4)/6 2x1+4x2-2x3-2x4=-4


x2= (-26+12x1+6x3+10x4)/8 x1+2x2+4x3-3x4=5
x3= (-19-3x1+13x2-3x4)/9 -3x1-3x2+8x3-2x4=7
x4= (34-6x1+4x2+x3)/18 -x1+x2+6x3-3x4=7
Assumption:
x1=0 x3=0
x2 x3 x4 x2=0 x4=0
0 0 0 n x1 x2 x3
0.75000 -1.91667 1.06019 0 0.00000 0.00000 0.00000
0.91088 -2.09838 1.02514 1 -2.00000 2.50000 0.87500
0.93713 -2.09463 0.98533 2 -8.45833 -1.75000 0.47917
0.94124 -2.08678 0.97533 3 2.89583 7.14583 -2.72396
0.94276 -2.08305 0.97407 4 -18.15451 7.79167 4.85590
0.94413 -2.08068 0.97447 5 -19.09201 -7.68142 -4.60221
0.94551 -2.07862 0.97509 6 24.78783 45.29123 -5.15824
0.94687 -2.07665 0.97572 7 -105.47493 -11.17878 25.22109
0.94820 -2.07473 0.97633 8 39.76331 -3.92773 -44.32398
0.94950 -2.07286 0.97692 9 41.07239 190.57766 34.19857
0.95077 -2.07104 0.97750 10 -454.50172 -244.75079 61.35752
0.95200 -2.06926 0.97806 11 664.75800 280.88415 -232.36997
0.95320 -2.06752 0.97861 12 -605.83959 420.30897 403.06548
0.95437 -2.06583 0.97915 13 -1034.58369 -1393.25801 -217.45679
0.95551 -2.06418 0.97967 14 3712.90637 2673.47614 -623.10385
0.95663 -2.06258 0.98018 15 -6528.86114 -1442.95300 2256.56719
0.95771 -2.06101 0.98068 16 3545.45541 -3638.73048 -3387.30975
0.95877 -2.05948 0.98116 17 10094.25497 14313.54743 1517.42279
0.95981 -2.05799 0.98163 18 -36283.35353 -21836.99525 6860.88054
0.96081 -2.05654 0.98209 19 54971.81410 11080.83028 -20684.52003
0.96179 -2.05513 0.98254 20 -24313.30009 41687.95375 29404.33677
Case Study

Rah Rav Rev


0 1 1
0 1
0 0 0
0 0 0
0 0 0
0 0 0
1 0 0

1
Gauss Seidel Method

x1= (-4-4x2+2x3+2x4)/2 2x1+4x2-2x3-2x4=-4


x2= (5-x1-4x3+3x4)/2 x1+2x2+4x3-3x4=5
x3= (7+3x1+3x2+2x4)/8 -3x1-3x2+8x3-2x4=7
x4= (-7-x1+x2+6x3)/3 -x1+x2+6x3-3x4=7
Assumption:
x1=0 x3=0
x2=0 x4=0 n x1 x2 x3
x4 0 0 0 0
0.00000 1 -2.00000 3.50000 1.43750
-2.33333 2 -5.18750 5.78125 1.69141
0.91667 3 -7.16536 9.75846 3.02384
0.86111 4 -9.13745 15.05448 5.43280
-6.36458 5 -10.07994 21.56873 9.33236
16.02720 6 -8.92417 28.61878 14.98071
-7.73423 7 -4.11444 34.80940 22.42121
-5.81535 8 6.28612 37.74008 31.25575
79.54090 9 24.43841 33.76352 40.36643
-105.54498 10 52.34728 17.85534 47.57796
115.89889 11 91.19254 -16.26428 49.30441
190.29869 12 140.28954 -75.56877 40.25943
-595.03123 13 195.62971 -166.48458 13.36261
1145.84715 14 248.01890 -292.70399 -39.96012
-556.80501 15 280.95330 -451.79822 -128.81548
-1595.01778 16 268.56616 -630.47432 -260.89426
6206.10376 17 174.25236 -798.54068 -439.18362
-9171.68146 18 -49.06719 -902.04655 -657.03388
4438.94306 19 -455.66833 -856.68940 -891.44104
18534.88050 20 -1098.95134 -543.47590 -1094.75749
-56001.70132 21 -2016.49552 190.22799 -1185.64778
AB AH HG BC BG CG CD
0 0 0 0 0 0 0
-0.9880316241 0 0 0 0 0 0
0.1542514499 1 0 0 0 0 0
0 0 0 0 0 0 0
0 0.9880316241 -0.9880316241 0.98803162409 0 0
0 0.9880316241 0.1542514499 0.15425144989 0 0
0 0 0 0 0 0 0
-1
-0.9880316241 -1 0.988031624
-0.1542514499 0.15425145
-0.9880316241 1
-0.1542514499

-0.988031624
-0.15425145
x1= (-4-4x2+2x3+2x4)/2
x2= (5-x1-4x3+3x4)/2
x3= (7+3x1+3x2+2x4)/8
x4= (-7-x1+x2+6x3)/3

x4
0
2.37500
4.70573
9.35563
16.59624
26.88094
40.14240
55.48371
70.66282
81.50790
81.32527
60.45655
6.23275
-96.31287
-262.49454
-504.21480
-823.80201
-1204.96493
-1600.72755
-1918.88910
-2006.68983
-1638.05440
DG FG DF BH Rah DE EF HG
0 0 0 0 0 0 0 0 6000
0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0
0 0 0 1
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
1
0 0 2000

-0.988031624
0.15425145 1 -1
1
-1 1
0.988031624 0.9880316241 3000
-0.15425145 0.1542514499

-0.9880316241
-0.1542514499 -1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

You might also like