function ce = clok1d(x,fe)

h = x(2)-x(1);
ce = (fe*h/2) * [1;1];