A product quadrature algorithm by Hermite interpolation