function prikladMKP1d()

h = 0.05;

l = 1;
k1 = 1;
k2 = 10;
f = 1;
tau = 1;

x = 0:h:1;
m = length(x)-1;
n = m+1;
vek_k = [k1*ones(1,m/2),k2*ones(1,m/2)];
vek_f = f*ones(1,m);
[A,b] = mkp1d(x,vek_k,vek_f,tau);

u = zeros(n,1);
idx = 2:n;
u(idx) = A\b;

plot(x,u')