Question about the # operator

Hi All,

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?

Thanks,

Mariano