function Ae = aloc_1dMass(Pe,qe)

Ae = qe*(Pe(2)-Pe(1)) * [1/3,1/6;1/6,1/3];