Newton Secany

Download as xlsx, pdf, or txt
Download as xlsx, pdf, or txt
You are on page 1of 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