function priklad3_elektrostatika()

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

% diskretizacni krok
h = 0.5; %0.25;

[Uzly,Trojuhelniky,HraniceMinus,HranicePlus,Oblasti] = diskretizace3_elektrostatika_MKP(h);
[A,b] = MKP3_elektrostatika(Uzly,Trojuhelniky,HraniceMinus,HranicePlus,Oblasti,epsilonr,U);
u = A\b;
h = trisurf(Trojuhelniky',Uzly(1,:),Uzly(2,:),u);
set(h,'LineStyle','none')
