add quadrature for path and alleventlist

[tlienart]

[migrated from private repo]

• For polynomials, we can do exact quadrature which can be directly implemented along a piecewise linear path
• For general functions, we can use a simple GL or GK quadrature with a number of points depending on the size of the interval (proportional) with a number of function evaluation as a budget
• If feeling fancy, we could do adaptive quadrature but it's kind of overkill

Note that this is all well and easy for reduced dimensionality, as usual.

[tlienart]

Working on this with the global BPS. The test should be that with the samplepath function and sufficiently many sampling points, the quadrature and the MC estimator should be approximately equal.

[tlienart]

this is now done in e3f3727 and tested in https://github.com/alan-turing-institute/PDMP.jl/blob/master/test/path_test.jl#L17_L27

It follows the trivial extension of formula below to any polynomial