Неявное решениеsolve_ivp сокращает количество точек, для которых мне нужно решение ⇐ Python

[Anonymous]

Я пробовал использовать как явные, так и неявные методы, где мой t_eval представляет собой массив размером 55. Явный метод (метод по умолчанию) выводит массив размера (1x55), который я и ожидал.


Там, где я пытаюсь использовать неявное Метод 'BDF'. Я вижу, что в моем решении сохранены только первые два значения:
[img]https://i.sstatic.net /QqoitBnZ.png[/img]

Не уверен, связано ли это с событиями или чем-то еще. Проблема, которую я решаю, кажется сложной, поэтому мне может помочь неявный метод. Я попытался изменить размер шага, используя scipy.integrate.BDF напрямую, но это все равно не сработало. Это пример кода, который я использую для обоих случаев:
sol = solve_ivp(F, [0.1, 1.1], [c_p - h_p[0]], t_eval=h_p)

sol = solve_ivp(F, [0.1, 1.1], [c_p - h_p[0]], t_eval=h_p,
method='BDF', jac=jac)


Подробнее здесь: https://stackoverflow.com/questions/793 ... lution-for