Compilations of Algorithms
Compilations of Algorithms
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-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
x3-4x2-1=0
f(x)=x3-4x2-1=0 xnew=x-((x3-4x2-1)/(3x2-8x))
f'(x)=3x2-8x=0
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
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
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
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
f(x)=0.471405*((20x)^(5/3)/(20+2x)^(2/3))-5
Xn=xL-f(xl)*(XR - XL)/(f(XR) - f(XL))