我正在尝试以数值方式求解以下微分方程组。有哪些可用的有限差分方法和 matlab 求解器来求解这样的系统？也欢迎其他解决问题的方法。前两个方程是准稳态的，其余方程在其中具有时间导数。我正在尝试使用 matlab 的 pdepe 求解器，但是由于即使对于具有时间导数的 ODE 也必须指定边界条件，所以我认为它不是一个可行的选择。

$$0= \rho \frac{\partial}{\partial z}\left(D \frac{\partial c_1}{\partial z} \right) - \it{source_1}\,$$

$$0= \rho \frac{\partial}{\partial z}\left(D \frac{\partial c_2}{\partial z} \right) - \it{source_2}\,$$

$$\frac{\partial c_3}{\partial t} = source_3$$

$$\frac{\partial c_4}{\partial t} = source_4$$

$$D = A_1 c_1 + A_2 c_2 + \dots + A_n c_n \,, \rho = 1 (say)\,.$$

任何特定方程（或浓度）的源项都包含对其他浓度的依赖性。