function priklad_2d()

% diskretizace
h = 0.1;
[P,T,ID,IN] = discretizeRectangle(1,1,h);
n = size(P,2);
m = size(T,2);
mN = size(IN,2);
nD = length(ID);
nonID = setdiff(1:n,ID);

% MKP
ke = ones(1,m); fe = ones(1,m); taue = ones(1,mN);
PD = P(:,ID);
ustriska = 1 + 4*PD(2,:).*(1-PD(2,:));
[A,b] = MKP_2d(P,T,{'alok_gradgrad_2d'},{ke},{'blok_id_2d'},{fe},ID,ustriska,IN,taue);

% reseni
u = zeros(n,1); u(ID) = ustriska; u(nonID) = A\b;
figure; trisurf(T',P(1,:),P(2,:),u')
figure; vykresliGradientReseni(P,T,u);