function prikladMKP1d_verze2(h)

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

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

f = ones(1,n);
alfa = 1; ustriska = 2;
r = 1;
ID = []; uD = [];
IN = []; tau = [];
IR = [1,n+1];
[A,b] = mkp1d_verze2(x,k,f,ID,uD,IN,tau,alfa,r,ustriska,IR);

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

plot(x,u')
axis equal