function r = trap_single( f, a, b )

r = 0.5 * (b - a) * ( f(b) + f(a) ) ;