function priklad2_elektrostatika()

epsilonr = 1;
U = 1; % [Volty]

% diskretizacni krok
h = 0.1; %0.25;

[Uzly,Trojuhelniky,Hranice,Oblasti] = diskretizace_elektrostatika_MKP(h);
[A,b] = MKP2_elektrostatika(Uzly,Trojuhelniky,Hranice,Oblasti,epsilonr,U);
u = A\b;
h = trisurf(Trojuhelniky',Uzly(1,:),Uzly(2,:),u);
set(h,'LineStyle','none')
