我试图在膜上求解 Maxwell Stefan 方程，以获得膜厚度“z”上的瞬态摩尔分数分布。但不知何故，我无法使用 ODE45 对其进行编码，更有可能我无法编写使用 ODE45 解决的系统。如果有人可以帮助我了解主要的语法和功能，那就太好了。我试图解决的方程是

$\frac{dy_{H_2}}{dz}={\frac{1}{C \times D_{H_2,{H_2O}}}}\left[y_{H_2}(N_{H_2}+N_{H_2O})-N_{H_2}\right]$

其中C是浓度，$D_{H_2,H_2O}$是二元扩散系数，N 是摩尔通量，y 是摩尔分数

先感谢您。