0% found this document useful (0 votes)
61 views28 pages

Compilations of Algorithms

The third document demonstrates the secant method to find the root of the function f(x)=x^2-x-1 by iteratively improving the

Uploaded by

Jessica Naces
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)
61 views28 pages

Compilations of Algorithms

The third document demonstrates the secant method to find the root of the function f(x)=x^2-x-1 by iteratively improving the

Uploaded by

Jessica Naces
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/ 28

Gauss Elimination With Row Pivoting

4x1+x2-x3+x4=-2
x1+4x2-x3-x4=-1
(-)x1-x2+5x3+x4=0
x1-x2+2x3+3x4=1

x1 x2 x3 x4 B
4 1 -1 1 -2 Pivot
1 4 -1 -1 -1 M21=1/4
-1 -1 5 1 0 M31=-1/4
1 -1 2 3 1 M41=1/4

x1 x2 x3 x4 B
4 1 -1 1 -2
0 3.75 -0.75 -1.25 -0.5 Pivot
0 -0.75 4.75 1.25 -0.5 M32=-1/5
0 -1.25 2.25 2.75 1.5 M42=-1/3

x1 x2 x3 x4 B
4 1 -1 1 -2
0 3.75 -0.75 -1.25 -0.5
0 0 4.6 1 -0.6 Pivot
0 0 2 2.333333 1.333333 M43=1/2.3

x1 x2 x3 x4 B
4 1 -1 1 -2
0 3.75 -0.75 -1.25 -0.5
0 0 4.6 1 -0.6
0 0 0 1.898551 1.594203

x4 = 0.839695
x3 = -0.312977
x2 = 0.083969
x1 = -0.80916
27x1+6x2-x3=85
6x1+15x2+2x3=72
x1+x2+54x3=110

x1 x2 x3 B
27 6 -1 85 Pivot
6 12 2 72 M21=2/9
1 1 54 110 M31=1/27

x1 x2 x3 B
27 6 -1 85
0 10.66667 2.222222 53.11111 Pivot
0 0.777778 54.03704 106.8519 M32=0.77778/10.66667

x1 x2 x3 B
27 6 -1 85
0 10.66667 2.222222 53.11111
0 0 54.03704 106.8519

x3 = 1.977382
x2 = 4.567212
x1 = 2.206448
Gauss Seidel Method

9x1+2x2+3x3=7 9 2 3 7 x1=(7-2x2-3x3)/9
x1+12x2+9x3=2 1 12 9 2 x2=(2-x1-9x3)/12
4x1+6x2+14x3=1 4 6 14 1 x3=(1-4x1-6x2)/14

i x1 x2 x3
0 0 0 0
1 0.777778 0.101852 -0.194444
2 0.819959 0.24417 -0.26749
3 0.812681 0.299561 -0.289149
4 0.807592 0.316229 -0.294839
5 0.805784 0.320647 -0.296216
6 0.805261 0.321723 -0.296528
7 0.805126 0.321968 -0.296594
8 0.805094 0.322021 -0.296607
9 0.805087 0.322032 -0.29661
10 0.805085 0.322034 -0.29661
11 0.805085 0.322034 -0.29661
12 0.805085 0.322034 -0.29661
13 0.805085 0.322034 -0.29661
14 0.805085 0.322034 -0.29661
15 0.805085 0.322034 -0.29661

27x1+6x2-x3=85 27 6 -1 85 x1=(85-6x2+x3)/27
6x1+15x2+2x3=72 6 15 2 72 x2=(72-6x1-2x3)/15
x1+x2+54x3=110 1 1 54 110 x3=(110-x1-x2)/54

i x1 x2 x3
0 0 0 0
1 3.148148 3.540741 7.379365
2 2.634627 2.762234 7.471653
3 2.811046 2.679361 7.464971
4 2.829215 2.672985 7.464129
5 2.830601 2.672543 7.464061
6 2.830697 2.672513 7.464056
7 2.830703 2.672511 7.464056
8 2.830703 2.672511 7.464056
9 2.830703 2.672511 7.464056
10 2.830703 2.672511 7.464056
11 2.830703 2.672511 7.464056
12 2.830703 2.672511 7.464056
13 2.830703 2.672511 7.464056
14 2.830703 2.672511 7.464056
15 2.830703 2.672511 7.464056
x1=(7-2x2-3x3)/9
x2=(2-x1-9x3)/12
x3=(1-4x1-6x2)/14

x1=(85-6x2+x3)/27
x2=(72-6x1-2x3)/15
x3=(110-x1-x2)/54
Jacobi Method

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

