function priklad1_elektrostatika()

epsilonr = 2;
q = 1;

% diskretizacni krok
h = 0.25;

[Uzly,Trojuhelniky,Hranice,Oblasti] = diskretizace_elektrostatika_MKP(h);
[A,b] = MKP_elektrostatika(Uzly,Trojuhelniky,Hranice,Oblasti,epsilonr,q);
u = A\b;
trisurf(Trojuhelniky',Uzly(1,:),Uzly(2,:),u)
