function priklad2dmkp()

h = 0.01;

[P,T,ID,IN] = diskretizujObdelnik(1,1,h);
m = size(T,2);
n = size(P,2);
nD = length(ID);
mN = size(IN,2);

f = ones(1,m);
k = ones(1,m);
x = 0:h:1;
uD = sin(pi*x);
tau = zeros(1,mN);

[A,b] = mkp2d(P,T,k,f,ID,uD',IN,tau);

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

trisurf(T',P(1,:),P(2,:),u')

