functionint_f=CompoundEchelon( f, a, b, m )% input : f : function handler% a : the lower limit of integral% b : the upper limit of integral% m : cut integral area into m peace% output : int_f : the answer of the integralh = (b - a) / m; int_f =0;ifm >=2fori=1:...