function be = bloc_1dsource(Pe,fe)

be = fe * (Pe(2)-Pe(1))/2 * [1;1];