I think h should equals to (b-a)/n under the index method of example.

Are the f(x_i)'s here considered as quadrature points?