i x1 x2 x3
0 1 1 1
1 1.25 1.66666666666667 -0.8571428571429
2 0.95238095238095 2.05952380952381 -0.9166666666667
3 1.03571428571429 1.9702380952381 -1.015306122449
4 0.98873299319728 2.01445578231292 -0.9906462585034
5 1.00595238095238 1.99468537414966 -1.0036746841594
6 0.99775267249757 2.00259657434402 -0.9983904275996
7 1.0010515366861 1.99898262876579 -1.0006919859781
8 0.99957266069693 2.00046584322504 -0.9997044417257
9 1.00019035037484 1.99980829385326 -1.000127597504
10 0.99992017408731 2.00008471637562 -0.9999454204969
11 1.00003482396968 1.99996429477859 -1.0000235060412
12 0.99998519718435 2.00001552566342 -0.9999899244013
13 1.00000640031554 1.99999338646166 -1.0000043326399
14 0.99999726345545 2.00000285554516 -0.999998140878
15 1.00000117866678 1.99999877796482 -1.00000079887
16 0.99999949477372 2.00000052603392 -0.9999996570426
17 1.00000021724784 1.99999977443167 -1.0000001473229
18 0.9999999067772 2.00000009696976 -0.9999999367405
19 1.0000000400573 1.99999995838249 -1.0000000271704
20 0.99999998280303 2.00000001788083 -0.9999999883322
21 1.00000000738716 1.99999999232304 -1.0000000050111
22 0.99999999682798 2.00000000329757 -0.999999997848
23 1.0000000013624 1.99999999858399 -1.0000000009242
24 0.99999999941494 2.00000000060817 -0.9999999996031
25 1.00000000025127 1.99999999973883 -1.0000000001705
26 0.99999999989209 2.00000000011217 -0.9999999999268

5x1-x2+2x3=12 5 -1 2
3x1+8x2-2x3=-25 3 8 -2
x1+x2+4x3=6 1 1 4

i x1 x2 x3
0 2 2 2
1 2 -3.375 0.5
2 1.525 -3.75 1.84375
3 0.9125 -3.2359375 2.05625
4 0.9303125 -2.953125 2.080859375
5 0.97703125 -2.95365234375 2.005703125
6 1.00698828125 -2.9899609375 1.9941552734375
7 1.004345703125 -3.0040817871094 1.9957431640625
8 1.00088637695313 -3.0026938476563 1.99993402099609
9 0.99948762207031 -3.0003488861084 2.00045186767578
10 0.99974947570801 -2.9996948913574 2.00021531600952
11 0.99997489532471 -2.9998522243881 1.99998635391235
12 1.00003501355743 -2.9999939972687 1.99996933226585
13 1.00001346763992 -3.0000207970176 1.99998974592781
14 0.99999994222536 -3.000007613883 2.00000183234441
15 0.99999774428563 -2.9999995202484 2.00000191791441
16 0.99999932878455 -2.9999986746285 2.00000044399069
17 1.00000008747802 -2.9999996372965 1.99999983646099
18 1.0000001379563 -3.000000073689 1.99999988745463
19 1.00000003028035 -3.00000007987 1.99999998393318
20 0.99999999045274 -3.0000000153718 2.0000000123974
21 0.99999999196667 -2.9999999933204 2.00000000622977
22 0.99999999884401 -2.9999999954301 2.00000000033844
23 1.00000000077861 -2.9999999994819 1.99999999914651
24 1.00000000044502 -3.0000000005054 1.99999999967582
25 1.0000000000286 -3.0000000002479 2.00000000001508
26 0.99999999994438 -3.000000000007 2.00000000005483
27 0.99999999997668 -2.9999999999654 2.00000000001564
3 x1=(3+x2+x3)/4
9 x2=(9+2x1-x3)/6
-6 x3=(-6+x1-x2)/7

12 x1=(12+x2-2x3)/5
25 x2=(-25-3x1+2x3)/8
6 x3=(6-x1-x2)/4
Secant Method

f(x)=x2-x-1 xa=3 xb=2 xnew=xa-((f(xa)(xa-xb))/(f(xa)-f(xb)))

i xa xb f(xa) f(xb) rel. error


1 3 2 5 1
2 2 1.75 1 0.3125 14.29%
3 1.75 1.636364 0.3125 0.041322 6.94%
4 1.636364 1.619048 0.041322 0.002268 1.07%
5 1.619048 1.618042 0.002268 1.842E-05 0.06%
6 1.618042 1.618034 1.842E-05 8.346E-09 0.00%
7 1.618034 1.618034 8.346E-09 3.109E-14 0.00%
8 1.618034 1.618034 3.109E-14 0 0.00%
9 1.618034 1.618034 0 0 0.00%

