function prikladRichardson_1d

h = 0.05;
[A,b] = mkd_1d(h);

% predpodmineni
[L,U] = luinc(A,'0');
M = inv(full(L*U));
%M = diag(1./diag(A)); %inv(full(A));
A = M*A;
b = M*b;

Lambda = Gershgorin(A);
omega = 1/(2*Lambda);

[u,k] = Richardson(A,b,zeros(size(b)),omega,inf,1e-3);
k