>> 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)

Osiagnieto wymagana dokladnosc.

x =

6.8083

2.4382

0.8892

2.6914

it =

40

blad =

7.1697e-009

>> h=A*x

h =

30.0000

-0.0000

10.0000

5.0000