%Jacobijeva iteracija - demonstracija algoritma 2.8.1

n = 4;
A = [-2 1 0 0; 1 -2 1 0; 0 1 -2 1; 0 0 1 -2]
b = [-1 0 0 -1]'
x = [0 0 0 0]'
X = x;

it = 20;
alg281
plot(X)