function Ae = aloc_1dLaplace(Pe,ke)

Ae = ke/(Pe(2)-Pe(1)) * [1,-1;-1,1];