f(x)=x3-4x2-1=0 xa=3 xb=2 xnew=xa-((f(xa)(xa-xb))/(f(xa)-f(xb)))

i xa xb f(xa) f(xb) rel. error


1 1 2 -4 -9
2 2 0.2 -9 -1.152 900.00%
3 0.2 -0.06422 -1.152 -1.016762 -411.43%
4 -0.06422 -2.050707 -1.016762 -26.44565 -96.87%
5 -2.050707 0.015209 -26.44565 -1.000922 ###
6 0.015209 0.096476 -1.000922 -1.036332 84.24%
7 0.096476 -2.281902 -1.036332 -33.71034 -104.23%
8 -2.281902 0.171911 -33.71034 -1.113134 1427.37%
9 0.171911 0.255705 -1.113134 -1.24482 32.77%
10 0.255705 -0.536384 -1.24482 -2.305152 -147.67%
11 -0.536384 1.185609 -2.305152 -4.956103 145.24%
12 1.185609 -2.033754 -4.956103 -25.95655 -158.30%
13 -2.033754 1.945379 -25.95655 -8.77571 204.54%
14 1.945379 3.977858 -8.77571 -1.350357 51.09%
15 3.977858 4.34748 -1.350357 5.56758 8.50%
16 4.34748 4.050007 5.56758 -0.179755 7.35%
17 4.050007 4.059311 -0.179755 -0.022673 0.23%
18 4.059311 4.060654 -0.022673 0.000117 0.03%
19 4.060654 4.060647 0.000117 -7.52E-08 0.00%
20 4.060647 4.060647 -7.52E-08 -2.42E-13 0.00%
21 4.060647 4.060647 -2.42E-13 0 0.00%
22 4.060647 4.060647 0 -2.42E-13 0.00%
Fixed Point Iteration

f(x)=x2-2x-3 x=g(x)
xnew=g(xold)
x2-2x-3=0 x2-2x-3=0
x2=2x+3 2x=x2-3 x2-x-1=0
x=√(2𝑥+3) x=(x2-3)/2 1.618 or x=-0.618
x -x=1
2

x(x-1)=1
x=1/(x-1)
i x g(x)
1 0 1.732051 i x g(x)
2 1.732051 2.54246 1 0 -1
3 2.54246 2.843399 2 -1 -0.5
4 2.843399 2.947338 3 -0.5 -0.666667
5 2.947338 2.982394 4 -0.666667 -0.6
6 2.982394 2.994126 5 -0.6 -0.625
7 2.994126 2.998041 6 -0.625 -0.615385
8 2.998041 2.999347 7 -0.615385 -0.619048
9 2.999347 2.999782 8 -0.619048 -0.617647
10 2.999782 2.999927 9 -0.617647 -0.618182
11 2.999927 2.999976 10 -0.618182 -0.617978
12 2.999976 2.999992 11 -0.617978 -0.618056
13 2.999992 2.999997 12 -0.618056 -0.618026
14 2.999997 2.999999 13 -0.618026 -0.618037
15 2.999999 3 14 -0.618037 -0.618033
16 3 3 15 -0.618033 -0.618034
17 3 3 16 -0.618034 -0.618034
18 3 3 17 -0.618034 -0.618034
18 -0.618034 -0.618034
x=3 19 -0.618034 -0.618034

x=-0.618
Newton Raphson Method xnew=xold-f(xold)/f'(xold)

x2-x-1=0
f(x)=x2-x-1=0 xnew=x-((x2-x-1)/(2x-1))
f'(x)=2x-1=0

i x f(x) f'(x) rel. error


1 1 -1 1
2 2 1 3 50.0%
3 1.666667 0.111111 2.333333 20.0%
4 1.619048 0.002268 2.238095 2.9%
5 1.618034 1.027E-06 2.236069 0.1%
6 1.618034 2.105E-13 2.236068 0.0%
7 1.618034 0 2.236068 0.0%
8 1.618034 0 2.236068 0.0%
9 1.618034 0 2.236068 0.0%
10 1.618034 0 2.236068 0.0%

x3-4x2-1=0
f(x)=x3-4x2-1=0 xnew=x-((x3-4x2-1)/(3x2-8x))
f'(x)=3x2-8x=0

i x f(x) f'(x) rel. error


