quadrature method