  Numerical integration is encountered in all fields of numerical analysis and
the engineering sciences. By now, various efficient and accurate quadrature
rules are known; for instance, Gauss-type quadrature rules. In many
applications, however, it might be impractical---if not even impossible---to
obtain data to fit known quadrature rules. Often, experimental measurements are
performed at equidistant or even scattered points in space or time. In this
work, we propose stable high order quadrature rules for experimental data,
which can accurately handle general weight functions.
