Using Thomas algorithm:
a1=2/3,b1=12/3=4
a2=2/(3−2⋅2/3)=6/5,b2=(17−2⋅4)/(3−2⋅2/3)=27/5
a3=2/(3−2⋅6/5)=10/3,b3=(14−2⋅27/5)/(3−2⋅6/5)=16/3
b4=(7−2⋅16/3)/(3−2⋅10/3)=1
x4=b4=1
x3=b3−a3x4=16/3−10/3=2
x2=b2−a2x3=27/5−6⋅2/5=3
x1=b1−a1x2=4−2⋅3/3=2
Using Jacobi method:
x1(1)=12/3=4,x2(1)=17/3,x3(1)=14/3,x4(1)=7/3
x1(2)=1/3(12−2⋅17/3)=2/9
x2(2)=1/3(17−2⋅4−2⋅14/3)=−1/9
x3(2)=1/3(14−2⋅17/3−2⋅7/3)=−2/3
x4(2)=1/3(7−2⋅14/3)=−7/9
x1(3)=1/3(12+2/9)=110/27
x2(3)=1/3(17−2⋅2/9+2⋅2/3)=161/27
x3(3)=1/3(14+2/9+2⋅7/9)=142/27
x4(3)=1/3(7+2⋅2/3)=25/9
x1(4)=1/3(12−2⋅161/27)=2/81
x2(4)=1/3(17−2⋅110/27−2⋅142/27)=−45/81=−5/9
x3(4)=1/3(14−2⋅161/27−2⋅25/9)=−94/81
x4(4)=1/3(7−2⋅142/27)=−95/81
x1(5)=1/3(12+2⋅5/9)=118/27=4.37
x2(5)=1/3(17−2⋅2/81+2⋅94/81)=1561/243=6.42
x3(5)=1/3(14+2⋅5/9+2⋅95/81)=1414/243=5.82
x4(5)=1/3(7+2⋅94/81)=755/243=3.11
Using Gauss-Seidel method:
x1(1)=12/3=4
x2(1)=1/3(17−2⋅4)=3
x3(1)=1/3(14−2⋅3)=8/3
x4(1)=1/3(7−2⋅8/3)=5/9
x1(2)=1/3(12−2⋅3)=2
x2(2)=1/3(17−2⋅2−2⋅8/3)=23/9
x3(2)=1/3(14−2⋅23/9−2⋅5/9)=90/27=10/3
x4(2)=1/3(7−2⋅10/3)=1/9
x1(3)=1/3(12−2⋅23/9)=62/27
x2(3)=1/3(17−2⋅62/27−2⋅10/3)=155/81
x3(3)=1/3(14−2⋅155/27−2⋅1/9)=90/27=62/81
x4(3)=1/3(7−2⋅62/81)=443/243
x1(4)=1/3(12−2⋅155/81)=662/243
x2(4)=1/3(17−2⋅662/243−2⋅62/81)=2435/729
x3(4)=1/3(14−2⋅2435/729−2⋅443/243)=90/27=1.22
x4(4)=1/3(7−2⋅1.22)=1.52
x1(5)=1/3(12−2⋅2435/729)=1.77
x2(5)=1/3(17−2⋅1.77−2⋅1.22)=3.67
x3(5)=1/3(14−2⋅3.67−2⋅1.52)=90/27=1.21
x4(5)=1/3(7−2⋅1.21)=1.53
Comments