function Ak = lokalni_matice_magnetostatika(Uzlyk,nuk)

Rk = [ Uzlyk(:,2)-Uzlyk(:,1) , Uzlyk(:,3)-Uzlyk(:,1) ];
Bk = inv(Rk') * [-1 1 0;-1 0 1];
Ak = nuk * Bk' * Bk * abs(det(Rk))/2;
