function [A,b] = mkp(k,f,U,h)

% komentar se pripravuje

n = max(size(k)) + 1;  %  n je pocet uzlu (pocet neznamych), n-1 pocet prvku (usecek)
A = zeros(n,n);
b = zeros(n,1);

for i=1:n-1
    ki = k(i);
    Ai = alok(ki,h);
    A(i:i+1,i:i+1) = A(i:i+1,i:i+1) + Ai;
    
    fi = f(i);
    bi = blok(fi,h);
    b(i:i+1) = b(i:i+1) + bi;
end

% okrajova podminka: -k(1)*u'(1) = u(1)-U
A(n,n) = A(n,n) - 1;
b(n) = b(n) - U;
