Каковы мои варианты решения в системе ODE:
Очевидно, что если A обратим, то я могу решить эквивалентно:
с использованием стандартных методов из scipy.integrate. Но все еще остаются следующие вопросы: как насчет эффективности и что, если A единственное число и не может быть обращено?