function de = dlok2d(Pe,taue)

delkae = norm(Pe(:,2)-Pe(:,1));
de = (1/2) * taue * delkae * [1;1];