>> A=[4 -1 -0.2 2; -1 5 0 -2; 0.2 1 10 -1; 0 -2 -1 4]
A =
4.0000 -1.0000 -0.2000 2.0000
-1.0000 5.0000 0 -2.0000
0.2000 1.0000 10.0000 -1.0000
0 -2.0000 -1.0000 4.0000
>> b=[30;0;10;5]
b =
30
0
10
5
>> [x, it, blad] = gauseidl(A, b, 200, 0.00000001)
Osiagnieto wymagana dokladnosc.
x =
6.8083
2.4382
0.8892
2.6914
it =
14
blad =
6.5911e-009
>> g=A*x
g =
30.0000
-0.0000
10.0000
5.0000
>> [x, it, blad] = jacobi(A, b, 200, 0.00000001)
40
7.1697e-009
>> h=A*x
h =
scrapek.crw