function prikladRichardson()


% N = [];
% K = [];
% 
% for i=3:30
i = 10;
    h = 1/i;
    [A,b] = priklad(h);
    
    n = length(b);
    [V,D] = eig(full(A));
    
    u = A\b;

    prec = 1e-3;
    Lambda = Gershgorin(A);
    omega = 1/2/Lambda;
    n = length(b);
    u0 = zeros(n,1);

    [u,k] = Richardson(A,b,u0,omega,10000,prec,[V(:,1),V(:,floor(n/2)),V(:,n)]);
    
%     [n,k]
%     N = [ N , n ];
%     K = [ K , k ];
% end
%
% plot(N,K)
