function ce = bloc_2dsource(Pe,fe)

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