function bk = lokalni_vektor_elektrostatika(Uzlyk,qk)

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