function prikladMKP1d(h)

l = 1;
x = 0:h:l;
n = length(x)-1;

% nehomogenni material
k1 = 1; k2 = 10;
k = k1*ones(1,n);
k((n/2+1):end) = k2;

f = ones(1,n);
% % nehomogenni zdroj tepla
% f(ceil(n/2):n) = 100;
tau = 0; %1;

ID = [1,n+1];
uD = [1,2];
[A,b] = mkp1d(x,k,f,ID,uD,[],[]);

idx = setdiff(1:(n+1),ID);
u = zeros(n+1,1);
u(ID) = uD;
u(idx) = A\b;

plot(x,u')