function ce = blok_id_2d(fe,Pe)

Re = [Pe(:,2)-Pe(:,1),Pe(:,3)-Pe(:,1)];
ce = fe * abs(det(Re)) / 6;