ENH: add cumulative_simpson integration to scipy.integrate

Hello everyone!

The cumulative integration and differentiation of numerical samples is a
common operation in the engineering domains. Scipy currently only provides
the integrate.cumulative_trapezoid function for cumulative integration.
However, this is insufficiently accurate for engineering purposes.
The integrate.simpson method is available but only as a complete integral,
not a cumulative integral.

I propose to add a cumulative integration function using Simpson's 1/3 rule
(https://en.wikipedia.org/wiki/Simpson%27s_rule).
CAE software such as ANSYS/LS-DYNA provide simpson's integration as a
default processing option. Adding this to Scipy.integrate will help expand
the utility of the library in the engineering domains.

I have added a PR for this here: https://github.com/scipy/scipy/pull/18151

I would really appreciate your reviews and opinions.

Thank you,
Praveer

_______________________________________________
SciPy-Dev mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://mail.python.org/mailman3/lists/scipy-dev.python.org/
Member address: [email protected]