function Alok_tau = alok_tau(tau,x,e)

% Alok_tau = alok_tau(tau,x,e)
%             Funkce realizuje lokalni prispevek do matice,
%             ktery odpovida bil. forme:
%              a(v,u) = tau * \int_{x1}^{x2} v'(x)*u'(x) dx

x1 = x(e(1));
x2 = x(e(2));

D = (x2-x1) * tau;
B = (1/(x2-x1)) * [-1 1];

Alok_tau = B' * D * B;
