Does anybody know if there is a way to use the # operator but with matrices instead of scalars?
I have a situation where 7 steady state ratios that appear in the equations of the model form a linear system of 7 equations. Something like Ax=b, where A and b are functions of deeper parameters: A(theta) and b(theta).
The # operator does not seem to allow to simultaneously solve for the ratios using # x = inv(A) b.
Is there a way to implement this?