1 1 -4 -5
2 0.2 -1.16 -0.6 400.0%
3 -1.733333 3.737778 -4.466667 111.5%
4 -0.896517 0.700261 -2.793035 93.3%
5 -0.645801 0.062859 -2.291601 38.8%
6 -0.61837 0.000752 -2.236741 4.4%
7 -0.618034 1.132E-07 -2.236068 0.1%
8 -0.618034 0 -2.236068 0.0%
9 -0.618034 0 -2.236068 0.0%
10 -0.618034 0 -2.236068 0.0%
Bisection Method

f(x)= x3-4x2-1=0 x= -0.618034

i xL xR f(xL) f(xR) xm f(xm)


1 -1 0 -6 -1 -0.5 -2.125
2 -1 -0.5 1 -0.25 -0.75 0.3125
3 -0.75 -0.5 0.3125 -0.25 -0.625 0.015625
4 -0.625 -0.5 0.015625 -0.25 -0.5625 -0.121094
5 -0.625 -0.5625 0.015625 -0.121094 -0.59375 -0.053711
6 -0.625 -0.59375 0.015625 -0.053711 -0.609375 -0.019287
7 -0.625 -0.609375 0.015625 -0.019287 -0.617188 -0.001892
8 -0.625 -0.617188 0.015625 -0.001892 -0.621094 0.006851
9 -0.621094 -0.617188 0.006851 -0.001892 -0.619141 0.002476
10 -0.619141 -0.617188 0.002476 -0.001892 -0.618164 0.000291
11 -0.618164 -0.617188 0.000291 -0.001892 -0.617676 -0.000801
12 -0.618164 -0.617676 0.000291 -0.000801 -0.61792 -0.000255
13 -0.618164 -0.61792 0.000291 -0.000255 -0.618042 1.79E-05
14 -0.618042 -0.61792 1.79E-05 -0.000255 -0.617981 -0.000119
15 -0.618042 -0.617981 1.79E-05 -0.000119 -0.618011 -5.03E-05
16 -0.618042 -0.618011 1.79E-05 -5.03E-05 -0.618027 -1.62E-05

f(x)= x2-x-1=0 x= 1.618034

i xL xR f(xL) f(xR) xm f(xm)


1 1 2 -1 1 1.5 -0.25
2 1.5 2 -0.25 1 1.75 0.3125
3 1.5 1.75 -0.25 0.3125 1.625 0.015625
4 1.5 1.625 -0.25 0.015625 1.5625 -0.121094
5 1.5625 1.625 -0.121094 0.015625 1.59375 -0.053711
6 1.59375 1.625 -0.053711 0.015625 1.609375 -0.019287
7 1.609375 1.625 -0.019287 0.015625 1.617188 -0.001892
8 1.617188 1.625 -0.001892 0.015625 1.621094 0.006851
9 1.617188 1.621094 -0.001892 0.006851 1.619141 0.002476
10 1.617188 1.619141 -0.001892 0.002476 1.618164 0.000291
11 1.617188 1.618164 -0.001892 0.000291 1.617676 -0.000801
12 1.617676 1.618164 -0.000801 0.000291 1.61792 -0.000255
13 1.61792 1.618164 -0.000255 0.000291 1.618042 1.79E-05
14 1.61792 1.618042 -0.000255 1.79E-05 1.617981 -0.000119
15 1.617981 1.618042 -0.000119 1.79E-05 1.618011 -5.03E-05
16 1.618011 1.618042 -5.03E-05 1.79E-05 1.618027 -1.62E-05
1.618027 1.618042 -1.62E-05 1.79E-05 1.618034 8.364E-07
False Position Method xn=xL+(xR-xL)(yL/(yL-yR))

f(x)= 4x2+3x-3=0 x= 0.568729

i xL xR yL yR xn yn rel. error
1 0 1 -3 4 0.428571 -0.979592
2 0.428571 1 -0.979592 4 0.540984 -0.206396 0.207792
3 0.540984 1 -0.206396 4 0.563506 -0.039324 0.039969
4 0.563506 1 -0.039324 4 0.567756 -0.007347 0.007485
5 0.567756 1 -0.007347 4 0.568548 -0.001368 0.001394
6 0.568548 1 -0.001368 4 0.568696 -0.000254 0.000259
7 0.568696 1 -0.000254 4 0.568723 -4.73E-05 4.823E-05
8 0.568723 1 -4.73E-05 4 0.568728 -8.8E-06 8.971E-06
9 0.568728 1 -8.8E-06 4 0.568729 -1.64E-06 1.669E-06
10 0.568729 1 -1.64E-06 4 0.568729 -3.04E-07 3.103E-07
11 0.568729 1 -3.04E-07 4 0.568729 -5.66E-08 5.772E-08
12 0.568729 1 -5.66E-08 4 0.568729 -1.05E-08 1.074E-08
13 0.568729 1 -1.05E-08 4 0.568729 -1.96E-09 2E-09
14 0.568729 1 -1.96E-09 4 0.568729 -3.64E-10 3.714E-10
15 0.568729 1 -3.64E-10 4 0.568729 -6.78E-11 6.908E-11
16 0.568729 1 -6.78E-11 4 0.568729 -1.26E-11 1.285E-11
17 0.568729 1 -1.26E-11 4 0.568729 -2.34E-12 2.39E-12
18 0.568729 1 -2.34E-12 4 0.568729 -4.36E-13 4.445E-13
19 0.568729 1 -4.36E-13 4 0.568729 -8.13E-14 8.257E-14
20 0.568729 1 -8.13E-14 4 0.568729 -1.51E-14 1.542E-14
21 0.568729 1 -1.51E-14 4 0.568729 0 0
22 0.568729 1 0 4 0.568729 0 0

