The matrix
A=
\begin{pmatrix} A & {\color{red}b}\\
C & D \end{pmatrix}
defines a system y' = A \cdot y
.
Determine the entry {\color{red}b}
such that the system has a steady state.
\color{red} b
=
A*D/C
In general, for an eigenvector v
with corresponding eigenvalue \lambda
, y(t) = e^{\lambda \cdot t}\cdot v
is a solution of the system.
In case \lambda =0
, we get e^{\lambda \cdot t} = e^{0} = 1
, i.e. a corresponding eigenvector
gives a stationary solution.
The eigenvalue zero exists if and only if \det(A) = 0
.
Note that \det(A) = negParens(A) \cdot negParens(D) -
{\color{red}b}
\cdot negParens(C).
Setting this equal to zero and solving for {\color{red}b}
gives:
{\color{red}b} =
fractionReduce(A*D,C).