Python — Какой математический метод используется odeint?

Я работаю с scipy.integrate.odeint и хочу понять это лучше. Для этого у меня есть два слегка связанных вопроса:

1. Какой математический метод он использует? Рунге-Кутта? Адамс-Bashforth? я нашел это сайт, но, похоже, это для C ++, но, насколько мне известно, функция python также использует версию C ++ … В ней говорится, что она автоматически переключается между неявным и явным решателем, кто-нибудь знает, как это происходит?

2. Чтобы понять / повторно использовать информацию, которую я хотел бы знать, в какие моменты времени она оценивает функцию и как именно она вычисляет решение ODE, но полный вывод не помогает / я не смог выяснить, как это сделать. Чтобы быть более точным, пример с Рунге-Кутта-Фельбергом: я хочу разные моменты времени, в которые он оценивал f, и веса, которые он использовал для его умножения.

Дополнительная информация (для чего нужна эта информация):
Я хочу повторно использовать эту информацию, чтобы использовать автоматическое разграничение. Поэтому я бы назвал odeint черным ящиком, выяснил все соответствующие шаги и использовал эту информацию для вычисления разницы. dx(T_end)/dx0,

Если вам известен какой-либо другой способ решения моей проблемы, пожалуйста, продолжайте. Также, если для этого подойдет другой решатель од.

PS: я новичок, так что будет лучше разделить этот вопрос на вопросы? То есть отдельно 1. и 2.?