f(x)=x2-x-1=0 x=1.61803398874989 xn=xL+(xR-xL)(yL/(yL-yR))

i xL xR yL yR xn yn rel. error
1 1 2 -1 1 1.5 0.25
2 1.5 2 -0.25 1 1.6 0.56 0.005882
3 1.5 1.6 -0.25 -0.04 1.619048 1.661315 0.00112
4 1.5 1.619048 -0.25 0.002268 1.617978 1.615584 6.295E-05
5 1.5 1.617978 -0.25 -0.000126 1.618037 1.61817 3.506E-06
6 1.5 1.618037 -0.25 7.036E-06 1.618034 1.618026 1.954E-07
7 1.5 1.618034 -0.25 -3.92E-07 1.618034 1.618034 1.089E-08
8 1.5 1.618034 -0.25 2.185E-08 1.618034 1.618034 6.069E-10
9 1.5 1.618034 -0.25 -1.22E-09 1.618034 1.618034 3.382E-11
10 1.5 1.618034 -0.25 6.786E-11 1.618034 1.618034 1.885E-12
11 1.5 1.618034 -0.25 -3.78E-12 1.618034 1.618034 1.05E-13
12 1.5 1.618034 -0.25 2.105E-13 1.618034 1.618034 5.852E-15
13 1.5 1.618034 -0.25 -1.2E-14 1.618034 1.618034 3.396E-16
14 1.5 1.618034 -0.25 0 1.618034 1.618034 0
15 1.5 1.618034 -0.25 0 1.618034 1.618034 0
xL)(yL/(yL-yR))
Least Square Line

Write a linear equation that "best fits" the data in the table shown below:

x 1 2 3 4 5 6 7
y 1.5 3.8 6.7 9 11.2 13.6 16

x y xy x2
1 1.5 1.5 1 y=mx+b m=
2 3.8 7.6 4 y=b0+b1x m=
3 6.7 20.1 9
4 9 36 16 slope=m=b1 b=
5 11.2 56 25 y-int=b=b0 b=
6 13.6 81.6 36
7 16 112 49 y=
∑ 28 61.8 314.8 140 y=

∑x= 28
∑y= 61.8 x y
∑xy= 314.8 1 1.5
∑x2= 140 2 3.8
n= 7 3 6.7
4 9
5 11.2
6 13.6
7 16
slope: 2.41428571
y intercept: -0.8285714

Find the least squares line given

x 0 2 4 6 8 10 12
y -0.8 -1 -0.2 0.2 -2 0.8 -0.6

x y xy x2
0 -0.8 0 0 y=mx+b m=
2 -1 -2 4 y=b0+b1x m=
4 -0.2 -0.8 16
6 0.2 1.2 36 slope=m=b1 b=
8 -2 -16 64 y-int=b=b0 b=
10 0.8 8 100
12 -0.6 -7.2 144 y=
∑ 42 -3.6 -16.8 364 y=

