Мне нужно решить дифференциальное уравнение в форме: dy/dt = H*y
, где H
- матрица n x n
, а y
- вектор-столбец n
.
Я могу переписать уравнениекак y = y_0 * Exp(H*t)
.
Возможны варианты запуска решателя ОДУ в первом уравнении или во второй форме ... и наивно я бы сказал, что экспоненциальная матрица будет более численно интенсивной?
Как масштабируются оба метода с точки зрения вычислительных операций?