function priklad_magnetostatika()

mur = 5000;
J = 0.25;

% diskretizacni krok
h = 0.25;

[Uzly,Trojuhelniky,Hranice,Oblasti] = diskretizace_magnetostatika_MKP(h);
[A,b] = MKP_magnetostatika(Uzly,Trojuhelniky,Hranice,Oblasti,mur,J);
u = A\b;
% trisurf(Trojuhelniky',Uzly(1,:),Uzly(2,:),u)
vykresli_mag_indukci(Uzly,Trojuhelniky,u);