function bk = lokalni_vektor_magnetostatika(Uzlyk,Jk)

Rk = [ Uzlyk(:,2)-Uzlyk(:,1) , Uzlyk(:,3)-Uzlyk(:,1) ];
bk = Jk * [1/3;1/3;1/3] * abs(det(Rk))/2;
