如何在Matlab中获得稳态ODE系统的非零符号解？

Question

下面是我为Matlab使用Solve函数编写的代码，以符号形式求解线性ODE的稳态解系统。但是，我只能得到零个解决方案（本文末尾的确切输出）。

syms c1 c2 o1 o2 %variables 
syms a1 a2 b1 b2 h1 h2 g1 g2 %parameters  

eq1 = -(a1 + g2)*c1 + b1*c2 + h1*o1 == 0; %steady state equations
eq2 = -(b1 + g1)*c2 + a1*c1 + h1*o2 == 0;
eq3 = -(a2 + h2)*o1 + g2*c1 + b2*o2 == 0;
eq4 = -(b2 + h1)*o2 + g1*c2 + a2*o1 == 0; 

sol = solve([eq1, eq2, eq3, eq4], [c1, c2, o1, o2]); 
%solve for steady state solutions 

sol1 = sol.c1
sol2 = sol.c2
sol3 = sol.o1
sol4 = sol.o2 
%display solutions 

以下是输出：

>>> source("linsolve.m")
>>> sol1 = (sym) 0
>>> sol2 = (sym) 0
>>> sol3 = (sym) 0
>>> sol4 = (sym) 0

我已经手动制定了该系统的解决方案，但希望对较大的系统（由8个方程和8个变量组成）执行相同的计算，因此希望使用Matlab或其他程序来代替它。

任何帮助将不胜感激！