function b = bloc_id_L1rect(f,x)

% b = bloc_id_L1rect(f,x)

area = (x(1,2)-x(1,1))*(x(2,4)-x(2,1));

b = (area/4) * f * ones(4,1);