function b = LoadAssembler1D(x,f)
n = length(x)-1;
b = zeros(n+1,1);
for i = 1:n
    h = x(i+1) - x(i);
    b(i) = b(i) + f(x(i))*h/2;
    b(i+1) = b(i+1) + f(x(i+1))*h/2;
end

>> LoadAssembler1D( [ -10, -4, 8 ], inline('abs(2 * x + 8)', 'x' ))
ans =
    36
     0
   144