Я использую scikits.odes.odeint для решения относительно большой системы ОДУ (40 переменных, 100 кинетических параметров). Я хочу ускорить процесс, поставляя якобиан решателю. Однако я не знаю, как наилучшим образом определить для этого якобиана и передать его решателю.