function q = composite_simpson13 ( m, a, b, f )

  breaks = linspace ( a, b, m + 1 );

  q = 0.0;
  for i = 1 : m
    ai = breaks(i);
    bi = breaks(i+1);
    q = q + simpson13 ( ai, bi, f );
  end

  return
end
