function prikladMKP1d_verze2()

h = 0.05;

l = 1;
k1 = 1;
k2 = 10;
f = 1;
alfa = 1;
r = 0.01;
ustriska = 10;

% diskretizace
x = 0:h:1;
m = length(x)-1;
n = m+1;

% Dirichletovy podminky
ID =[]; % [1,n];
uD = []; %[1,2];

% Neumannovy podminky
IN = [];  tau = [];

% Robinovy podminky
IR = [1,n];

vek_k = [k1*ones(1,m/2),k2*ones(1,m/2)];
vek_f = f*ones(1,m);

[A,b] = mkp1d_verze2(x,vek_k,vek_f,ID,uD,...
    IN,tau,alfa,r,ustriska,IR);

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

plot(x,u')