Найти, когда интеграл интерполированной функции равен определенному значению (python)

Question

У меня есть массивы t_array иточек x и y. Назовем M = trapz(dMdt_array, t_array). Я хочу найти, при каком значении t интеграл от dM/dt vs t равен некоторому значению - скажем 0.05*M. В python есть ли хороший способ сделать это?

Я думал что-то вроде F = interp1d(t_array, dMdt_array). Тогда какой-то корень найдет, где интеграл от F равен 0,05 * M. Могу ли я сделать это в python?

Answer 1

Один простой способ - использовать класс CubicSpline. Тогда это CubicSpline(x, y).antiderivative().solve(0.05*M) или около этого.