∑x= 42
∑y= -3.6 x y
∑xy= -16.8 0 -0.8
∑x2= 364 2 -1
n= 7 4 -0.2
6 0.2
8 -2
10 0.8
12 -0.6
slope: 0.04285714
y intercept: -0.7714286
(n∑xy-∑x∑y)/((n∑x2-(∑x)2)
2.41428571

(∑y-m∑x)/n
-0.8285714

mx+b
2.414x+(-0.828)

(n∑xy-∑x∑y)/((n∑x2-(∑x)2)
0.04285714
(∑y-m∑x)/n
-0.7714286

mx+b
0.0428571x+(-0.771429)
SECANT METHOD
f(x)=0.471405*((20x)^(5/3)/(20+2x)^(2/3))-5

i xa xb f(xa) f(xb)
0.0000 2.0000 3.0000 21.5066 44.3928
1.0000 3.0000 1.0603 44.3928 4.7189
2.0000 1.0603 0.8296 4.7189 1.5480
3.0000 0.8296 0.7169 1.5480 0.1703
4.0000 0.7169 0.7030 0.1703 0.0084
5.0000 0.7030 0.7023 0.0084 0.0001
6.0000 0.7023 0.7023 0.0001 0.0000
7.0000 0.7023 0.7023 0.0000 0.0000
8.0000 0.7023 0.7023 0.0000 0.0000

FIXED POINT ITERATION METHOD


f(x)=0.471405*((20x)^(5/3)/(20+2x)^(2/3))-5
f'(x)=(69.466098*((5x^(2/3)*(20+2x)^(2/3))/3-(4x^(5/3))/(3*(20+2x)^(1/3))))/((20+2x)^(4/3))

i x g(x) x-g(x) g'(x)


0 2 21.50656 -19.50656 20.61599
1 21.50656 724.6396 -703.133 41.10463
2 724.6396 31417.79 -30693.15 43.75638
3 31417.79 1374589 -1343171 43.7609
4 1374589 60153576-58778988 43.7609
5 60153576 2.63E+09 -2.57E+09 43.7609
6 2.63E+09 1.15E+11 -1.13E+11 43.7609
7 1.15E+11 5.04E+12 -4.93E+12 43.7609

The equation don't have solution therefore, it DIVERGE.


NEWTON RAPHSON METHOD
F(X)=0.471405*((20x)^(5/3)/(20+2x)^(2/3))-5
f'(x)=(69.466098*((5x^(2/3)*(20+2x)^(2/3))/3-(4x^(5/3))/(3*(20+2x)^(1/3))))/((20+2x)^(4/3))

i x f(x) f'(x)
0.0000 2.0000 21.5066 20.6160
1.0000 0.9568 3.2415 13.8543
2.0000 0.7228 0.2394 11.7548
3.0000 0.7025 0.0020 11.5560
4.0000 0.7023 0.0000 11.5543
5.0000 0.7023 0.0000 11.5543
BISECTION METHOD

f(x)=0.471405*((20x)^(5/3)/(20+2x)^(2/3))-5
Xm=(Xl+xR)/2

i XL XR F(XL) F(XR) Xm F(Xm)


1.0000 0.6000 0.8000 -1.1292 1.1748 0.7000 -0.0265
2.0000 0.7000 0.8000 -0.0265 1.1748 0.7500 0.5623
3.0000 0.7000 0.7500 -0.0265 0.5623 0.7250 0.2649
4.0000 0.7000 0.7250 -0.0265 0.2649 0.7125 0.1184
5.0000 0.7000 0.7125 -0.0265 0.1184 0.7063 0.0458
6.0000 0.7000 0.7063 -0.0265 0.0458 0.7031 0.0096
7.0000 0.7000 0.7031 -0.0265 0.0096 0.7016 -0.0084
8.0000 0.7016 0.7031 -0.0084 0.0096 0.7023 0.0006
FALSE POSITION METHOD

f(x)=0.471405*((20x)^(5/3)/(20+2x)^(2/3))-5
Xn=xL-f(xl)*(XR - XL)/(f(XR) - f(XL))

i XL XR F(XL) F(XR) Xn F(Xn)


1.0000 0.6000 0.8000 -1.1292 1.1748 0.6980 -0.0493
2.0000 0.6980 0.8000 -0.0493 1.1748 0.7021 -0.0019
3.0000 0.7021 0.8000 -0.0019 1.1748 0.7023 -0.0001

Answer: The Depth is 0.7023m


12x1-8x2+6x3+10x4=26
3x1-13x2+9x3+3x4=-19
6x1-2x2+2x3+4x4=16
-6x1+4x2+x3-18x4=-34
JACOBI METHOD
x1 x2 x3 x4
k
0 2 2 2 2
1 0.83333 3.76923 0.00000 1.77778
2 3.19801 2.06410 5.71368 2.44872
3 -1.35470 6.72025 -4.42735 1.59900
4 7.52801 -1.54717 15.58635 3.58788
5 -9.64786 14.81729 -23.30695 -0.09835
6 23.78030 -16.92317 51.95757 7.10274
7 -41.01318 44.55902 -94.46955 -6.91205
8 84.86750 -75.00013 189.42266 20.21365
9 -159.38945 156.84980 -362.02990 -32.54349
10 314.86772 -293.46676 708.10514 69.76145
11 -605.66495 580.44875 -1369.59283 -128.94268
12 1181.38115 -1116.24295 2663.32896 256.67733
13 -2287.55755 2177.16431 -5165.74105 -491.99610
14 4446.47681 -4216.25618 10031.82917 961.23675
15 -8625.61600 8194.50794 -19470.16013 -1859.89203
16 16750.16205 -15897.61271 37799.13998 3616.42042
17 -32509.49548 30870.07748 -73372.93969 -7014.34906
18 63113.97905 -59916.07622 142435.26206 13622.13017
19 -122511.29032 116318.82202 -276494.27373 -26437.71762
20 237826.61622 -225790.49903 536736.12823 51326.81979
21 -461665.24663 438315.72785 -1041915.98729 -99630.64251
22 896196.18098 -850854.81177 2022580.75276 193409.57799
23 -1739699.39922 1651697.15783 -3926254.51069 -375443.42112
24 3377130.37816 -3206284.61979 7621690.19774 728820.25089
25 -6555716.22113 6224083.28206 -14795308.25605 -1414788.69724
26 12726035.73044 -12082251.54304 28720817.33995 2746408.67772
27 -24703914.76344 23454208.32954 -55753168.08980 -5331353.84527
28 47955519.96899 -45529561.51033 108228668.31041 10349288.21161
29 -93091779.83844 88382496.79504 -210094689.84052 -20090147.97471
30 180710801.59588 -171569074.84648 407838140.25979 38999222.79835
31 -350797803.52617 333051796.50160 -791699917.23082 -75705718.59452
32 680972590.61191 -646523830.95688 1536856652.26914 146960784.66292
33 ### 1255039231.94206 ### ###
34 2566109321.31827 ### 5791337633.33125 553792388.25552
CONCLUSION: THE EQUATION DIVERGE
12x1=26+8x2-6x3-10x4 x1=(13/6)+(2x2/3)-(1x3/2)-(5x4/6)
-13x2=-19-3x1-9x3-3x4 x2=(19/13)+(3x1/13)+(9x3/13)+(3x4/13)
2x3=16-6x1+2x2-4x4 x3=8-(3x1)+(x2)-(2x4)
-18x4=-34+6x1-4x2-x3 x4=(17/9)-(1x1/3)+(2x2/9)+(1x3/18)
GAUSS-SIEDEL METHOD
x1 x2 x3 x4
k
0 2 2 2 2
1 0.83333333 3.5 5 2.66666667
2 -0.22222222 5.48717949 8.820512821 3.67236467
3 -1.64577398 8.03572211 13.62831471 4.98032484
4 -3.44061331 11.2518436 19.61303386 6.62578267
5 -5.66010675 15.262641 26.99139587 8.66681112
6 -8.3762799 20.2149351 36.01015256 11.173754
7 -11.6732479 26.2763763 46.94861188 14.2274225
8 -15.6462405 33.6370041 60.12088059 17.9192411
9 -20.4018052 42.5107872 75.87772046 22.3517611
10 -26.0581364 53.1369506 94.60783756 27.6391366
11 -32.7452322 65.7809424 116.7383658 33.9074182
12 -40.6047365 80.7348721 142.734245 41.294564
13 -49.7893445 98.3172203 173.0961257 49.9500597
14 -60.4616324 118.87157 208.356348 60.0340235
15 -72.7921468 142.764059 249.0724519 71.7156426
16 -86.9565557 170.379179 295.817561 85.1707562
17 -103.131625 202.113496 349.1668571 100.578366
18 -121.489736 238.366739 409.6792159 118.11581
19 -142.191624 279.529654 477.8729058 137.952293
20 -165.376927 325.967865 554.1940623 160.240394
21 -191.152116 378.000876 638.9764364 185.105147
22 -219.575256 435.8752 732.3906747 212.630167
23 -250.63701 499.730426 834.3811223 242.840272
24 -284.23717 569.556877 944.5878437 275.67991
25 -320.155928 645.143272 1062.251238 310.986661
26 -358.018988 726.012628 1186.09627 348.458928
27 -397.25549 811.344365 1314.19298 387.616855
28 -437.047625 899.880347 1443.789513 427.75537
29 -476.270667 989.812287 1571.113548 467.88815
30 -513.422041 1078.64771 1691.137536 506.681146
31 -546.537915 1163.05135 1797.3028 542.374204
32 -573.095672 1238.65852 1881.19713 572.689181
33 -589.900534 1299.85693 1932.180172 594.722839
34 -592.954497 1339.53281 1936.950627 604.822715
35 -577.3057 1348.77744 1879.049108 598.444059
36 -536.876312 1316.54963 1738.290449 569.985936
37 -464.267083 1229.29005 1490.119424 512.604562
38 -350.536816 1070.48293 1104.884252 418.002048
39 -184.955214 820.161443 547.0229903 276.188892
40 45.2720576 454.353058 -225.8408974 75.2188328
41 355.306794 -55.5377844 -1263.895831 -199.104874
42 763.010121 -743.411287 -2626.231902 -563.552099
43 1289.30184 -1649.2183 -4382.019627 -1037.81688
44 1958.54502 -2819.76863 -6611.769916 -1644.89525
45 2798.95191 -4309.59687 -9408.6621 -2411.48672
46 3843.00541 -6181.87714 -12879.91992 -3368.4145
47 5127.88728 -8509.37392 -17148.20676 -4551.05701
48 6695.90161 -11375.4098 -22353.00059 -5999.78052
49 8594.87754 -14874.8242 -28651.89576 -7760.35876
50 10878.5307 -19114.8881 -36221.76283 -9884.3611
51 13606.7569 -24216.1291 -45259.67758 -12429.4853
52 16845.8238 -30313.0064 -55983.50721 -15459.8042
53 20668.4195 -37554.363 -68632.01319 -19045.8879
54 25153.5045 -46103.5591 -83464.2969 -23264.7534
55 30385.9035 -56138.1709 -100758.3747 -28199.5822
56 36455.5585 -67849.111 -120808.6222 -33939.1343
57 43456.3491 -81438.9966 -143921.7752 -40576.7698
58 51484.3647 -97119.5532 -170411.1076 -48208.9727
59 60635.4956 -115107.8 -200588.3413 -56933.2508
60 71002.1797 -135620.714 -234752.7517 -66845.2604
61 82669.1167 -158868.015 -273176.8443 -78034.978
62 95707.7271 -185042.642 -316087.8677 -90581.7111
63 110169.098 -214308.434 -363644.3067 -104547.702
64 126075.115 -246784.425 -415906.3667 -119970.042
65 143407.435 -282525.086 -472799.308 -136850.57
66 162093.905 -321495.751 -534068.3261 -155143.376
67 181991.975 -363542.318 -599223.4924 -174739.479
68 202868.6 -408354.236 -667473.0777 -195448.201
69 224376.049 -455419.627 -737643.3723 -216974.676
70 246022.998 -503971.26 -808082.9024 -238892.885
71 267140.182 -552921.941 -876548.7163 -260613.532
72 286839.841 -600787.732 -940072.1909 -281346.009
73 303968.115 -645597.262 -994801.5879 -300053.629
74 317049.477 -684785.211 -1035818.384 -315401.227
75 324222.241 -715067.955 -1056924.223 -325694.194
76 323164.136 -732299.168 -1050395.189 -328807.927
77 311006.921 -731303.132 -1006700.043 -322106.672
78 284238.993 -705683.417 -914179.0531 -302350.705
79 238595.002 -647604.583 -758680.1815 -265590.807
80 168931.541 -547544.649 -523149.6577 -207049.083
81 69088.1325 -394016.213 -187174.4442 -120985.227
82 -68267.0639 -173254.452 273525.1947 -547.56825
83 -251807.092 131129.367 887653.7785 162391.544
84 -491731.431 538529.488 1688948.694 377416.069
85 -799965.912 1071762.13 2716835.73 655761.874
86 -1190375.84 1757515.28 4017127.052 1010526.51
87 -1678989.93 2626828.63 5642753.399 1456891.19
88 -2284231.44 3715598.37 7654518.307 2012351.91
89 -3027151.33 5065098.89 10121857.06 2696955.26
90 -3931656 6722510 13123577 3533533
91 -5024724 8741434 16748547 4547925
92 -6336587 11182381 21096299 5769188
93 -7900883 14113202 26277482 7229757
94 -9754735 17609417 32414117 8965568
95 -11938752 21754425 39639552 11016100
96 -14496907 26639505 48098034 13424307
97 -17476268 32363573 57943770 16236428
98 -20926524 39032589 69339313 19501603
99 -24899264 46758527 82453122 23271269
100 -29446932 55657779 97456043 27598266
101 -34621388 65848849 114516491 32535570
102 -40471985 77449168 133793992 38134589
103 -47043039 90570815 155430762 44442905
104 -54370590 105314910 179540877 51501338
105 -62478279 121764392 206196559 59340213
106 -71372194 139974854 235411017 67974646
107 -81034475 159963053 267117194 77398683
108 -91415462 181692648 301141677 87578060
109 -102424121 205056687 337172937 98441358
110 -113916473 229856239 374722951 109869266

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy